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DISCLAIMER 


The  Mass  Flux  Toolkit  is  made  available  on  an  as-is  basis  without  guarantee  or  warranty  of  any 
kind,  expressed  or  implied.  The  United  States  Government,  Groundwater  Services,  Inc.,  the 
authors  and  reviewers  accept  no  liability  resulting  from  the  use  of  the  Mass  Flux  Toolkit  or  its 
documentation.  Implementation  of  the  Mass  Flux  Toolkit  and  interpretation  of  the  predictions  of 
the  model  are  the  sole  responsibility  of  the  user. 
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INTRODUCTION 


Traditional  site  investigations  have  focused  on  collecting  both  chemical  concentrations  and 
subsurface  flow  characteristics.  These  data  have  typically  been  analyzed  separately,  such  as 
when  chemical  concentrations  are  compared  to  regulatory  standards.  When  the  two  types  of 
data  are  combined,  it  is  typically  for  use  in  solute  transport  models. 

Recently  there  has  been  increased  emphasis  on  combining  chemical  and  flow  data  into  a  single 
metric  to  provide  estimates  of  contaminant  mass  flux.  This  flux  is  the  rate  that  a  contaminant 
passes  through  a  defined  cross  section  in  units  of  mass  per  time.  Although  the  strict  definition  of 
mass  flux  is  mass  discharge  per  unit  area  (in  units  of  mass  per  time  per  area),  many  researchers 
and  practitioners  use  the  term  mass  flux  to  denote  the  rate  of  “mass  discharge”  or  “total  mass 
discharge".  This  document  will  use  the  term  “mass  flux”  to  represent  data  in  units  of  mass  per 
time  passing  by  a  plume  transect  perpendicular  to  groundwater  flow. 

Regulatory  acceptance  of  this  approach,  and  the  number  of  sites  where  mass  flux  has  been 
used,  are  rapidly  increasing.  More  research  and  industry  groups  (U.S.  EPA,  1998;  Devlin  et  al., 
2002;  Connecticut,  2002;  Rao  et  al.,  2002;  Newell  et  al.,  2003;  ASTM,  1998,  Kavanaugh  et  al., 
2003)  suggest  that  a  mass  flux  approach  may  provide  key  insights  about  how  groundwater 
plumes  behave  and  how  they  should  be  managed.  For  example: 

The  USEPA’s  Scientific  Advisory  Board  (U.S.  EPA,  2001)  recommended  that: 

“In  summary,  measurements  of  mass  flux  of  the  contaminants  and  footprint  parameters  -  not 
just  concentrations  -  are  necessary  to  document  cause-and-effect  and  to  assess  long-term 
sustainability/permanence.  Site-characterization  and  monitoring  plans  should  be  proactively 
designed  to  accommodate  mass-flux  estimates.  ” 

Pankow  and  Cherry  (1996)  state  that: 

"Therefore,  the  ultimate  impact  of  plumes  emanating  from  solvent  DNAPL  source  zones  can 
be  evaluated  in  terms  of  impact  of  relatively  small  annual  mass  fluxes  to  the  receptor  such  as 
water-supply  wells  or  surface  waters.  In  some  cases,  the  fluxes  present  significant  risk  to 
human  health  and/or  the  environment,  and  extensive  remedial  action  is  warranted.  In  other 
cases,  the  fluxes  are  insignificant,  and  remedial  action  would  provide  little  or  no  actual 
environmental  risk  reduction. " 

Despite  the  interest  in  using  the  mass  flux  approach,  calculating  mass  flux  from  field  data  is  not 
easy.  There  are  several  methods  to  choose  from,  each  with  its  own  strengths  and  limitations. 
The  most  commonly  used  method,  calculating  mass  flux  from  transect  data,  remains  a  difficult 
and  laborious  process.  The  concentration  data  is  multiplied  against  the  corresponding  flow  area 
data  to  give  flux  for  a  specific  cross-sectional  area  of  the  plume.  The  results  from  all  the  flow 
areas  are  summed  together  to  obtain  a  mass  flux  across  the  transect.  However,  few  groundwater 
monitoring  networks  provide  data  that  offer  a  straightforward  calculation  of  mass  flux  (Figure  1). 
Rather,  a  transect  drawn  across  a  plume  often  intersects  monitoring  locations  at  different  lateral 
spacings,  and  each  well  may  have  a  unique  screened  interval.  Hydraulic  characteristics 
determined  by  aquifer  tests  are  rarely  defined  at  all  wells,  and  can  vary  from  one  well  to  the  next. 
Thus,  a  consistent  evaluation  of  mass  flux  may  require  a  significant  amount  of  analysis  time  and 
effort. 

Other  methods,  such  as  extraction  well  capture  and  in-situ  flux  meters  are  less  common.  There 
is  a  need  for  increased  tech  transfer  information  on  how  to  apply  these  new  technologies. 

There  are  several  methods  that  use  mass  flux  data  for  site  assessment  and  plume  management. 
Einarson  and  MacKay  (2001a)  present  simple  dilution  calculations  based  on  mass  flux  estimates 
for  plumes  and  pumping  rates  from  wells  that  capture  this  mass  flux.  The  American  Petroleum 
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Institute’s  (API)  “Groundwater  Remediation  Strategies  Tool”  (Newell  et  al.,  2003)  provide  a  useful 
framework  for  using  mass  flux  for  making  site  management  decisions,  but  the  key  worksheets 
from  this  product  will  only  be  available  in  a  paper  form. 

Finally,  there  is  a  concern  that  mass  flux  estimates  can  have  significant  uncertainty  (Kavanaugh 
et  al.,  2003).  For  example,  Fraser  et  at.,  (2005)  evaluated  mass  flux  vs.  sampling  density  for  a 
naphthalene  plume  at  the  Borden  research  aquifer.  When  the  sampling  grid  density  was  reduced 
from  1.7  points  per  meter  squared  to  0.7  points  per  meter  squared,  the  range  (as  a  standard 
deviation)  in  mass  discharge  increased  to  more  than  50%.  Guilbeault  et  al.  (2005)  showed  that 
75%  of  the  mass  flux  occurred  within  5%  to  10%  of  the  transect  cross  section  area  for  three 
plumes  in  Ontario,  New  Hampshire  and  Florida,  and  that  a  spacing  no  larger  than  15  to  30  cm 
was  needed  at  some  locations  to  identify  high  concentration  zones. 

To  help  site  managers  and  site  consultants  estimate  mass  flux  and  understand  the  uncertainty  in 
those  estimates,  the  Environmental  Security  Technology  Certification  Program  (ESTCP)  of  the 
U.S.  Department  of  Defense  (DoD)  has  funded  the  development  of  a  computerized  Mass  Flux 
Toolkit  (Toolkit). 

The  Mass  Flux  Toolkit  is  an  easy-to-use,  free  software  tool  that  gives  site  personnel  the  capability 
to  compare  different  mass  flux  approaches,  calculate  mass  flux  from  transect  data,  and  apply 
mass  flux  to  manage  groundwater  plumes.  In  this  software,  the  term  mass  flux  is  used  to 
describe  the  mass  discharge  rate  in  a  groundwater  plume  in  units  of  mass  per  time  passing 
across  a  plume  transect. 

The  software,  programmed  in  the  Microsoft®  Excel  spreadsheet  environment,  provides  the 
following  tools  to  calculate  and  evaluate  mass  flux  of  contaminants  in  groundwater: 

1)  Calculation  of  Mass  Flux 

A  module  to  calculate  the  total  mass  flux  across  one  or  more  transects  of  a  plume  and 
plot  mass  flux  vs.  distance  to  show  the  effect  of  remediation/impact  of  natural  attenuation 
processes.  Three  types  of  uncertainty  analysis  are  included:  uncertainty  range  due  to 
interpolation;  uncertainty  due  to  the  variability  in  the  input  data  using  a  Monte-Carlo  like 
analysis;  and  an  uncertainty  analysis  that  shows  how  dependent  the  mass  flux  estimate 
is  on  data  from  each  monitoring  point. 

2)  Impact  of  Mass  Flux 

A  module  allowing  users  to  perform  critical  dilution  calculations  for  plumes  approaching 
production  wells  or  streams.  An  additional  feature  calculates  the  capture  zone  of  the 
supply  well,  directing  the  user  to  alter  the  transect  dimensions  if  the  transect  does  not 
encompass  the  capture  zone. 

3)  Learn  About  Mass  Flux 

A  module  that  provides  a  review  of  theory  and  methods  of  estimating  mass  flux. 

a.  Overview  of  Mass  Flux  Calculations:  a  module  that  provides  information  about 
different  mass  flux  calculation  methods; 

b.  Uncertainties  in  Mass  Flux  Estimates:  a  module  that  provides  information  about  the 
main  sources  of  uncertainty  associated  with  mass  flux  estimates; 

c.  Emerging  Mass  Flux  Resources:  a  review  of  the  emerging  mass  flux  methods; 

d.  Extraction  Well  Capture  Summary:  a  review  of  theory  and  methods  of  estimating 
mass  flux  from  a  pumping  well  or  wells  that  fully  capture  a  contaminant  plume; 

e.  In-Situ  Flux  Meter  Summary:  a  detailed  explanation  of  this  emerging  mass  flux 
method; 

f.  Computer  Model  Resources:  an  outline  of  how  commonly  available,  public  domain 
models  such  as  BIOSCREEN  and  BIOCHLOR  can  be  used  to  obtain  mass  flux 
estimates  at  several  points  in  a  plume; 
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g.  API  Groundwater  Remediation  Strategies  Tool:  an  electronic  version  of  the  American 
Petroleum  Institute’s  (API)  Groundwater  Remediation  Strategies  Tool  to  help  users 
apply  mass  flux  to  groundwater  remediation  decision  making. 


The  Mass  Flux  Toolkit  was  developed  for  the  Environmental  Security  Technology  Certification 
Program  (ESTCP)  by  Groundwater  Services,  Inc.,  Houston,  Texas. 
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INTENDED  USES  FOR  MASS  FLUX  TOOLKIT  AND 
LIMITATIONS 


The  Mass  Flux  Toolkit  attempts  to  assist  site  managers  and  site  consultants  perform  mass  flux 
calculations  more  quickly  and  cheaply,  permitting  their  inclusion  in  more  evaluations  of 
groundwater  plumes.  With  mass  flux  data,  the  progress  of  natural  attenuation  and  remediation 
systems  can  be  demonstrated  more  vividly  and  directly  to  regulators,  making  the  entire 
remediation  process  clear,  effective,  and  efficient. 

The  Mass  Flux  Toolkit  is  intended  to  be  used  in  four  ways: 

1.  As  a  tool  for  calculating  mass  flux  and  performing  dilution  calculations.  The 

Toolkit  can  be  used  to  help  estimate  the  mass  flux  using  the  transect  method. 
Although  other  mass  flux  estimation  techniques  are  emerging,  the  transect  method  is 
an  established  technology.  The  Toolkit  brings  together  key  technical  resources,  easy- 
to  use  calculation  worksheets,  and  case  studies  together  in  one  easy-to-access 
platform. 

a.  In  addition,  the  Toolkit  provides  three  methods  for  analyzing  uncertainty  in  the 
total  mass  flux  estimates  derived  from  the  transect  method.  One  option  utilizes 
the  Monte  Carlo  type  approach  to  analyze  uncertainty  in  the  actual  concentration, 
hydraulic  conductivity,  and  gradient  measurements.  With  this  tool,  groundwater 
practitioners  can  estimate  the  accuracy  of  the  hydrologic  measurements  that  are 
being  used  for  the  mass  flux  calculation. 

The  second  option  provides  a  tool  for  estimating  the  contribution  of  each 
individual  observation  to  the  total  mass  flux.  If  a  single  monitoring  point 
represents  a  high  percentage  of  the  total  mass  flux,  then  the  uncertainty  in  the 
calculation  is  high  and  additional  monitoring  points  should  be  added  to  reduce 
the  uncertainty. 

The  third  method  shows  the  uncertainty  involved  in  the  interpolation  scheme  that 
is  used  to  calculate  mass  flux. 

b.  The  Toolkit  can  also  be  used  to  perform  critical  dilution  calculations  for  plumes 
approaching  production  wells  or  streams.  An  additional  feature  calculates  the 
capture  zone  of  the  supply  well  and  compares  it  to  the  transect  used  to  calculate 
the  mass  flux,  directing  the  user  to  alter  the  transect  dimensions  if  the  transect 
does  not  encompass  the  capture  zone. 

2.  As  a  data  management  system  for  storing  and  interpreting  mass  flux  data. 

Many  sites  will  continue  to  be  monitored  for  a  long  period  of  time.  Data  from  this 
monitoring  program  can  be  stored  in  MAROS  (Monitoring  and  Remediation 
Optimization  System,  a  data  management  system,  http://www.gsi- 
net.com/software/MAROS.htm)  and  then  used  to  update  remediation  timeframe  data. 
This  database  function  can  be  used  to  calculate,  store,  and  compare  mass  fluxes  for 
different  time  periods. 

3.  As  a  tool  for  identifying  what  information  is  needed  to  reduce  the  uncertainty  in 
mass  flux  estimates.  Applying  the  Toolkit  to  a  site  can  help  users  identify  what 
types  of  data  are  required  to  improve  the  mass  flux  estimates  and  reduce  the 
uncertainty  in  the  estimate. 
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4.  As  a  tool  for  learning  about  and  comparing  various  methods  of  calculating 
mass  flux.  The  Toolkit  reviews  emerging  mass  flux  methodologies  with  an  emphasis 
on  data  requirements,  applications,  advantages,  and  disadvantages  of  each 
approach. 

In  summary,  the  software  will  give  site  personnel  the  capability  to  compare  different  mass  flux 
approaches,  calculate  mass  flux  from  transect  data,  and  apply  mass  flux  to  manage  groundwater 
plumes.  With  this  tool,  site  personnel  will  be  able  to  perform  mass  flux  calculations  more  quickly 
and  cheaply,  permitting  their  inclusion  in  more  evaluations  of  groundwater  plumes.  With  the 
mass  flux  data,  the  progress  of  natural  attenuation  and  remediation  systems  can  be 
demonstrated  more  vividly  and  directly  to  regulators,  making  the  entire  remediation  process  clear, 
effective,  and  efficient. 

The  Toolkit  has  the  following  assumptions  and  limitations: 

•  The  Toolkit  assumes  the  user  is  familiar  with  basic  groundwater  transport  and  mass 
balance  concepts. 

•  The  Toolkit  assumes  uniform  flow  and  concentration  throughout  each  grid  cell. 

•  The  Toolkit  assumes  that  the  plume  top,  bottom,  and  sides  are  known. 

•  The  Toolkit  assumes  that  the  sampling  network  is  dense  enough  to  characterize  the 
structure  of  the  plume.  Fraser  et  al.,  (2005)  and  Guilbeault  et  al.  (2005)  show  that 
monitoring  networks  need  to  have  a  very  extensive  system  of  monitoring  points  to 
capture  the  heterogeneity  of  the  plume  in  the  Y-Z  plane.  This  means  that  mass  flux 
estimates  can  have  significant  uncertainty. 

•  To  run  the  Monte  Carlo  analysis,  users  need  to  estimate  what  type  of  statistical 
distribution  best  fits  the  input  data  and  what  values  best  describe  the  distribution.  In 
many  cases  there  will  not  be  data  to  make  these  estimates,  so  the  user  may  have  to  rely 
on  scientific/engineering  judgment  to  use  the  Monte  Carlo  analysis. 

•  The  Monte  Carlo  analysis  cannot  account  for  plume  data  that  are  not  part  of  the 
monitoring  system.  Actual  mass  flux  values  can  be  outside  the  range  of  the  reported 
range  of  mass  flux  values  from  the  Monte  Carlo  analysis  (for  example,  if  new  data  shows 
high  concentration  zones  that  were  not  captured  by  the  original  monitoring  network). 

•  The  capture  zone  module  does  not  adjust  the  mass  flux  based  on  the  degree  of  capture, 
but  only  provides  information  on  the  estimated  size  of  the  capture  zone. 
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MASS  FLUX  TOOLKIT  MODELS 


The  Mass  Flux  Toolkit  software  has  models  for  calculating  mass  flux  and  determining  its  impact 
on  plumes  approaching  a  stream  or  supply  well. 


Mass  Flux  Transect  Calculator 

The  Toolkit  calculates  the  total  mass  flux  across  one  or  more  transects  of  a  plume,  using  the 
transect  approach,  and  plots  mass  flux  vs.  distance  to  show  the  effect  of  remediation/impact  of 
natural  attenuation  processes.  The  information  presented  in  this  section  is  obtained  directly  from 
the  API’s  “Groundwater  Remediation  Strategies  Tool”  (Newell  et  al.,  2003): 

Estimation  of  mass  flux  across  one  or  more  transects  through  an  affected  groundwater  plume 
involves  the  following  principal  steps: 

1.  Characterize  Plume  Concentrations:  For  each  selected  plume  transect,  sufficient 
groundwater  sampling  points  must  be  available  to  define  i)  the  full  width  and  thickness  of  the 
plume  and  ii)  the  distribution  of  contaminant  concentrations  within  the  plume.  Either  single- 
level  or  multilevel  groundwater  monitoring  points  may  be  used  for  this  purpose.  Multilevel 
monitoring  points  can  provide  a  more  detailed  three-dimensional  characterization  of 
contaminant  concentrations  in  groundwater.  However,  single-level  groundwater  monitoring 
networks,  while  less  accurate  than  multi-level  networks,  can  still  provide  sufficient  accuracy  to 
support  a  mass  flux  analysis  at  many  sites. 

(NOTE:  an  alternative  method  is  to  use  the  concentrations  from  contour  lines  as 

measurement  points.  The  locations  where  the  transect  intersects  contour  lines  from  plume 
maps  can  be  used  to  construct  flow  areas  for  the  mass  flux  calculations). 

2.  Characterize  Groundwater  Flow:  To  characterize  the  specific  discharge  (q)  across  each 
plume  transect,  representative  measurements  are  required  for  both  the  hydraulic  flow 
gradient  (i)  and  the  hydraulic  conductivity  (K)  of  the  flow  system  (where  q  =  K  x  i).  The 
groundwater  flow  direction  and  hydraulic  gradient  for  each  segment  of  a  transect  line  can  be 
determined  from  a  potentiometric  surface  contour  map  based  on  static  water  level 
measurements  of  available  sampling  points.  Representative  measurements  of  the  hydraulic 
conductivity  of  the  groundwater-bearing  unit  should  be  obtained  at  one  or  more  locations, 
using  appropriate  slug  test  or  pumping  test  methods. 

3.  Select  Plume  Transects:  To  characterize  mass  flux,  transects  should  be  located  at  points 
where  sufficient  data  are  available  to  define  affected  groundwater  concentrations  and  specific 
discharge,  as  defined  in  Steps  1  and  2  above.  For  two-dimensional  data  (i.e. ,  from  single¬ 
screen  monitoring  wells),  the  transect  will  represent  a  line  extending  across  the  full  width  of 
the  plume,  perpendicular  to  the  direction  of  groundwater  flow  (see  Figures  1  and  2).  For 
uniform  flow  fields,  this  transect  will  be  a  straight  line,  but,  for  converging  or  diverging 
flowlines,  the  transect  will  be  curvilinear  in  shape.  For  three-dimensional  data  (i.e.,  from 
multilevel  monitoring  wells),  the  transect  line  will  represent  a  vertical  plane  through  the 
groundwater  plume,  positioned  perpendicular  to  groundwater  flow  (see  Figures  3  and  4). 
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Figure  1.  Example  Transects  through  2-D  Plume  Contour  Map,  Dover  AFB,  Delaware 
(Adapted  from:  Einarson,  2001) 


Transect  1:  Intersection  with  Contour  Lines 
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Md  =  Z  C,A,qj  ~  630  grams/dav 


Figure  2.  Concentration  Profile  for  2-D  Transect  No.  1  Based  on  Concentration  Contours  Shown 
in  Figure  2  (Adapted  from:  Einarson,  2001) 


Transect  1  Transect  2 


Figure  3.  Example  Transects  through  3-D  Plume  Delineation 
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Figure  4.  MTBE  Concentration  Profile  for  3-D  Transect  No.  1 
(Adapted  from:  Einarson,  2001) 


4.  Subdivide  Transects  Into  Subareas:  Each  transect  should  be  divided  into  subareas. 
Typically,  each  subarea  represents  a  different  concentration  value.  Two  methods  are 
commonly  applied.  For  Method  1,  subareas  are  divided  to  represent  the  area  between 
concentration  isopleths  on  a  contour  map  of  the  plume.  The  concentration  is  assumed  to  be 
the  geometric  mean  of  the  two  contour  values.  For  Method  2,  sufficient  monitoring  points  are 
located  directly  on  the  transect  to  construct  transect  subareas.  The  dividing  line  between 
subareas  is  typically  halfway  between  the  measurement  points.  In  some  cases,  a 
combination  of  Method  1  and  Method  2  can  be  applied.  For  three-dimensional  transects,  the 
transect  plume  should  be  subdivided  into  polygons  bounded  by  contour  data  (Method  1)  or 
centered  on  available  measuring  points  (Method  2)  (See  Figures  3  and  4  for  an  example  of 
Method  2). 

5.  Calculate  Cumulative  Mass  Flux  Across  Transect:  The  total  contaminant  mass  flux 
across  the  transect  is  calculated  as  follows: 


w 


Zi=n 

,„C,gACF 


where: 

w  =  total  mass  flux  from  source  zone  (g/day)  (also  called  mass  discharge) 

C,  =  concentration  of  constituent  at  flow  area  in  transect  (mg/L) 

q,  =  specific  discharge  (also  called  Darcy  velocity)  through  flow  area  associated  with  an 
individual  constituent  measurement  i  (cm/sec),  q,  can  be  calculated  using: 

cji  =  K-i 


where: 

K  =  hydraulic  conductivity  (cm/sec) 
i  =  hydraulic  gradient  (cm/cm) 

A,  =Flow  area  associated  with  an  individual  constituent  measurement  (ft2) 

CF=  conversion  factor  =  80.3  (ft/cm)/(sec/day)(L/ft3)(g/mg) 

Guidelines  for  selecting  key  input  parameters  for  the  model  are  outlined  in  Mass  Flux  Transect 
Data  Entry.  For  help  on  results  see  Mass  Flux  Transect  Results. 
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Uncertainty  Analysis 

Uncertainty  in  mass  flux  estimates  is  a  key  issue  in  using  mass  flux  as  a  metric.  There  are  three 
main  sources  of  uncertainty  in  a  mass  flux  estimated  from  transect  data: 

Type  1:  Uncertainty  in  the  actual  concentration,  hydraulic  conductivity,  and  gradient 
measurements.  Uncertainty  in  the  water  quality  data  (concentration)  can  be  estimated  by 
evaluating  temporal  trends  over  time  and  evaluating  the  amount  of  data  scatter.  Groundwater 
practitioners  can  estimate  the  accuracy  of  the  hydrologic  measurements  (hydraulic  conductivity 
and  hydraulic  gradient)  that  are  being  used  for  the  mass  flux  calculation. 

Type  2:  Uncertainty  in  the  interpolation  scheme.  Different  interpolation  schemes  will  result  in 
different  mass  flux  estimates.  Some  interpolation  schemes,  such  as  kriging,  provide  local 
estimates  of  uncertainty,  however,  kriging  is  not  be  used  in  the  Toolkit.  The  effect  of  the  different 
interpolation  schemes  can  be  addressed  by  performing  the  mass  flux  calculation  with  the  different 
schemes  included  in  the  Toolkit. 

Type  3:  Uncertainty  associated  with  unmeasured  values.  This  type  of  uncertainty  is  related  to 
Type  2  uncertainty.  However,  the  uncertainty  associated  with  areas  of  high  mass  flux  that  are 
missed  by  the  monitoring  scheme  is  difficult  to  assess. 

The  Toolkit  provides  options  for  characterizing  Type  1  uncertainty  quantitatively  using  two 
features:  Monte  Carlo  analysis  and  crossvalidation.  Guidelines  for  selecting  key  input  parameters 
for  the  model  are  outlined  in  Uncertainty  Analysis.  Type  2  uncertainty  is  handled  by  calculating 
mass  flux  with  all  possible  interpolation  schemes,  and  showing  the  range  in  the  resulting  mass 
flux  estimates. 

The  Toolkit  calculates  Type  2  uncertainty  in  mass  flux  for  the  Nearest  Neighbor,  Linear 
Interpolation,  and  Log  Transformation  for  any  grid  that  requires  interpolation  (except  User 
Interpolated)  and  plots  the  minimum  and  maximum  values  as  error  bars  on  the  Mass  Flux 
Summary  sheet. 

Simple  Uncertainty  Analysis  (Spatial  Uncertainty  Analysis) 

The  Toolkit  applies  crossvalidation  to  evaluate  the  uncertainty  in  the  sampling  point  location. 
Crossvalidation  is  a  method  of  removing  one  observation  at  a  time,  estimating  a  value  for  the 
removed  observation  using  the  remaining  observations,  and  calculating  the  output  (Olea,  1999). 
The  removed  value  is  then  replaced  and  the  next  observation  removed.  This  process  is 
continued  until  all  the  observations  have  been  removed  one  at  a  time.  The  Toolkit  will  remove 
each  observed  value  of  concentration,  hydraulic  conductivity,  and  gradient  one  at  a  time  during 
crossvalidation. 

Advanced  Uncertainty  Analysis  (Input  Uncertainty  Analysis) 

Monte  Carlo  analysis  is  a  method  of  analyzing  and  quantifying  uncertainties  in  model  outputs  due 
to  the  uncertainties  in  the  input  parameters  (Rong  et  a!.,  1998).  Monte  Carlo  analysis  refers  to  a 
computer  based  system  that  uses  random  numbers  from  a  probability  distribution  to  obtain  an 
approximation  for  the  parameter  of  interest  (U.S.  EPA,  1997;  Bergin  and  Milford,  2000). 

In  the  standard  Monte  Carlo  approach,  simple  random  sampling,  a  large  number  of  runs  (typically 
100  to  1000)  are  required  to  obtain  a  meaningful  probability  distribution  for  the  parameter.  For 
each  run  of  the  standard  approach,  a  random  number  is  generated  for  every  value  of 
concentration,  hydraulic  conductivity,  and  hydraulic  gradient  entered  by  the  user.  This  set  of 
random  inputs  is  then  used  to  create  and  interpolate  the  respective  grids  and  mass  flux 
calculated  from  those  grids.  Repeating  this  procedure  a  large  number  times  yields  a  probability 
distribution  of  the  mass  flux  from  which  statistical  characteristics  such  as  mean,  percentile,  and 


MASS  FLUX  TOOLKIT 
▼  USER’S  MANUAL  ▼ 


9 


MODEL  TYPES 


variance  can  be  obtained.  However,  the  standard  Monte  Carlo  technique  tends  to  be  time 
intensive  because  of  the  large  number  of  iterations  requires. 

Latin  hypercube  sampling  is  a  type  of  stratified  Monte  Carlo  sampling.  It  can  be  viewed  as  a 
method  for  controlling  the  upper  and  lower  ends  of  the  distributions  used  in  the  analysis  such  that 
all  portions  of  the  distribution  are  well  represented  (U.S.  EPA,  1997;  Iman  and  Helton,  1988;  Gwo 
et  al.,  1996;  Helton  and  Davis,  2002).  In  LHS,  the  given  range  of  each  variable  in  Xi,...,Xn  is 
divided  into  j  number  of  equal  non-overlapping  probability  intervals.  A  random  number  with 
respect  to  the  probability  density  is  then  obtained  for  X )  for  each  j  interval.  The  j  values  of  X-i  that 
are  obtained  are  then  randomly  paired  with  the  j  values  obtained  for  X2  and  so  on  until  all  the  j 
values  in  Xn  are  paired  (Iman  and  Helton,  1988).  LHS  ensures  that  a  minimal  number  of 
probability  intervals  are  required  to  represent  the  full  range  of  the  distribution  (Bergin  etal.,  1999). 
Additionally,  LHS  requires  fewer  simulations  than  simple  random  sampling  to  achieve  the  same 
level  of  precision  (U.S.  EPA,  1997;  Helton  and  Davis,  2002),  hence  a  fraction  of  the  time  involved 
with  simple  random  sampling. 

LHS  sampling  was  retained  as  the  methodology  of  choice  in  the  Mass  Flux  Toolkit  for 
computational  efficiency. 


Impact  of  Mass  Flux 

The  Mass  Flux  Toolkit  software  has  models  for  determining  the  impact  of  mass  flux  for  plumes 
approaching  a  stream  or  supply  well. 

Receptor  impact  Worksheet  -  Wells 

The  Toolkit  determines  the  maximum  concentration  in  water  extracted  from  a  supply  well  using 
the  relationship  described  by  Einarson  and  Mackay  (2001a): 

Csw  =  Md  /  Qsw 

where: 

Csw  =  the  maximum  concentration  of  a  contaminant  in  water  extracted  from  a  supply  well 
(mass/volume) 

Qsw  =  is  the  pumping  rate  from  the  supply  well  (volume/time),  and 

Md  =  the  total  mass  flux  or  mass  discharge  of  the  contaminant  from  the  well  (mass/time), 
which  can  be  predicted  before  the  plume  reaches  the  well  using  the  entire  mass  flux 
through  upgradient  transects. 

Additionally,  the  Toolkit  calculates  the  capture  zone  of  the  supply  well  using  the  relationship 
described  by  Javandel  and  Tsang  (1986),  directing  the  user  to  alter  the  transect  dimensions  if  the 
transect  does  not  encompass  the  capture  zone,  thus  ensuring  accuracy  in  the  model. 

Receptor  Impact  Worksheet  -  Streams 

The  Toolkit  applies  the  dilution  calculation  used  for  water  supply  wells  to  plumes  that  are 
discharging  to  streams. 

Guidelines  for  selecting  key  input  parameters  for  the  model  are  outlined  in  Receptor  Impact  Data 
Entry.  For  help  on  results  see  Impact  of  Receptor  -  Wells  Results  and  Impact  of  Receptor  - 
Streams  Results. 
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MASS  FLUX  TRANSECT  CALCULATOR 
MAIN  SCREEN 
INPUT  DATA  AND  GRID 

GRID  COMPLETION 

MASS  FLUX  RESULTS 

FINAL  CONCENTRATION  GRID 

MASS  FLUX  SUMMARYUNCERTAINTY  ANALYSIS 


Three  important  considerations  regarding  data  input  are: 

1 )  To  see  the  example  data  set  in  the  input  screen  of  the  software,  click  on  the 
Paste  Example  button  on  the  lower  right  portion  of  the  input  screen. 

2)  Because  the  Mass  Flux  Toolkit  is  based  on  the  Excel  spreadsheet,  you  have  to  click 
outside  of  the  cell  where  you  just  entered  data  or  hit  Return  before  any  of  the 
buttons  will  work.  Additionally,  REMOVING  rows  or  columns  from  input  screens  may 
cause  the  program  to  crash. 

3)  Parameters  used  in  the  model  are  to  be  entered  directly  into  the  white  cells. 


NOTE:  Although  literature  values  are  provided,  it  is  strongly  recommended  that  the  user  employ 
measured  hydrogeological  and  source  characteristic  values  whenever  possible.  If  literature 
values  are  used  and  there  is  uncertainty  in  the  value  chosen,  sensitivity  analyses  should  be 
conducted  to  determine  the  effects  of  the  uncertainty  on  model  predictions. 


MASS  FLUX  TOOLKIT 
▼  USER’S  MANUAL  ▼ 


11 


TRANSECT  MAIN  SCREEN  DATA  ENTRY 


Mass  Flux  Transect  Calculator  Data  Entry 


PARAMETER 

UNIT  SYSTEM 

Description 

Unit  system  to  perform  mass  flux  calculations  in. 

How  to  Enter  Data 

Choose  the  appropriate  radio  button. 

PARAMETER 

CONCENTRATION  UNITS 

Units 

mg/L  or  ug/L 

Description 

Concentration  measurement  units. 

How  to  Enter  Data 

Choose  from  drop  down  list. 

PARAMETER 

ELEVATION  UNITS 

Units 

Below  ground  surface  (bgs)  or  Mean  Sea  Level  (MSL) 

Description 

Elevation  measurement  units. 

How  to  Enter  Data 

Choose  from  drop  down  list.  If  MSL  is  selected  then  a  value  for  the  Ground 
Surface  Elevation  must  be  entered. 

PARAMETER 

CONTINUE  INPUT  OF  TRANSECT  DATA 

Description 

Proceeds  to  the  next  step  in  entering  data  for  mass  flux  calculations. 

PARAMETER 

RETURN  TO  MAIN  SCREEN 

Description 

Returns  to  the  Mass  Flux  Toolkit  Main  Screen. 

PARAMETER 

NEW  SITE 

Description 

Clears  ALL  data  in  the  Toolkit  memory  banks.  Use  this  button  to  start  a  new 
project. 

PARAMETER 

PASTE  EXAMPLE 

Description 

Clears  ALL  data  in  the  Toolkit  memory  banks  and  pastes  an  example 
dataset.  The  example  dataset  used  in  the  Toolkit  is  obtained  from  the  API’s 
“Groundwater  Remediation  Strategies  Tool”  (Newell  et  a!.,  2003). 

PARAMETER 

LOAD  FILE 

Description 

Loads  data  files  saved  through  the  Toolkit.  DO  NOT  EDIT  ANY  TOOLKIT 
FILES.  Editing  files  may  cause  the  Toolkit  to  crash. 
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PARAMETER 

SAVE  FILE 

Description 

Saves  all  Transect  and  Receptor  Impact  data  entered  into  the  Toolkit.  DO 

NOT  ADD  ANY  EXTENSTIONS  TO  FILE  NAME  WHEN  SAVING. 

PARAMETER 

PRINT 

Description 

Prints  the  data  shown  on  the  screen  on  the  default  printer.  To  print  on  a 
different  printer,  select  the  printer  in  the  “Print”  options  in  Excel  and  then 
press  the  “Print”  button  in  the  Toolkit. 
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Input  Data  and  Grid:  Choose  Transect 


PARAMETER  TRANSECT 


Description  The  plane  perpendicular  to  groundwater  flow  representing  vertical  distribution 

of  plume  concentration. 

To  characterize  mass  flux,  transects  should  be  located  at  points  where 
sufficient  data  are  available  to  define  affected  groundwater  concentrations 
and  specific  discharge.  For  two-dimensional  data  (i.e. ,  from  single-screen 
monitoring  wells),  the  transect  will  represent  a  line  extending  across  the  full 
width  of  the  plume,  perpendicular  to  the  direction  of  groundwater.  For  uniform 
flow  fields,  this  transect  will  be  a  straight  line,  but,  for  converging  or  diverging 
flowlines,  the  transect  will  be  curvilinear  in  shape.  For  three-dimensional  data 
(i.e.,  from  multilevel  monitoring  wells),  the  transect  line  will  represent  a 
vertical  plane  through  the  groundwater  plume,  positioned  perpendicular  to 
groundwater  flow. 


Transects  through  2-D  Plume. 


Transect  1  Transect  2 


Transects  through  3-D  Plume  Delineation. 


How  to  Enter  Data 


Choose  the  appropriate  transect  from  the  drop  down  list.  Data  for  up  to  five 
different  plume  transects  can  be  entered  into  the  Toolkit. 


Input  Data  and  Grid:  Choose  Time  Period 


PARAMETER 

TIME  PERIOD 

Description 

Data  collection  time  period. 
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How  to  Enter  Data  Choose  the  appropriate  time  period  from  the  drop  down  list.  At  least  one 

concentration  value  must  be  entered  for  the  data  to  be  saved.  Data  for  up  to 
five  different  time  periods  can  be  entered  into  the  Toolkit 

When  entering  data  for  an  additional  time  period  for  a  particular  transect,  the 
user  has  the  option  of  automatically  keeping  previously  entered  data  (such  as 
the  well  name,  sampling  interval,  etc).  Select  the  parameter  information  to 
keep  from  the  “Data  to  Keep”  options  window. 


Input  Data  and  Grid:  Enter  Transect  Data 


PARAMETER 

DISTANCE  OF  TRANSECT  FROM  SOURCE 

Units 

ft  (or  m) 

Description 

Distance  of  the  plume  transect  of  interest  from  the  source.  The  source  zone 
is  defined  as  the  zone  that  includes  the  affected  soils  in  both  the  vadose  and 
smear  zones.  This  value  is  only  used  in  a  summary  plot  that  shows  mass 
flux  vs.  distance  from  the  plume  (maximum  of  five  different  transects) 

How  to  Enter  Data 

Enter  directly. 

PARAMETER 

DARCY  VELOCITY 

Units 

cm/sec,  ft(or  m)/day,  ft(or  m)/yr 

Description 

Groundwater  Darcy  velocity  through  the  transect  of  interest. 

To  characterize  the  specific  discharge  across  each  plume  transect, 
representative  measurements  of  the  Darcy  velocity  should  be  obtained  at  one 
or  more  locations,  using  appropriate  slug  or  pumping  test  methods. 

Typical  Values 

0.2  -  200  ft/yr  (Newell  et  at.,  1 996) 

Source  of  Data 

Calculated  by  multiplying  hydraulic  conductivity  by  hydraulic  gradient  (Vd  =  K 
*  i).  Use  of  actual  site  data  for  hydraulic  conductivity  and  hydraulic  gradient 
parameters  is  strongly  recommended. 

How  to  Enter  Data 

1)  Select  appropriate  radio  button, 

2)  Select  units,  and 

3)  Enter  directly: 

a)  As  a  uniform  value  for  the  entire  transect,  or 

b)  Individually  at  each  sampling  depth  for  each  monitoring  point. 

PARAMETER 

HYDRAULIC  CONDUCTIVITY 

Units 

cm/sec,  ft(or  m)/day,  ft(or  m)/yr 

Description 

Measure  of  the  permeability  of  the  saturated  porous  medium. 

To  characterize  the  specific  discharge  across  each  plume  transect, 
representative  measurements  are  required  for  both  the  hydraulic  flow 
gradient  and  the  hydraulic  conductivity  of  the  flow  system.  The  groundwater 
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flow  direction  and  hydraulic  gradient  for  each  segment  of  a  transect  line  can 
be  determined  from  a  potentiometric  surface  contour  map  based  on  static 
water  level  measurements  of  available  sampling  points.  Representative 
measurements  of  the  hydraulic  conductivity  of  the  groundwater-bearing  unit 
should  be  obtained  at  one  or  more  locations,  using  appropriate  slug  test  or 
pumping  test  methods  (Newell  et  al.,  2003). 

Typical  Values  Clays:  <1x10"®  cm/s 

Silts:  IxlO"6  -1x10"3  cm/s 

Siltysands:  1x10"^  -  1x10"  ^  cm/s 

Clean  sands:  1x10"  ^  -1  cm/s 

Gravels:  >  1  cm/s 

(Newell  ef  al.,  1996) 

Source  of  Data  Pump  tests  or  slug  tests  at  the  site.  It  is  strongly  recommended  that  actual 

site  data  be  used  for  all  mass  flux  studies. 

How  to  Enter  Data  1 )  Select  appropriate  radio  button, 

2)  Select  units,  and 

3)  Enter  directly: 

a)  As  a  uniform  value  for  the  entire  transect,  or 

b)  Individually  at  each  sampling  depth  for  each  monitoring  point. 


PARAMETER 

HYDRAULIC  GRADIENT 

Units 

ft/ft  (or  m/m) 

Description 

The  slope  of  the  potentiometric  surface.  In  unconfined  aquifers,  this  is 
equivalent  to  the  slope  of  the  water  table. 

Typical  Values 

0.0001-  0.1  ft/ft 

Source  of  Data 

Calculated  by  constructing  potentiometric  surface  maps  using  static  water 
level  data  from  monitoring  wells  and  estimating  the  slope  of  the 
potentiometric  surface. 

How  to  Enter  Data 

Enter  directly: 

a)  As  a  uniform  value  for  the  entire  transect,  or 

b)  Individually  at  each  sampling  depth  for  each  monitoring  point. 

PARAMETER 

SAMPLING  INTERVAL 

Units 

ft  bgs,  ft  MSL,  m  bgs,  or  m  MSL 

Description 

The  vertical  interval  from  which  concentration  data  is  obtained.  This  option  is 
used  when  most  or  all  of  the  monitoring  points  have  defined  screen  intervals. 
The  sampling  interval  is  the  screening  interval  of  the  monitoring  well  in  units 
of  distance  below  ground  surface  or  elevation  in  MSL.  (For  direct  push-type 
samples,  used  the  Mid  Point  of  Sampling  Interval  option). 

(Note:  if  you  have  mixtures  of  monitoring  wells  (with  screened  intervals)  or 
monitoring  points  (just  a  single  sampling  point  depth),  you  need  to  have  all 
points  be  represented  with  either  the  Sampling  Interval  method  or  the  Mid 
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Point  of  Sampling  Interval  method). 

For  each  selected  plume  transect,  sufficient  groundwater  sampling  points 
must  be  available  to  define  i)  the  full  width  and  thickness  of  the  plume  and  ii) 
the  distribution  of  contaminant  concentrations  within  the  plume.  Either  single- 
level  or  multilevel  groundwater  monitoring  points  may  be  used  for  this 
purpose.  Multilevel  monitoring  points  can  provide  a  more  detailed  three- 
dimensional  characterization  of  contaminant  concentrations  in  groundwater. 
However,  single-level  groundwater  monitoring  networks,  while  less  accurate 
than  multi-level  networks,  can  still  provide  sufficient  accuracy  to  support  a 
mass  flux  analysis  at  many  sites  (Newell  et  at.,  2003). 

Typical  Values  3 -20  ft 

Source  of  Data  Well  installation  logs. 

How  to  Enter  Data  1 )  Select  appropriate  radio  button,  and 

2)  Enter  directly. 


PARAMETER  MID  POINT  OF  SAMPLING  INTERVAL 

Units  ft  bgs  (or  m  bgs),  ft  MSL  (or  m  MSL) 

Description  The  mid  point  of  the  interval  from  which  concentration  data  is  obtained.  This 

option  is  used  when  you  have  a  direct  push  type  sample  with  a  very  small 
screen  that  can  be  considered  a  single  value.  Typically  this  is  the  mid  point 
of  screening  interval  of  the  monitoring  point.  For  monitoring  wells  with 
screened  intervals,  you  can  use  the  mid-  point  of  the  screened  interval,  or 
use  the  Sampling  Interval  method. 

(Note:  if  you  have  mixtures  of  monitoring  wells  (with  screened  intervals)  or 
monitoring  points  (just  a  single  sampling  point  depth),  you  need  to  have  all 
points  be  represented  with  either  the  Sampling  Interval  method  or  the  Mid 
Point  of  Sampling  Interval  method). 

For  each  selected  plume  transect,  sufficient  groundwater  sampling  points 
must  be  available  to  define  i)  the  full  width  and  thickness  of  the  plume  and  ii) 
the  distribution  of  contaminant  concentrations  within  the  plume.  Either  single- 
level  or  multilevel  groundwater  monitoring  points  may  be  used  for  this 
purpose.  Multilevel  monitoring  points  can  provide  a  more  detailed  three- 
dimensional  characterization  of  contaminant  concentrations  in  groundwater. 
However,  single-level  groundwater  monitoring  networks,  while  less  accurate 
than  multi-level  networks,  can  still  provide  sufficient  accuracy  to  support  a 
mass  flux  analysis  at  many  sites  (Newell  et  at.,  2003). 

Source  of  Data  Well  installation  logs. 

How  to  Enter  Data  1 )  Select  appropriate  radio  button,  and 

2)  Enter  directly. 

The  maximum  number  of  lines  of  data  that  can  be  entered  is  250. 


PARAMETER 

MONITORING  POINT 

Description 

For  each  selected  plume  transect,  sufficient  groundwater  sampling  points 
must  be  available  to  define  i)  the  full  width  and  thickness  of  the  plume  and  ii) 
the  distribution  of  contaminant  concentrations  within  the  plume.  Either  single- 
level  or  multilevel  groundwater  monitoring  points  may  be  used  for  this 
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purpose.  Multilevel  monitoring  points  can  provide  a  more  detailed  three- 
dimensional  characterization  of  contaminant  concentrations  in  groundwater. 
However,  single-level  groundwater  monitoring  networks,  while  less  accurate 
than  multi-level  networks,  can  still  provide  sufficient  accuracy  to  support  a 
mass  flux  analysis  at  many  sites. 

(NOTE:  an  alternative  method  is  to  use  the  concentrations  from  contour  lines 
as  measurement  points.  The  locations  where  the  transect  intersects  contour 
lines  from  plume  maps  can  be  used  to  construct  flow  areas  for  the  mass  flux 
calculations)  (Newell  et  a/.,  2003). 

Source  of  Data 

Plume  maps. 

How  to  Enter  Data 

Enter  directly. 

The  maximum  number  of  monitoring  points  that  can  be  entered  is  1 00,  while 
the  maximum  number  of  lines  of  data  that  can  be  entered  is  250. 

PARAMETER 

START  OF  TRANSECT 

Units 

ft  (or  m) 

Description 

The  Mass  Flux  Toolkit  calculates  mass  flux  based  on  a  two-dimensional  grid 
obtained  by  dividing  the  transect  into  subareas.  For  this  purpose,  an  edge  of 
the  transect  is  defined  as  the  origin  and  the  locations  of  all  monitoring  points 
based  on  that  origin.  The  start  of  the  transect  is  automatically  designated  as 
the  origin,  the  user  must  provide  a  distance  representing  the  total  length  of 
the  transect. 

q  start  of  transect 

/  _ 

MW-1  /  O 

- - MW-5 

^ - &  MW-6 

°  - - 

end  of  O 

transect  " 

Transects  through  2-D  Plume. 

PARAMETER 

END  OF  TRANSECT 

Units 

ft  (or  m) 

Description 

Total  length  of  transect. 

The  Mass  Flux  Toolkit  calculates  mass  flux  based  on  a  two-dimensional  grid 
obtained  by  dividing  the  transect  into  subareas.  For  this  purpose,  an  edge  of 
the  transect  is  defined  as  the  origin  and  the  locations  of  all  monitoring  points 
based  on  that  origin.  The  start  of  the  transect  is  automatically  designated  as 
the  origin,  the  user  must  provide  a  distance  representing  the  total  length  of 
the  transect. 
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q  start  of  transect 

/  _ 

MW-1  /  O 

MW.5 

Gl/1/p.  ^ - &  MW-6 

0  / - 

end  of  O 

transect  ~~ 

Transects  through  2-D  Plume. 

Typical  Values 

10-500  ft 

Source  of  Data 

To  determine  a  transect,  draw  a  line  perpendicular  to  the  direction  of 
groundwater  flow  in  the  plume  area.  Define  either  the  left  or  the  right  edge  of 
the  transect  as  the  origin.  Calculate  the  distances  of  all  monitoring  points 
defining  the  transect  from  that  origin. 

How  to  Enter  Data 

Enter  directly.  Note  that  the  ‘End  of  Transect’  distance  must  be  greater  than 
the  distance  of  the  farthest  well  from  the  start  of  the  transect. 

PARAMETER 

DISTANCE  OF  MONITORING  POINT  FROM  START  OF  TRANSECT 

Units 

ft  (or  m) 

Description 

Distance  of  monitoring  point  from  start  of  transect. 

The  Mass  Flux  Toolkit  calculates  mass  flux  based  on  a  two-dimensional  grid 
obtained  by  dividing  the  transect  into  subareas.  For  this  purpose,  an  edge  of 
the  transect  is  defined  as  the  origin  and  the  locations  of  all  monitoring  points 
based  on  that  origin.  The  start  of  the  transect  is  automatically  designated  as 
the  origin,  the  user  must  provide  a  distance  representing  the  total  length  of 
the  transect. 

q  start  of  transect 

MW-1  /  O 

- MW.5 

GlVp.  ^ - &  MW-6 

°  - 

end  of  O 

transect  " 

Transects  through  2-D  Plume. 

Typical  Values 

10-500  ft 

Source  of  Data 

To  determine  a  transect,  draw  a  line  perpendicular  to  the  direction  of 
groundwater  flow  in  the  plume  area.  Define  either  the  left  or  the  right  edge  of 
the  transect  as  the  origin.  Calculate  the  distances  of  all  monitoring  points 
defining  the  transect  from  that  origin. 
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How  to  Enter  Data  Enter  directly.  Note  that  the  distance  of  the  first  well  from  the  “Start  of 

Transect”  must  be  >  0  while  the  ‘End  of  Transect’  must  be  greater  than  the 
distance  of  the  farthest  well  from  the  start  of  the  transect. 


PARAMETER 

PLUME  TOP 

Units 

ft  bgs  (or  m  bgs),  ft  MSL  (or  m  MSL) 

Description 

Top  of  contaminant  plume. 

For  each  selected  plume  transect,  sufficient  groundwater  sampling  points 
must  be  available  to  define  the  thickness  of  the  plume. 

Typical  Values 

5-500  ft  bgs 

Source  of  Data 

Multilevel  groundwater  monitoring  points. 

How  to  Enter  Data 

Enter  directly. 

The  maximum  number  of  lines  of  data  that  can  be  entered  is  250. 

PARAMETER 

PLUME  BOTTOM 

Units 

ft  bgs  (or  m  bgs),  ft  MSL  (or  m  MSL) 

Description 

Bottom  of  contaminant  plume. 

For  each  selected  plume  transect,  sufficient  groundwater  sampling  points 
must  be  available  to  define  the  thickness  of  the  plume. 

Typical  Values 

5-500  ft  bgs 

Source  of  Data 

Multilevel  groundwater  monitoring  points. 

How  to  Enter  Data 

Enter  directly. 

The  maximum  number  of  lines  of  data  that  can  be  entered  is  250. 

PARAMETER  CONCENTRATION 

Units  mg/L  or  ug/L 

Description  Aqueous  phase  concentration  of  constituents  in  the  transect  of  interest. 

For  each  selected  plume  transect,  sufficient  groundwater  sampling  points 
must  be  available  to  define  i)  the  full  width  and  thickness  of  the  plume  and  ii) 
the  distribution  of  contaminant  concentrations  within  the  plume.  Either  single- 
level  or  multilevel  groundwater  monitoring  points  may  be  used  for  this 
purpose.  Multilevel  monitoring  points  can  provide  a  more  detailed  three- 
dimensional  characterization  of  contaminant  concentrations  in  groundwater. 
However,  single-level  groundwater  monitoring  networks,  while  less  accurate 
than  multi-level  networks,  can  still  provide  sufficient  accuracy  to  support  a 
mass  flux  analysis  at  many  sites. 

(NOTE:  an  alternative  method  is  to  use  the  concentrations  from  contour  lines 
as  measurement  points.  The  locations  where  the  transect  intersects  contour 
lines  from  plume  maps  can  be  used  to  construct  flow  areas  for  the  mass  flux 
calculations)  (Newell  et  al.,  2003). 
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NOTE:  The  Toolkit  assumes  the  concentrations  at  the  start  and  end  of 
the  transect  are  zero. 

Typical  Values 

0.0001  -  500  mg/L 

Source  of  Data 

Monitoring  well  data. 

How  to  Enter  Data 

1 )  Enter  name  of  constituent. 

2)  Enter  concentration  values. 

The  maximum  number  of  lines  of  data  that  can  be  entered  is  250. 

PARAMETER  CHOOSE  GRID 


Description  Subareas  the  transect  plane  is  divided  into  for  mass  flux  calculations.  The 

width  of  the  grid  cells  is  obtained  from  the  distance  between  monitoring 
points.  The  thickness  of  the  cells  is  based  on  the  plume  width  and  length  of 
sampling  intervals. 

The  columns  in  the  default  grid  extend  half  way  on  either  side  between 
adjacent  monitoring  points.  The  width  of  the  columns  and  thickness  of  the 
ceils  can  be  decreased  by  refining  the  grid. 


Distance  of  first  monitoring 

point  from  edge  of  transect  Distance  from  Edge  of  Transect  (ft) 

MW-1  MW-2  MW-3  MW -4  MW-5 

0.0  10.0 *  27.5  45.0  62.5  60.0 


W 

O) 

-Q 


Cell  width 


5.0 

6.5 
8.0 

9.5 
11.0 

12.5 
14.0 

15.5 
17.0 

18.5 
20.0 


2.300 

2.300 


Data  from  monitoring 
well  screening  interval 
of  5  to  8  ft  bgs 


j  Cell  thickness  7.200 

7.200 


7.200 


35.600 

35.600 

35.600 

35.600 

9.500 

9.500 

9.500 


15.300 

15.300 

15.300 


Bottom  of  plume 


4.500 

4.500 

4.500 

4.500 


Typical  Values 


Source  of  Data 


How  to  Enter  Data 


1  x  1  to  3  x  3  (for  10  or  less  monitoring  points  per  transect) 
.Monitoring  well  data  and  plume  maps. 


1.  The  default  refinement  is  1  by  1.  To  refine  the  default  grid,  enter  the 
factors  to  divide  the  cell  thickness  and  width  with. 


2.  Enter  directly. 

NOTE:  Although  the  maximum  grid  refinement  is  10  for  rows  and  100/total 
number  of  sampling  points  for  columns  (resulting  in  a  maximum  grid  of  1 00 
rows  and  100  columns),  the  greater  the  grid  refinement  the  greater  the 
computational  time.  For  10  monitoring  points  or  less,  a  maximum  grid 
refinement  of  3  x  3  is  recommended  for  most  computers. 
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PARAMETER 

SELECT  CONSTITUENT  OF  INTEREST 

Description 

Select  the  constituent  to  use  in  the  mass  flux  calculations. 

How  to  Enter  Data 

Select  radio  button. 

PARAMETER 

CONTINUE  DATA  INPUT 

Description 

Creates  the  concentration  grid  to  be  used  in  the  mass  flux  calculations. 

PARAMETER 

EXPORT  MW  DATA 

Description 

The  Mass  Flux  Toolkit  allows  the  monitoring  well  data  entered  on  the  Input 

Data  and  Grid  screen  to  be  exported  into  a  text  file  for  use  in  other  programs. 
The  Toolkit  exports  data  one  transect  and  one  time  period  at  a  time  into  tab 
delimited  text  file. 

Import/Export  File  Format 


Transect 

Number 

(enter  vabe) 

Time  Poriod 
Number 

(enter  value) 

Syelem  Unit* 

1  =  En^ish 

2  =  SI 

Elevation  Unite 

1  =  bgs 

2  =  MSI 

Ground 

Surface 

Elevation 

(errter  value  if 

MSL) 

Concentration 

Units 

1  = 

2  =  u^L 

Dercy  Or 
Conductivity? 

1  =  Darcy 

2  =  Conductivity 

Darcy/ 

Conductivity 

Unite 

1  =  cm/sec 

2  =  ft/d 

3  =  ft/yr 

Uniform  Darcy1 
Conductivity? 

1  =  yes 

2  «  no 

Uniform  Darcy/ 

Conductivity 

Value 

(erter  value  if 
uniform) 

Uniform 

Gradient? 

1  =  yes 

2  *  no 

(leave  Wank  if 
Darcy) 

Uniform 

Gradient 

Value 

(erter  value  if 
uniform  or 
leave  Wank  if 
Darcy) 

Distance  From 
Source 

(enter  vn'ue) 

Sampling 

Interval  or 
Midpoint? 

1  =  Interval 

2  =  Midpoint 

End  of 

Traneeet 

(enter  value) 

Name 

Dietance  from 
Start  of 

Traneeet 

Sampling 
Interval  Top 

Sampling 

Interval 

Bottom 

Midpoint 

of 

Sampling 

Interval 

Plume 

Top 

Plume 

Bottom 

Leave 

Column 

Blank 

Leave 

Column 

Blank 

Darcy 

Velocity 

Hydraulic 

Conductivity 

Hydraulic 

Gradient 

COC_A 

COC_B 

Sample  File 


|  4 

1  1  1 

1 

1 

1 

2 

3 

1 

rO032 

1  1 

0.002 

864  3 

1 

60 

Name 

Distance 

InlervalTcp 

Interval  Bottom 

FlumeTcp 

FlumeS  ol 

Conductivity  Gradient 

COC.A 

MW -2 

10 

6 

id  I 

6 

20 

0032  0002 

MW -2 

Id  1 

10 

15 

5 

20 

!  0  082  0  002 

005 

MW -2 

10 

IS 

20 

6 

20 

0  CG2  I  0.002 

005 

7RZ-8 

as _ ~~ 

6 

10  ! 

6 

20 

)  03? 

2  67 

TRZ-8 

35 

10 

MS 

5 

20 

~ 0082  fOCO? 

142 

TftZ-6 

36 

16  1 

20 

6 

20 

0  032  0002 

005 

TRZ-19 

67  6 

6 

id 

6 

20 

0032  i  0002 

1  13 

TRZ-19 

57  5 

10 

15 

5 

20 

0  032  0.002 

0  68 

TRZ-19 

67  6 

16 

20 

5 

20 

0  032  0002 

006 

PARAMETER 

IMPORT  MW  DATA 

Description 

The  Mass  Flux  Toolkit  allows  monitoring  well  data  to  be  imported  from  a  tab- 
delimited  text  file.  The  Toolkit  imports  data  one  transect  and  one  time  period 
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at  a  time. 


Imported  File  Format 

See  Exported  File  Format  above. 

How  to  Enter  Data 

Enter  values  in  cells  where  data  is  available.  Leave  cells  with  no  data  blank. 

PARAMETER 

BACK  TO  TRANSECT  MAIN 

Description 

Returns  to  the  Mass  Flux  Transect  Calculator  main  screen. 

PARAMETER 

CLEAR  SCREEN 

Description 

Clears  the  data  in  the  Toolkit  memory  banks  for  the  selected  transect  and 
time  period.  Use  this  button  to  clear  the  screen. 

PARAMETER 

PASTE  EXAMPLE 

Description 

Clears  ALL  data  in  the  Toolkit  memory  banks  and  pastes  an  example 
dataset.  The  example  dataset  used  in  the  Toolkit  is  obtained  from  the 

American  Petroleum  Institute’s  (API)  “Groundwater  Remediation  Strategies 
Tool”  (Newell  et  al„  2003). 

PARAMETER 

RESTORE  TABLE  FORMATTING 

Description 

If  you  pasted  data  into  the  table  portion  of  the  screen  from  another  excel 
sheet,  then  this  button  will  restore  the  Toolkit  formatting  to  the  table. 

PARAMETER 

SEE  SAVED  GRIDS 

Description 

Proceeds  to  the  Final  Concentration  Grid  screen  if  any  mass  flux  has  been 
calculated  previously. 

PARAMETER 

PRINT 

Description 

Prints  the  data  shown  on  the  screen  on  the  default  printer.  To  print  on  a 
different  printer,  select  the  printer  in  the  “Print”  options  in  Excel  and  then 
press  the  “Print”  button. 

PARAMETER 

SELECT  DATA  TO  KEEP  FROM  PREVIOUS  TIME  PERIOD 

Description 

Gives  the  user  the  option  of  keeping  data  already  entered  for  a  Transect 

Time  Period  for  the  next  Time  Period.  This  user  can  select  to  copy  over  all 
the  listed  parameters,  none  of  the  parameters,  or  customize  the  selection  of 
parameters  to  copy  from  the  preceding  Time  Period. 
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Method  for  Interpolating  Grid 

PARAMETER  CHOOSE  A  METHOD  FOR  INTERPOLATING  GRID 

Description  The  user  specifies  which  of  three  averaging  methods  will  be  used  to 

interpolate  the  concentration,  conductivity/Darcy  velocity,  and  gradient  grids: 
user-defined,  nearest  neighbor,  or  interpolation  (linear  interpolation  or  log 
transformation).  Interpolated  concentration  data  can  be  edited  by  directly 
selecting  the  desired  cell  in  the  grid. 

User  Defined 

When  the  user-defined  method  is  chosen,  the  user  must  interpolate  the 
concentrations  manually. 

Nearest  Neighbor 

When  the  nearest  neighbor  option  is  chosen,  the  model  assigns  values  to 
each  cell  based  on  the  closest  input  data  point. 

Linear  Interpolation 

When  the  linear  interpolation  option  is  chosen,  the  grid  is  first  interpolated 
vertically  and  then  horizontally.  The  user  has  to  further  select  whether  to  use 
a)  linear  interpolation  or  b)  log  transformation  to  interpolate  the  horizontal 
cells  in  the  grid. 

Log  Transformation  Interpolation 

Under  this  option,  the  natural  logarithms  of  the  values  are  obtained, 
interpolations  performed,  and  the  interpolated  data  re  transformed. 

When  the  log  transformation  option  is  chosen,  the  grid  is  first  interpolated 
vertically  then  horizontally.  The  user  has  to  further  select  whether  to  use  a) 
linear  interpolation  or  b)  log  transformation  to  interpolate  the  horizontal  cells 
in  the  grid. 

NOTE:  The  Toolkit  also  allows  the  grid  to  be  interpolated  outside  of  the 
Toolkit.  For  this  purpose,  a)  the  “User  Defined”  option  should  be 
chosen,  b)  the  grid  “exported”  using  the  Export  Grid  function  in  the 
Toolkit,  c)  interpolated  by  user  outside  the  Toolkit,  and  d)  imported  back 
into  the  Toolkit. 

How  to  Enter  Data  1.  Select  radio  button  (Step  9). 

2.  If  the  “interpolation”  option  is  selected  then 

a.  Select  the  interpolation  scheme  for  vertical  interpolation  of  the  grid 

(Step  9A). 

b.  If  desired,  edit  interpolated  values. 

c.  Select  the  interpolation  scheme  for  horizontal  interpolation  of  the  grid 

(Step  10). 

3.  Interpolated  data  can  be  edited  by  directly  selecting  the  desired  cell  in  the 
grid. 
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Grid  Completion:  Interpolate  Concentration,  Conductivity,  and  Gradient 

Before  total  mass  flux  can  be  calculated,  the  following  grids  must  be  interpolated: 

1 .  Concentration 

2.  Darcy  Velocity/Hydraulic  Conductivity 

3.  Hydraulic  Gradient  (if  necessary) 

To  get  back  to  the  previously  calculated  mass  flux  grid  do  the  following: 

a.  Press  Back  to  Data  Input  button 

b.  Press  See  Conc/Flux  Grids'  button 

c.  Press  View  Mass  Flux  Result’  button 


PARAMETER 

INTERPOLATE  HORIZONTALLY 

Description 

When  the  linear  or  log  transformation  interpolation  option  is  chosen  in  Step  9, 
the  user  has  to  select  whether  to  use  a)  linear  interpolation  or  b)  log 
transformation  to  interpolate  the  horizontal  rows  of  the  grid. 

How  to  Enter  Data 

1 .  Select  the  interpolation  scheme  for  horizontal  interpolation  of  the  grid  (Step 

10). 

2.  Interpolated  data  can  be  edited  by  directly  selecting  the  desired  cell  in  the 
grid. 

PARAMETER 

EXPORT  GRID 

Description 

The  Mass  Flux  Toolkit  allows  the  exporting  of  the 
concentration/conductivity/gradient  grids  into  a  text  file  for  external 
interpolation. 

Note  that  the  Toolkit  only  exports  the  grid  currently  shown  on  the  screen. 

PARAMETER 

IMPORT  GRID 

Description 

The  Mass  Flux  Toolkit  allows  the  exporting/importing  of  the 
concentration/conductivity/gradient  grids  into  a  text  file  for  external 
interpolation. 

NOTE:  Any  externally  interpolated  grid  being  imported  into  the  Toolkit 
MUST  have  the  same  number  of  rows  and  columns  as  the  Toolkit  grid. 
Additionally,  uncertainty  analysis  is  not  available  for  imported  grids. 

PARAMETER 

CALCULATE  MASS  FLUX 

Description 

The  total  contaminant  mass  flux  across  the  transect  is  calculated  as  follows: 

w  =  H«QiACF 

where: 
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w  =  total  mass  flux  from  source  zone  (g/day)  (also  called  mass  discharge) 
C,  =  concentration  of  constituent  at  flow  area  in  transect  (mg/L) 

qt  =  specific  discharge  (also  called  Darcy  velocity)  through  flow  area 
associated  with  an  individual  constituent  measurement  i  (cm/sec),  qi  can 
be  calculated  using 

q,  —  K  -  i 

where: 

K  =  hydraulic  conductivity  (cm/sec) 
i  =  hydraulic  gradient  (cm/cm) 

A i  =  Flow  area  associated  with  an  individual  constituent  measurement  (ft2) 
CF  =  conversion  factor  =  80.3  (ft/cm )/(sec/day)(l_/ft3)(g/mg) 


PARAMETER 

BACK  TO  DATA  INPUT 

Description 

Returns  to  the  Input  Data  and  Grid  screen. 

PARAMETER 

BACK  TO  CONC  GRID 

Description 

Returns  to  the  Grid  Completion:  Concentration  screen. 

PARAMETER 

BACK  TO  K/DARCY  GRID 

Description 

Returns  to  the  Grid  Completion:  Hydraulic  Conductivity/Darcy  Velocity 
screen. 

PARAMETER 

PRINT 

Description 

Prints  the  data  shown  on  the  screen  on  the  default  printer.  To  print  on  a 
different  printer,  select  the  printer  in  the  “Print'’  options  in  Excel  and  then 
press  the  “Print”  button. 
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Mass  Flux  Results 


PARAMETER 

TOTAL  MASS  FLUX 

Description 

The  total  contaminant  mass  flux  across  the  transect  in  g/day  and  kg/yr. 

PARAMETER 

MASS  FLUX  SUMMARY 

Description 

Takes  the  user  to  Mass  Flux  Summary  screen. 

PARAMETER 

RUN/VIEW  UNCERTAINTY  ANALYSIS 

Description 

Uncertainty  in  mass  flux  estimates  is  a  key  issue  in  using  mass  flux  as  a 
metric.  The  Toolkit  provides  two  options  for  analyzing  uncertainty  in  the  total 
mass  flux  estimates  derived  from  the  transect  method. 

One  option  (Simple  Uncertainty  Analysis)  provides  a  tool  for  estimating  the 
contribution  of  each  individual  observation  to  the  total  mass  flux.  If  a  single 
monitoring  point  represents  a  high  percentage  of  the  total  mass  flux,  then  the 
uncertainty  in  the  calculation  is  high  and  additional  monitoring  points  should 
be  added  to  reduce  the  uncertainty. 

The  second  option  (Advanced  Uncertainty  Analysis)  utilizes  a  Monte  Carlo 
type  approach  to  analyze  uncertainty  in  the  actual  concentration,  hydraulic 
conductivity,  and  gradient  measurements.  With  this  tool,  groundwater 
practitioners  can  estimate  the  accuracy  of  the  hydrologic  measurements  that 
are  being  used  for  the  mass  flux  calculation. 

PARAMETER 

SELECT  TRANSECT  TO  VIEW 

Description 

Select  a  different  transect  to  view.  This  option  allows  the  user  to  see 
previously  calculated  flux  grids. 

PARAMETER 

SELECT  TIME  PERIOD  TO  VIEW 

Description 

Select  a  different  time  period  to  view.  This  option  allows  the  user  to  see 
previously  calculated  flux  grids. 

PARAMETER 

VIEW  CONCENTRATION  GRID 

Description 

Takes  the  user  to  the  Final  Concentration  Grid  screen.  The  concentration 
values  in  this  grid  cannot  be  edited.  To  edit  concentration  values  return  to 
the  “Grid  Completion”  screen. 

PARAMETER 

BACK  TO  DATA  GRID 

Description 

Returns  to  the  Input  Data  and  Grid  screen. 
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Final  Concentration  Grid 


PARAMETER 

SELECT  TRANSECT  TO  VIEW 

Description 

Select  a  different  transect  to  view.  This  option  allows  the  user  to  see 
previous  grids. 

PARAMETER 

SELECT  TIME  PERIOD  TO  VIEW 

Description 

Select  a  different  time  period  to  view.  This  option  allows  the  user  to  see 
previous  grids. 

PARAMETER 

MASS  FLUX  SUMMARY 

Description 

Takes  the  user  to  Mass  Flux  Summary  screen. 

PARAMETER 

VIEW  MASS  FLUX  RESULTS 

Description 

Takes  the  user  to  the  Mass  Flux  Results  screen. 

PARAMETER 

VIEW  FINAL  CONC  GRID 

Description 

Shows  the  user  the  Final  Concentration  grid.  This  screen  is  for  viewing 
purposes  only  and  data  cannot  be  edited.  To  edit  grid  data  use  the  ‘Back  to 
Data  Grid’  button;  to  edit  monitoring  well  input  data  use  the  ‘Back  to  Data 

Input’  button. 

PARAMETER 

VIEW  FINAL  K/DARCY  GRID 

Description 

Shows  the  user  the  Final  Hydraulic  Conductivity  or  Darcy  Velocity  grid.  This 
screen  is  for  viewing  purposes  only  and  data  cannot  be  edited.  To  edit  grid 
data  use  the  ‘Back  to  Data  Grid’  button;  to  edit  monitoring  well  input  data  use 
the  ‘Back  to  Data  Input’  button. 

PARAMETER 

VIEW  FINAL  1  GRID 

Description 

Shows  the  user  the  Final  Hydraulic  Conductivity  grid.  This  screen  is  for 
viewing  purposes  only  and  data  cannot  be  edited.  To  edit  grid  data  use  the 
‘Back  to  Data  Grid’  button;  to  edit  monitoring  well  input  data  use  the  ‘Back  to 
Data  Input’  button. 

PARAMETER 

PRINT 

Description 

Prints  the  data  shown  on  the  screen  on  the  default  printer.  To  print  on  a 
different  printer,  select  the  printer  in  the  “Print”  options  in  Excel  and  then 
press  the  “Print”  button. 
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Mass  Flux  Summary 


PARAMETER 

SUMMARY  GRAPH 

Description 

The  summary  graph  plots  mass  flux  vs.  distance  to  show  the  effect  of 
remediation/impact  of  natural  attenuation  processes. 

A  maximum  of  two  types  of  error  bars  are  shown  on  the  graph.  The  first  one 
(interpolation  error  bar)  represents  the  minimum  and  maximum  mass  fluxes 
calculated  using  all  the  three  available  interpolation  schemes  in  the  Toolkit. 

The  other  set  of  error  bars  (input  interpolation)  represent  the  15th  and  85th 
percentiles  from  the  Advanced  Uncertainty  analysis  (for  the  particular 
interpolation  scheme(s)  chosen  by  the  user  to  interpolate  concentration, 
conductivity,  and  gradient).  This  set  of  error  bars  will  only  be  shown  when  an 
Advanced  Uncertainty  analysis  is  performed  by  the  user. 

PARAMETER 

RETURN  TO  TOOLKIT  MAIN  SCREEN 

Description 

Returns  to  the  Mass  Flux  Toolkit  main  screen. 

PARAMETER 

RETURN  TO  TRANSECT  MAIN  SCREEEN 

Description 

Returns  to  the  main  screen  of  the  Mass  Flux  Transect  Calculator. 

PARAMETER 

VIEW  MASS  FLUX  RESULTS 

Description 

Returns  to  the  Mass  Flux  Results  screen. 

PARAMETER 

VIEW  CONC  GRID 

Description 

Returns  to  the  Final  Concentration  Grid  screen.  The  concentration  values  in 
this  grid  cannot  be  edited.  To  edit  concentration  values  return  to  the  “Grid 
Completion”  screen. 

PARAMETER 

PRINT 

Description 

Prints  the  data  shown  on  the  screen  on  the  default  printer.  To  print  on  a 
different  printer,  select  the  printer  in  the  “Print”  options  in  Excel  and  then 
press  the  “Print”  button. 
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Uncertainty  Analysis:  Perform  Uncertainty  Analysis 

Uncertainty  in  mass  flux  estimates  is  a  key  issue  in  using  mass  flux  as  a  metric.  The  Toolkit 
provides  two  options  for  analyzing  Type  1  uncertainty  in  the  total  mass  flux  estimates  derived 
from  the  transect  method. 

One  option  (Simple  Uncertainty  Analysis)  provides  a  tool  for  estimating  the  contribution  of  each 
individual  observation  to  the  total  mass  flux.  If  a  single  monitoring  point  represents  a  high 
percentage  of  the  total  mass  flux,  then  the  uncertainty  in  the  calculation  is  high  and  additional 
monitoring  points  should  be  added  to  reduce  the  uncertainty. 

The  second  option  (Advanced  Uncertainty  Analysis)  utilizes  the  Monte  Carlo  approach  to  analyze 
uncertainty  in  the  actual  concentration,  hydraulic  conductivity,  and  gradient  measurements.  With 
this  tool,  groundwater  practitioners  can  estimate  the  accuracy  of  the  hydrologic  measurements 
that  are  being  used  for  the  mass  flux  calculation. 


PARAMETER 

SIMPLE  UNCERTAINTY  ANALYSIS  (EVALUATE  UNCERTAINTY  IN 
INTERPOLATION  SCHEME) 

Description 

This  module  applies  crossvalidation  to  evaluate  the  uncertainty  in  the 
interpolation  scheme.  Crossvalidation  is  a  method  of  removing  one 
observation  at  a  time,  estimating  a  value  for  the  removed  observation  using 
the  remaining  observations,  and  calculating  the  output  (Olea,  1999).  The 
removed  value  is  then  replaced  and  the  next  observation  removed.  This 
process  is  continued  until  all  the  observations  have  been  removed  one  at  a 
time. 

The  Toolkit  removes  each  observed  value  of  concentration,  hydraulic 
conductivity,  and  gradient  one  at  a  time  during  crossvalidation,  creates  and 
interpolates  the  grid  using  the  same  options  chosen  in  the  ‘Grid  Completion’ 
screen,  and  calculates  mass  flux. 

How  to  Enter  Data 

Press  ‘Evaluate  Uncertainty  in  Interpolation  Scheme’. 

NOTE:  The  default  interpolation  scheme  that  the  Toolkit  uses  is  the  nearest 
neighbor  interpolation.. 

Uncertainty  analysis  assumes  that  the  Toolkit  interpolated  concentrations 
have  not  been  manually  edited.  Uncertainty  analysis  is  not  available  for 
manually  created  grids. 

PARAMETER 

ADVANCED  UNCERTAINTY  ANALYSIS  (EVALUATE  HOW 

UNCERTAINTY  IN  INPUT  DATA  AFFECTS  TOTAL  MASS  FLUX) 

Description 

This  module  uses  Latin  hypercube  sampling  (LHS),  a  modified  Monte  Carlo 
approach,  to  analyze  uncertainty  in  the  actual  concentration,  hydraulic 
conductivity,  and  gradient  measurements. 

In  the  Monte  Carlo  type  approach,  a  random  number  is  generated  for  every 
value  of  concentration,  hydraulic  conductivity,  and  hydraulic  gradient  entered 
by  the  user.  This  set  of  random  inputs  is  then  used  to  create  and  interpolate 
the  respective  grids  and  mass  flux  calculated  from  those  grids.  Repeating 
this  procedure  a  large  number  times  yields  a  probability  distribution  of  the 
mass  flux  from  which  statistical  characteristics  such  as  mean,  percentile,  and 
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variance  can  be  obtained. 

The  Mass  Flux  Toolkit  performs  10  LHS  iterations  of  20  equal  probability 
intervals  (See  Appendix  1  for  details  on  Latin  hypercube  sampling). 

How  to  Enter  Data 

1)  Specify  global  parameters  for  the  concentration,  Darcy/Conductivity,  and 
gradient.  The  Toolkit  assumes  that  the  values  entered  in  the  Input 
screen  are  the  mean  values: 

a.  Choose  a  probability  distribution  for  each  parameter  (see  Appendix 

A  of  the  User’s  Manual  for  details  on  probability  distributions). 

b.  For  the  normal  distribution,  specify  the  standard  deviation  as  a 
percent  of  the  mean.  For  lognormal  distributions,  specify  the 
error  factor  (EF)  (the  ratio  of  the  95th  percentile  to  the  median 
of  the  lognormal  data  or  the  ratio  of  the  median  to  the  5th 
percentile)  (NOTE:  the  error  factor  MUST  be  greater  than  one). 
For  uniform  distribution,  specify  the  lower  and  upper  limits  as 
percentages  of  the  mean. 

2)  If  desired,  Monte  Carlo  parameters  for  each  individual  value  of 
concentration,  Darcy/Conductivity,  and  gradient  can  be  specified  using 
the  ‘Edit  Individual  Values  Manually’  option. 

3)  Perform  Input  Uncertainty  Analysis. 

NOTE:  Uncertainty  analysis  is  performed  on  the  transect  and  time  period 
selected  in  the  Mass  Flux  Results  screen. 

Uncertainty  analysis  assumes  that  the  Toolkit  interpolated  concentrations 
have  not  been  manually  edited.  Uncertainty  analysis  is  not  available  for 
manually  created  grids. 

Uncertainty  Analysis:  Results 


PARAMETER 

INTERPOLATION  ERROR  RESULTS  (SIMPLE  UNCERTAINTY 

ANALYSIS) 

Description 

Comparison  of  the  total  mass  flux  calculated  using  all  the  observations  and 
the  contribution  of  each  removed  observation  to  the  total  mass  flux. 

If  a  single  monitoring  point  represents  a  high  percentage  of  the  total  mass 
flux,  then  the  uncertainty  in  the  calculation  is  high  and  additional  monitoring 
points  should  be  added  to  reduce  the  uncertainty. 

PARAMETER 

INPUT  DATA  UNCERTAINTY  ANALYSIS  RESULTS  (ADVANCED 
UNCERTAINTY  ANALYSIS) 

Description 

Summary  statistics  (minimum,  25th  percentile,  mean,  85th  percentile  , 
maximum,  and  variance)  on  the  estimate  of  uncertainty  in  the  mass  flux 
calculation  based  on  the  user’s  choice  of  interpolation  method  and 
uncertainty  in  the  input  variables  as  defined  by  their  probability  distributions, 
means,  variances,  and  ranges. 
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PARAMETER 

BACK  TO  MASS  FLUX  RESULTS 

Description 

Returns  to  the  Mass  Flux  Results  screen. 

PARAMETER 

BACK  TO  DATA  INPUT 

Description 

Returns  to  the  Data  Input  screen. 

PARAMETER 

PRINT 

Description 

Prints  the  data  shown  on  the  screen  on  the  default  printer.  To  print  on  a 
different  printer,  select  the  printer  in  the  “Print”  options  in  Excel  and  then 
press  the  “Print”  button. 
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Wells 

Streams 

API  Strategies  Tool 


Three  important  considerations  regarding  data  input  are: 

1 )  To  see  the  example  data  set  in  the  input  screen  of  the  software,  click  on  the 
Paste  Example  button  on  the  lower  right  portion  of  the  input  screen. 

2)  Because  the  Mass  Flux  Toolkit  is  based  on  the  Excel  spreadsheet,  you  have  to  click 
outside  of  the  cell  where  you  just  entered  data  or  hit  Return  before  any  of  the 
buttons  will  work.  Additionally,  REMOVING  rows  or  columns  from  input  screens  may 
cause  the  program  to  crash. 

3)  Parameters  used  in  the  model  are  to  be  entered  directly  into  the  white  cells. 

NOTE:  These  worksheets  only  work  if  mass  flux  has  been  calculated  for  at  least  one  transect  in 
the  Mass  Flux  Transect  Calculator. 
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Receptor  Impact  -  Wells:  Input  Data 


PARAMETER 

SELECT  CONSTITUENT  OF  INTEREST 

Description 

The  constituent  to  use  in  the  dilution  calculations. 

How  to  Enter  Data 

Select  radio  button. 

PARAMETER 

TRANSECT  OF  INTEREST 

Description 

The  transect  to  use  in  the  dilution  calculations 

Source  of  Data 

Plume  maps. 

How  to  Enter  Data 

Select  one  of  the  transects  for  which  mass  flux  has  been  calculated  in  the 
Toolkit. 

PARAMETER 

TIME  PERIOD  OF  INTEREST 

Description 

The  time  period  of  data  entered  in  the  Mass  Flux  Transect  Calculator. 

How  to  Enter  Data 

Select  one  of  the  time  periods  for  which  mass  flux  has  been  calculated  in  the 
Toolkit. 

PARAMETER 

DISCHARGE  RATE  OF  EACH  SUPPLY  WELL 

Units 

gpm,  gpd,  ft3(m3)/min,  ft3  (m3)/d 

Description 

Discharge  rate  of  each  supply  well. 

Typical  Values 

0.1  to  1000  gpm 

How  to  Enter  Data 

Enter  directly. 

PARAMETER 

NUMBER  OF  SUPPLY  WELLS 

Description 

Number  of  supply  wells  potentially  affected  by  the  contaminant  plume. 

Typical  Values 

1-5 

How  to  Enter  Data 

Enter  directly. 

PARAMETER 

AQUIFER  THICKNESS 

Units 

ft  (or  m) 

Description 

Estimated  thickness  of  the  aquifer  containing  the  supply  well  and  the 
contaminant  plume. 

Typical  Values 

10  to  100  ft 
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WELLS  RESULTS 


Receptor  Impact  -  Wells:  Results 


PARAMETER 

CONSTITUENT  CONCENTRATION 

Description 

The  maximum  concentration  of  a  contaminant  in  water  extracted  from  a 
supply  well  using  the  relationship  described  in  Einarson  and  Mackay  (2001): 

Osw  —  Mc/Qsw 

where 

Md  is  the  total  mass  flux  or  mass  discharge  of  the  contaminant  from  the  well, 
which  can  be  predicted  before  the  plume  reaches  the  well  using  the  entire 
mass  flux  through  upgradient  transects  and  Qsw  is  the  pumping  rate  from  the 
supply  well. 

PARAMETER 

CAPTURE  ZONE  OF  SUPPLY  WELL 

Description 

The  capture  zone  of  the  supply  well  using  the  relationship  described  by 
Javandel  and  Tsang  (1986). 

PARAMETER 

DISTANCE  OF  SUPPLY  WELL  FROM  STAGNATION  POINT 

Description 

The  distance  of  the  supply  well  from  the  stagnation  point  using  the 
relationship  described  by  Javandel  and  Tsang  (1986). 
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STREAMS  DATA  INPUT 


Receptor  Impact  -  Streams:  Input  Data 


PARAMETER 

CONTAMINANT  OF  INTEREST 

Description 

The  constituent  to  use  in  the  dilution  calculations. 

How  to  Enter  Data 

Select  from  radio  buttons.  . 

PARAMETER 

TRANSECT  OF  INTEREST 

Description 

The  transect  to  use  in  the  dilution  calculations 

Source  of  Data 

Plume  maps. 

How  to  Enter  Data 

Select  one  of  the  transects  for  which  mass  flux  has  been  calculated  in  the 
Toolkit. 

PARAMETER 

TIME  PERIOD  OF  INTEREST 

Description 

The  time  period  of  data  entered  in  the  Mass  Flux  Transect  Calculator. 

How  to  Enter  Data 

Select  one  of  the  time  periods  for  which  mass  flux  has  been  calculated  in  the 
Toolkit. 

PARAMETER 

STREAM  FLOW  RATE 

Units 

gpm,  gpd,  ft3(m3)/sec,  ft3(m3)/min,  ft3  (m3)/d 

Description 

The  flowrate  of  the  stream  approached  by  the  plume. 

Typical  Values 

0.1  -106  cfm 

Source  of  Data 

U.S.  geological  stream  gage  information 

How  to  Enter  Data 

Enter  directly. 

MASS  FLUX  TOOLKIT 
T  USER’S  MANUAL  ▼ 


39 


RECEPTOR  I  MPACT 


STREAMS  RESU  LTS 


Receptor  Impact  -  Streams:  Results 

PARAMETER  CONSTITUENT  CONCENTRATION 

Description  The  maximum  concentration  of  a  contaminant  in  water  extracted  from  a 

stream  downgradient  of  plume  using  the  relationship  described  in  Einarson 
and  Mackay  (2001): 

Osw  —  Mc/Qsw 

where 

Md  =  the  total  mass  flux  or  mass  discharge  of  the  contaminant  from  the  well, 
which  can  be  predicted  before  the  plume  reaches  the  stream  using  the  entire 
mass  flux  through  upgradient  transects,  and 

Qsw  =  the  flow  rate  of  the  stream. 
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API  STRATEGIES  TOOL 

The  Mass  Flux  Toolkit  provides  an  electronic  version  of  the  American  Petroleum  Institute’s  (API) 
Groundwater  Remediation  Strategies  Guide  (Newell  et  at.,  2003)  to  help  users  apply  mass  flux  to 
groundwater  remediation  decision-making.  Please  refer  to  the  above  document  (“API 
Groundwater  Remediation  Strategies  Tool"  tab  in  “Learn  About  Mass  Flux")  for  help. 
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CURRENT  ESTCP/SERDP  MASS  FLUX  PROJECTS 


Several  ESTCP  and  SERDP  projects  dealing  with  mass  flux  are  currently  ongoing.  The  following 
is  a  list  of  such  projects  by  lead  principal  investigator: 

•  Linda  Abriola:  Development  of  Assessment  Tools  for  Evaluation  of  the  Benefits  of  DNAPL 
Source  Zone  Treatment. 

This  project  deals  with  the  development  of  methodologies  for  reliable  prediction  or 
monitoring  of  plume  development  or  reduction  of  organic  mass  flux  under  heterogeneous 
field  conditions  after  the  source  zone  treatment. 

Details  on  this  project  can  be  obtained  at: 
http://www.serdp.org/research/CU/CU-1293.pdf 

•  Kirk  Hatfield:  Demonstration  and  Validation  of  a  Water  and  Solute  Flux  Measuring  Device 

This  project  involves  the  development  of  the  flux-meter  for  the  quantification  of 
contaminant  flux  leaving  a  site.  The  objective  being  that  risk  posed  by  a  site  can  be 
better  evaluated  and  remedial  efforts  made  by  measuring  cumulative  flux  rather  than 
contaminant  concentrations. 

Details  on  this  project  can  be  obtained  at: 
http://www.estcp.org/proiects/cleanup/2001 14o.cfm 

•  Tissa  lllangasekare:  Mass  Transfer  from  Entrapped  DNAPL  Sources  Undergoing 
Remediation:  Characterization  Methods  and  Prediction  Tools. 

This  project  seeks  to  develop  prediction  tools  and  site  characterization  methods  to  assist 
site  managers  1)  make  decisions  on  managing  sites  and  2)  implement  cost  and  benefit 
effective  remediation  technologies.  The  objective  being  to  understand,  quantify,  and 
model  the  transfer  of  mass  across  heterogeneous  DNAPL  source  zones  undergoing 
remediation. 

Details  on  this  project  can  be  obtained  at: 
http://www.serdp.org/research/cu/cu-1294.pdf 

•  Michael  Kavanaugh:  Diagnostic  Tools  for  Performance  Evaluation  of  Innovative  In-Situ 
Remediation  Technologies  at  Chlorinated  Solvent-Contaminated  Sites 

The  objective  of  this  project  is  to  evaluate  in-situ  remediation  technologies  at  three 
hydrogeologically  different  sites  by  applying  mass  flux  as  the  metric  of  overall  system 
performance.  Furthermore,  the  project  compares  various  methods  of  calculating  mass 
flux. 

Details  on  this  project  can  be  obtained  at: 
http://www.estcp.org/proiects/cleanup/CU-0318.cfm 

•  Lynn  Wood:  Impact  of  DNAPL  Source  Zone  Treatment:  Experimental  and  Modeling 
Assessment  of  Benefits  of  Partial  Source  Removal. 

This  project  seeks  to  develop  an  understanding  of  the  relationship  between  source  zone 
remediation  and  dissolved  plume  behavior  by  balancing  mass  removal  with  plume 
attenuation  for  the  purpose  of  optimizing  the  remedial  process.  The  objective  being  the 
development  of  cost-effective  remediation  endpoints  based  on  mass  flux. 

Details  on  this  project  can  be  obtained  at: 
http://www.serdp.org/research/CU/CU-1295.pdf 
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Minimum  System  Requirements 

The  Mass  Flux  Toolkit  model  requires  a  computer  system  capable  of  running  Microsoft®  Excel 
(2000/XP)  for  Windows  (2000/XP).  Operation  requires  an  IBM-compatible  PC  equipped  with  a 
Pentium  or  later  processor  running  at  a  minimum  of  450  MHz.  A  minimum  of  256  MB  of  system 
memory  (RAM)  is  strongly  recommended.  Computers  not  meeting  these  recommendations  will 
experience  slow  running  times  and/or  problems  with  memory. 

Installation  and  Start-Up 

The  software  is  installed  by  unzipping  the  Toolkit  model  file  (MassFluxToolkit.zip)  and  keeping  all 
the  unzipped  files  in  the  same  folder  on  your  computer  hard  drive.  To  use  the  software,  start 
Excel  and  load  the  MassFluxToolkit.xls  model  file  from  the  File  /  Open  menu.  If  you  are  using 
Excel  2000,  you  may  see  a  message  box  that  asks  you  whether  you  want  to  disable  or  enable  the 
macros.  For  the  Toolkit  to  operate  effectively,  you  must  enable  the  macros. 

NOTE:  Although  the  Transect  Calculator’  module  in  the  Toolkit  uses  Microsoft  Excel,  the  ‘Learn 
About  Flux’  module  calls  Microsoft  Word  and  Adobe  Acrobat  pdf  documents.  Some  features  in 
the  “Learn  about  Flux”  module  may  not  work  unless  you  have  these  programs  installed  on  your 
computer. 


Spreadsheet-Related  Problems 

Backspace  doesn’t  clear  cell.  Use  the  delete  key  on  the  keyboard  or  the  mouse  to  clear  data. 

The  buttons  won’t  work.  The  Mass  Flux  Toolkit  is  built  in  the  Excel  spreadsheet  environment, 
and  to  enter  data  one  must  click  anywhere  outside  the  cell  where  data  was  just  entered.  If  you 
can  see  the  numbers  you  just  entered  in  the  data  entry  part  of  Excel  above  the  spreadsheet,  the 
data  have  not  yet  been  entered.  Click  on  another  cell  to  enter  the  data. 

####  is  displayed  in  a  number  box.  The  cell  format  is  not  compatible  with  the  value,  (e.g.,  the 
number  is  too  big  to  fit  into  the  window).  To  fix  this,  unprotect  the  sheet.  Then,  select  the  cell, 
pull  down  the  format  menu,  select  Cells  and  click  on  the  Number  tab.  Change  the  format  of  the 
cell  until  the  value  is  visible.  If  the  values  still  cannot  be  read,  select  the  format  menu,  select 
Cells  and  click  on  the  Font  tab.  Reduce  the  font  size  until  the  value  can  be  read. 

#DIV/0!  is  displayed  in  a  number  box.  The  most  common  cause  of  this  problem  is  that  some 
input  data  are  missing.  In  some  cases,  entering  a  zero  in  a  box  will  cause  this  problem.  Double 
check  to  make  certain  that  data  required  for  your  run  have  been  entered  in  all  of  the  input  cells. 

#VALUE!  is  displayed  in  a  number  box.  The  most  common  cause  of  this  problem  is  that  some 
input  data  are  missing.  Double  check  to  make  certain  that  data  required  for  your  run  have  been 
entered  in  all  of  the  input  cells  and  all  options  have  been  selected. 


Common  Error  Messages 

Unable  to  Load  Help  File:  The  most  common  error  message  encountered  with  Toolkit  is  the 
message  ‘Unable  to  Open  Help  File’  after  clicking  on  a  Help  button.  Depending  on  the  version  of 
Windows  you  are  using,  you  may  get  an  Excel  Dialog  Box,  a  Windows  Dialog  Box,  or  you  may 
see  Windows  Help  load  and  display  the  error.  This  problem  is  related  to  the  ease  with  which  the 
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Windows  Help  Engine  can  find  the  data  file,  MassFluxToolkit.HLP.  Here  are  some  suggestions 
(in  decreasing  order  of  preference)  for  helping  WinHelp  find  it: 

•  If  you  are  asked  to  find  the  requested  file,  do  so.  The  file  is  called  MassFluxToolkit.HLP, 
and  it  was  installed  in  the  same  directory/folder  as  the  Mass  Flux  Toolkit  model  file 
(MassFluxToolkit.xls). 

•  Use  the  File/Open  menus  from  within  Excel  instead  of  double-clicking  on  the  filename  or 
Program  Manager  icon  to  open  the  Mass  Flux  Toolkit  model  file.  This  sets  the  current 
directory  to  the  directory  containing  the  Excel  file  you  just  opened. 
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APPENDIX  A.1 .  PROBABILITY  DISTRIBUTIONS 


This  section  describes  in  greater  detail  the  probability  distributions  employed  in  the  LHS 
sampling.  The  Mass  Flux  Toolkit  offers  the  user  three  distribution  options:  normal,  lognormal, 
and  uniform. 


A.1.1  Normal  Distributions 


Normal  distributions  are  defined  by  the  density  function 


fix)  = 


la2 


-oo  <  x  <  oo 


where  a  is  the  standard  deviation  and  p  the  mean  of  the  distribution.  The  Toolkit  assumes  that 
the  values  entered  in  the  “Input  Data  and  Grid”  screen  are  the  means.  The  uncertainty  analysis 
requires  the  user  to  specify  a  a  as  a  percentage  of  the  mean. 


A.1. 2  Lognormal  Distributions 


A  lognormal  distribution  is  a  distribution  whose  logarithms  are  normally  distributed.  The 
lognormal  density  function  is 


i  ( (In (x)~n)2  ) 

fix)  = - =e  {  2a1  J  x,  a>0 

xo-42k 

where  o  is  the  standard  deviation  and  pthe  mean  of  the  underlying  normal  distribution. 

Lognormal  distributions  are  typically  specified  in  two  ways  throughout  literature  (Swiler  and  Wyss, 
2004).  One  way,  as  described  above,  is  to  use  the  mean  and  standard  deviation  of  the 
underlying  normal  distribution.  The  other  way  is  to  by  using  the  mean  of  the  lognormal 
distribution  (a)  and  a  term  called  the  “Error  Factor”.  For  a  lognormal  distribution,  the  error  factor 
is  the  ratio  of  the  95th  percentile  to  the  median,  or  equivalently,  the  ratio  of  the  median  to  the  5th 
percentile.  Therefore,  the  error  factor  represents  the  width  of  a  90%  confidence  interval  around 
the  median. 


In  terms  of  the  error  factor,  the  relationship  between  the  underlying  normal  distribution  and  the 
lognormal  distribution  can  be  described  by 


a  =  ln(error  factor)/1 .645 


and 
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A  =  ln  (a) — — 

where  a  is  the  mean  of  the  lognormal  distribution,  and  a  and  p  the  standard  deviation  and  mean 
of  the  underlying  normal  distribution,  respectively. 

The  Mass  Flux  Toolkit  describes  the  lognormal  distribution  using  the  error  factor. 


A.1.3  Uniform  Distributions 


A  uniform  distribution  is  specified  over  a  particular  interval  and  implies  that  all  the  points  within 
that  interval  have  equal  probability  of  occurring.  The  uniform  probability  distribution  function  is 

f(x)  =  — A  <  x  <  B 
B-A 

where  A  and  B  are  the  lower  and  upper  bounds,  respectively. 
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APPENDIX  A.2.  LATIN  HYPERCUBE  SAMPLING 


A.2.1  Background 

In  the  standard  Monte  Carlo  approach,  simple  random  sampling,  a  large  number  of  runs  (typically 
100  to  1000)  are  required  to  obtain  a  meaningful  probability  distribution  for  the  parameter.  For 
each  run  of  the  standard  approach,  a  random  number  is  generated  for  every  value  of 
concentration,  hydraulic  conductivity,  and  hydraulic  gradient  entered  by  the  user.  This  set  of 
random  inputs  is  then  used  to  create  and  interpolate  the  respective  grids  and  mass  flux 
calculated  from  those  grids.  Repeating  this  procedure  a  large  number  times  yields  a  probability 
distribution  of  the  mass  flux  from  which  statistical  characteristics  such  as  mean,  percentile,  and 
variance  can  be  obtained.  However,  the  standard  Monte  Carlo  technique  tends  to  be  time 
intensive  because  of  the  large  number  of  iterations  requires. 

Latin  hypercube  sampling  (LHS)  is  a  type  of  stratified  Monte  Carlo  sampling.  It  can  be  viewed  as 
a  method  for  controlling  the  upper  and  lower  ends  of  the  distributions  used  in  the  analysis  such 
that  all  portions  of  the  distribution  are  well  represented  (U.S.  EPA,  1997;  Iman  and  Helton,  1988; 
Gwo  et  al.,  1996;  Helton  and  Davis,  2002).  In  LHS,  the  given  range  of  each  variable  in  X-i,...,Xn  is 
divided  into  j  number  of  equal  non-overlapping  probability  intervals.  A  random  number  with 
respect  to  the  probability  density  is  then  obtained  for  X )  for  each  j  interval.  The  j  values  of  X-i  that 
are  obtained  are  then  randomly  paired  with  the  j  values  obtained  for  X2  and  so  on  until  all  the  j 
values  in  Xn  are  paired  (Iman  and  Helton,  1988).  LHS  ensures  that  a  minimal  number  of 
probability  intervals  are  required  to  represent  the  full  range  of  the  distribution  (Bergin  et  al.,  1999). 
Additionally,  LHS  requires  fewer  simulations  than  simple  random  sampling  to  achieve  the  same 
level  of  precision  (U.S.  EPA,  1997;  Helton  and  Davis,  2002),  hence  a  fraction  of  the  time  involved 
with  simple  random  sampling. 

A.2.2  Mass  Flux  Toolkit  Methodology 

For  normal  and  lognormal  distributions,  the  Mass  Flux  Toolkit  applies  LHS  in  evaluating  input 
uncertainty  using  the  following  process: 

1 .  The  normal  probability  distribution  described  by  the  users  input  options  (or  the  underlying 
normal  distribution  for  the  lognormal  distribution)  is  divided  into  20  intervals  of  equal 
probability. 

2.  A  random  number  is  generated  from  each  interval.  The  numbers  are  sorted  by  magnitude 
and  each  assigned  a  rank.  Let  each  set  of  20  random  numbers  represent  a  LHS  set. 

3.  Step  2  is  repeated  for  each  value  of  concentration,  hydraulic  conductivity/Darcy  velocity, 
and  gradient  (if  required)  provided  by  the  user. 

4.  Each  concentration  LHS  set  is  then  randomly  paired  with  each  Hydraulic 
Conductivity/Darcy  velocity  set  and  a  Spearman  rank  correlation  obtained.  This  process 
is  repeated  until  a  correlation  coefficient  of  <  0.2  is  obtained. 

5.  If  hydraulic  conductivity  is  used  to  calculate  the  mass  flux  then,  in  addition  to  the  above 
pairing,  correlation  coefficients  of  the  pairings  of  concentration  and  gradient  and 
conductivity  and  gradient  are  also  obtained.  The  pairings  are  continuously  randomized 
until  all  the  correlation  coefficients  are  <0.2. 

6.  For  example,  let  concentration  =  A,  conductivity/Darcy  velocity  =  B,  and  gradient  =  C 
then,  the  parings  of  A,  B,  and  C  are  randomized  until 

Correlation  (A,  B)  <  0.2 
Correlation  (A,  C)  <  0.2 
Correlation  (A,  B)  <  0.2 

7.  Correlation  coefficients  are  obtained  to  minimize  any  undesired  pairings. 
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8.  Each  of  the  randomly  paired  values  of  concentration,  conductivity/Darcy  velocity,  and 
gradient  is  then  used  to  create  a  grid. 

9.  Each  grid  is  interpolated  (if  required)  by  the  interpolation  scheme  previously  chosen  by 
the  user  and  mass  flux  calculated.  This  results  in  20  mass  flux  measurements. 

10.  Steps  1  through  9  are  repeated  10  times  for  the  entire  LHS  application. 

11.  The  minimum,  15th  percentile,  50th  percentile,  85th  percentile,  maximum,  and  variance  are 
the  calculated  from  the  200  mass  fluxes  calculated. 

For  uniform  distributions: 

1.  Twenty  random  numbers  are  generated  between  the  lower  and  upper  bounds  specified 
by  the  user. 

2.  The  numbers  are  sorted  by  magnitude  and  each  assigned  a  rank. 

3.  Step  2  is  repeated  for  each  value  of  concentration,  hydraulic  conductivity/Darcy  velocity, 
and  gradient  (if  required)  provided  by  the  user. 

4.  Steps  4  through  1 1  for  the  normal/lognormal  distribution  are  then  followed. 
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EXAMPLE  1 :  MASS  FLUX  CALCULATION 


The  following  example  is  obtained  from  the  API’s  “Groundwater  Remediation  Strategies  Tool” 
(Newell  ef  at.,  2003). 


NOTES: 

—  0.1 -  Concentration  contour  in  mg/L 

0  100  200 

Horizontal  Scale  (ft) 

Source  zone 
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Mass  Flux  Toolkit  Transect  Calculator  Input 


Data  Type 

Parameter 

Value 

Source  of  Data 

Transect  Data 

•  Transect  1  From  Source: 

193  (ft) 

•  Plume  maps 

Hydrogeology 

•  Hydraulic  Conductivity: 

•  Hydraulic  Gradient: 

0.032  (cm/sec) 

0.002  (cm/cm) 

•  Slug-tests  results 

•  Static  water  level 
measurements 

End  of  Transect 

•  Transect  1  From  Source 
(start  of  transect  assumed  to 
be  left  facing  downstream): 

90  (ft) 

•  Plume  maps 

Concentration 

•  MTBE  Concentration 

Distance 

of  Point 

•  Plume  top  based  on  top  of 
water  table 

•  Plume  bottom  based  on  top  of 
clay  layer  below  aquifer 

from  Start 

Sampling  Plume 

Plume 

Name 

of  Transect 

Interval  Top 

Bottom 

Cone 

(ft) 

(ft  bgs)  (ft  bgs) 

(ft  bgs) 

(mg/L) 

TRI-2 

10 

5-10 

5 

15 

2.3 

TRI-2 

10 

10-15 

5 

15 

0.47 

TRI-4 

27.5 

5-10 

5 

20 

19.7 

TRI-4 

27.5 

10-15 

5 

20 

7.2 

TRI-4 

27.5 

15-20 

5 

20 

0.34 

TRI-6 

45 

5-10 

5 

20 

87.2 

TRI-6 

45 

10-15 

5 

20 

35.6 

TRI-6 

45 

15-20 

5 

20 

9.5 

TRI-8 

62.5 

5-10 

5 

20 

54.1 

TRI-8 

62.5 

10-15 

5 

20 

15.3 

TRI-8 

62.5 

15-20 

5 

20 

0.67 

TRI-12 

80 

5-10 

5 

15 

4.5 

TRI-12 

80 

10-15 

5 

15 

5.6 

Mass  Flux  Toolkit  Transect  Calculator  Output 


Parameter 

Value 

Source  of  Data 

Total  Mass  Flux  for  Transect  1 

105  g/day 

See  Figure  1-5 

Mass  Flux  Toolkit  Modeling  Summary: 

•  The  Mass  Flux  Toolkit  was  used  to  estimate  the  total  mass  flux  and  uncertainties 
associated  with  that  estimate  for  MTBE  at  a  petroleum  refinery. 

•  The  Mass  Flux  Transect  Calculator  was  used  as  the  primary  model  to  predict  the  total 
mass  flux  because  of  the  availability  of  site  hydrogeological  and  transect  concentration 
data. 
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•  The  Distance  of  Transect  from  Source  was  entered  in  Section  6.1,  hydraulic  conductivity 
in  Section  6.4,  hydraulic  gradient  in  Section  6.5,  and  individual  transect  data  in  the  table. 
Uniform  hydraulic  conductivity  and  gradient  fields  were  chosen  across  the  entire  transect. 
Because  the  data  was  collected  at  various  intervals,  the  Sampling  Interval  option  was 
selected  in  Section  6.6. 

•  A  grid  refinement  of  1  x  1  was  used  in  Section  7  as  a  conservative  measure. 

•  No  interpolation  was  required  for  the  concentration  grid  and  since  uniform  hydraulic 
conductivity  and  gradient  were  used,  no  interpolation  was  required  for  these  grids  either. 

•  Figure  1-5  shows  the  total  MTBE  mass  flux  calculated  across  Transect  1 . 

•  A  simple  uncertainty  analysis  was  performed  on  the  data  to  determine  potential  areas  of 
high  uncertainty.  Figure  1-6  shows  the  results  of  the  uncertainty  analysis. 

KEY  POINT: 


The  Mass  Flux  Toolkit  indicates  that  the  total  mass  groundwater  flux  across  transect  1  is  105 
g/day.  An  uncertainty  analysis  on  the  input  data  shows  that  of  all  the  sampling  points  collected, 
TRI-6  sampling  interval  15-20  ft  has  the  largest  influence  on  the  total  mass  flux  (23  %  contribution 
to  total  mass  flux). 
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Input  Data  and  Grid 


Site  Location  and  I.D.: 
Description: 


Data  Input  Instructions 

Enter  value  directly. 

Value  calculated  by  model 
(Don't  enter  any  data) 


4.  CHOOSE  TRANSECT  |  Transect  1  _▼] 

6.  ENTER  TRANSECT  DATA 

G.1  Distance  of  Transect  1  from  Source  - 


5.  CHOOSE  TIME  PERIOD  n~3 


I93 \fw 


6.2 


Q  Darcy  Velocity 


9i  Hydraul  i  c  Conducti  vi  ty 


6.6 


®  Sampling  Interval 


C  Mid  Poi nt  of  Sampl i ng  I nterval 


6. 3  Hydraul  i  c  Conducti  vi  ty  Uni  ts 

6. 4  Uni  form  Hydraul  i  c  Conducti  vi  ty? 

6. 5  Uni  form  Hydraul  i  c  G  radi  ent? 


Hydraul  i  c  Conducti  vi  ty 
->•  Hydraul  i  c  G  radi  ent 


Monitoring  Point 

Distance  of 
Monitoring 
Point  from  Start 

of  Transect 

(ft) 

Sampling  Interval 
(ft  bgs) 

Plume  Top 
(ft  bgs) 

Plume  Bottom 
(ft  bgs) 

Concentration  (mgfL) 

Constituent  A 

Constituent  B 

Top 

Bottom 

MTBE 

Start  of  Transect 

0 

0 

0 

End  of  Transect 

90 

0 

0 

TR1-2 

10 

5 

10 

5 

15 

2.3 

TR1-2 

10 

10 

15 

5 

15 

0.47 

TR1-4 

27.5 

5 

10 

5 

20 

19.7 

TR1-4 

27.5 

10 

15 

5 

20 

7.2 

TR1-4 

27.5 

15 

20 

5 

20 

0.34 

TR1-6 

45 

5 

10 

5 

20 

87.2 

TR1-6 

45 

10 

15 

5 

20 

35.6 

TR1-6 

45 

15 

20 

5 

20 

9.5 

TR1-8 

62.5 

5 

10 

5 

20 

54.1 

TR1-8 

62.5 

10 

15 

5 

20 

15.3 

TR1-8 

62.5 

15 

20 

5 

20 

0.67 

TR1-12 

80 

5 

10 

5 

15 

4.5 

TR1-12 

80 

10 

15 

5 

15 

5.6 

IOOSE  GRID  (OPTIONAL) 

Current  Grid 

Refine  Grid  By 

Refined  Grid 

8.  SELECT  CONSTITUENT  FOR  CALCULATIONS 

Number  of  rows 

10 

/ 

10 

$  MTBE 

O  Constituent  B 

Number  of  columns] 

7 

/ 

7 

1=1 


Next  Step: 

Back  to  Transect  Calculator  Screen 

import  MW  Dots  j  Export  MW  Dots 

See  Conc/Fiux  Br/ds  | 

Continue  Data  Input 

C/eer  Screen  1  Bests  Exempt 'e 

Besfore  Tebte  Formatting  Print 

HELP 

Figure  1-1.  Input  Data  and  Grid  Screen. 
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Data  Input  Instructions 

Enter  value  directly. 

1 10.80 

^  (Values  in 
v  italics  represent 

interpolated  values.) 

m 

.  T op/Bottom  of  pi  ume. 

v  (Don't  enter  any 

data). 

Grid  Completion:  Concentration 

Transect  1  Time  Period  1.  No  Interpolation  Required  (No  Space  Available  in  Grid  for  Interpolation). 


Step  10:  Interpola 

e  Cone  Horizontally 

Export  Cone  Grid 

Import  Cone  Grid 

_ 

Back  to  Data  Input 

Back  to  Cone  Grid  1 

LsJ 

HELP 

Start  of  Transect 
0.0 


£ 

a 

a> 

G 


5.0 

0.5 

0.0 

05 

no 

12.5 

no 

75.5 
77.0 

70.5 
20.0 


Next  Step: 
View  Hydraulic 
Conductivity  Grid 


MTBE  Concentrations  (mg/L) 

Distance  from  Edge  of  Transect  (ft) 


TR1-2 

70.0 


TR1-4 

27.5 


TR1-6 

45.0 


TR1-8 

02.5 


TR1-12 

00.0 


End  of  Transect 


00.0 


2.30E+00 


2.30E+00 


2.30E  +  00 


4.70E-01 


4.70E-01 


4.70E-01 


1.97E+01 


1.97E+01 


1.97E+01 


7.20E+00 


7.20E+00 


7.20E+00 


3.40E-01 


3.40E-01 


3.40E-01 


8.72E+01 


8.72E+01 


8.72E+01 


3.56E+01 


3.56E+01 


3.56E+01 


3.56E+01 


9.50E  +  00 


9.50E  +  00 


9.50E  +  00 


5.41E+01 


5.41E+01 


5.41E+01 


1.53E+01 


1.53E+01 


1.53E+01 


4.50E+00 


4.50E+00 


4.50E+00 


5.G0E+00 


5.G0E+00 


5.G0E+00 


G.70E-01 


G.70E-01 


G.70E-01 


Figure  1-2.  Concentration  Grid 
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Grid  Completion:  Hydraulic  Conductivity 

Transect  1  Time  Period  1.  Constant  Conductivity:  No  Interpolation  Required 


Hydraulic  Conductivity  (cm/sec) 


Data  Input  Instructions 

Enter  value  directly. 

10.80 

^  (Values  in 
r  italics  represent 

interpolated  values.) 

■i 

.  T op/Bottom  of  pi  ume. 

v  (Don't  enter  any 

data). 

Distance  from  Edge  of  Transect  (ft) 


Start  of  Transect 
0.0 


TR1-2 

mo 


TR1-4 

27.5 


TRIG 

45.0 


TR1-8 

62.5 


TR1-12 

00.0 


End  of  Transect 


90.0 


f 

o. 


5.0 


6.5 

0.0 


9.5 

no 


{2.5 


3.20E-Q2 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


3.20E-02 


Figure  1-3.  Hydraulic  Conductivity  Grid 
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Grid  Completion:  Hydraulic  Gradient 

Transect  1  Time  Period  1.  Constant  Gradient:  No  interpolation  Required 


Step  10:  Interpo 

late  i  Horizontally 

Export  i  Grid 

Import  i  Grid 

Next  Step: 

Back  to  Data  Input 

Back  to  Cone  Grid 

Back  to  K  Grid 

Print  |  HELP 

Calculate  Mass  Flux 

Data  Input  Instructions 

Enter  value  directly. 

1 10.80 

.  (Values  in 

r  italics  represent 

interpolated  values.) 

m 

.  T op/Bottom  of  pi  ume. 

r  (Don't  enter  any 

data). 

Hydraulic  Gradient  (cm/cm) 

Distance  from  Edge  of  Transect  (ft) 


Start  of  Transect 
0.0 


TR1-2 

too 


TR1-4 

27.5 


TRIG 

450 


TR1-8 

02.5 


TR1-12 

000 


End  of  Transect 


90.0 


c 

£ 

CL 

0) 

O 


50 

0.5 

00 

9.5 

no 

12.5 
74.0 

15.5 
17.0 

70.5 
20.0 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


2.00E-03 


Figure  1-4.  Hydraulic  Gradient  Grid 
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Mass  Flux  Result 


TOTAL  MASS  FLUX 


1.05E+02 


(g/day) 


3.85E+01 


(kg/yr) 


Next  Step:  Mass  1 

Run/View  Uncertainty  1 

View  Final  Concentration  Grid 

Flux  Summary 

Analysis  (Optional) 

Back  to  Data  Grid  j  Print  |  HELP  |  | 

Data  Representation 

1.  Bold  values  represent  calculations  based  on  given  values. 

2.  Values  in  italics  represent  calculations  based  on  interpolation. 

3.  Black  shaded  cells  represent  the  top  and  bottom  of  the  plume. 


SELECT  TRANSECT  TO  VIEW  |  Transect  1  _▼] 

SELECT  TIME  PERIOD  TO  VIEW  |l  _+J 
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Figure  1-5.  Mass  Flux  Results  Grid 
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EXAMPLE  1:  MASS  FLUX  CALCULATION 


Figure  1-6.  Simple  Uncertainty  Analysis 
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EXAMPLE  1:  MASS  FLUX  CALCULATION 


Transect  Calculator:  Mass  Flux  Summary 


Site  Location  and  I  D.: 
Description: 


Mass  Flux  Summary  for  MTBE 


Distance  from 

Mass  Flux  (gfday) 

Source  (ft) 

Time  Period  1 

Time  Period  2 

Time  Period  3 

Time  Period  4 

Time  Period  5 

Transect  1 

Transect  2 

Transect  3 

Transect  4 

Transect  5 

193 

1.05E  +  02 

85th  percentile  from  Advanced  Uncertainty  Analysis 
Maximum  of  the  three  Toolkit  interpolation  schemes 

Value  using  selected  interpolation  schemes 
Minimum  of  the  three  Toolkit  interpolation  schemes 
15th  percentile  from  Advanced  Uncertainty  Analysis 


-♦—Time  Period  1  | 

♦ 

loo.o  mo 

Distance  of  Transect  From  Source  (ft) 


Next  Step: 
Save  File 


View  Flux  Result 

View  Cone  Grid 

Print 

r  HELP 

Return  To  Toolkit  Main  Screen 

Return  To  Transect  Main  Screen 

Figure  1-7.  Mass  Flux  Summary 
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EXAMPLE  2:  MTBE  CONCENTRATION  AT  CONTROL  POINT 


EXAMPLE  2:  MTBE  CONCENTRATION  AT  CONTROL 
POINT 

The  following  example  is  obtained  from  the  API’s  “Groundwater  Remediation  Strategies  Tool” 
(Newell  et  al.,  2003). 

Total  Mass  Flux  is  Calculated  Using  Data  Provided  in  Example  1.  Transect  1  was  used  as  a 
conservative  pick  for  the  calculation  of  the  control  point  concentration  and  does  not  assume  any 
natural  attenuation. 


Mass  Flux  Toolkit  Receptor  Impact  -  Wells  Input 


Data  Type 

Parameter 

Value 

Supply  Well 

•  Flow  rate: 

10  (gpm) 

Mass  Flux  Toolkit  Receptor  Impact  -  Wells  Output 


Parameter 

Value 

Source  of  Data 

•  Total  Mass  Flux  for  Transect  1 

•  Concentration  at  Control  Point 

•  105  g/day 

•  1.94  mg/L 

See  Figure  1-5  of  Example  1 

See  Figure  2-1 
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EXAMPLE  2:  MTBE  CONCENTRATION  AT  CONTROL  POINT 


Mass  Flux  Toolkit  Modeling  Summary: 

•  The  Mass  Flux  Toolkit  was  used  to  estimate  the  total  mass  flux  and  uncertainties 
associated  with  that  estimate  for  MTBE  at  a  petroleum  refinery  (See  Example  1). 

•  The  Mass  Flux  Transect  Calculator  was  used  as  the  primary  model  to  predict  the  total 
mass  flux  because  of  the  availability  of  site  hydrogeological  and  transect  concentration 
data  (See  Example  1). 

•  The  Distance  of  Transect  from  Source  was  entered  in  Section  6.1 ,  hydraulic  conductivity 
in  Section  6.4,  hydraulic  gradient  in  Section  6.5,  and  individual  transect  data  in  the  table. 
Uniform  hydraulic  conductivity  and  gradient  fields  were  chosen  across  the  entire  transect. 
Because  the  data  was  collected  at  various  intervals,  the  Sampling  Interval  option  was 
selected  in  Section  6.6  (See  Example  1). 

•  A  grid  refinement  of  1  x  1  was  used  in  Section  7  (See  Example  1). 

•  No  interpolation  was  required  for  the  concentration  grid  and  since  uniform  hydraulic 
conductivity  and  gradient  were  used,  no  interpolation  was  required  for  these  grids  either 
(See  Example  1). 

•  Figure  1-5  shows  the  total  MTBE  mass  flux  calculated  across  Transect  1  (See  Example 
1)- 

•  The  Impact  to  Receptor-Wells  model  was  used  to  calculate  the  concentration  of  MTBE  at 
the  control  point. 

•  Transect  1  and  Time  Period  1  were  selected  in  Section  1  of  the  Impact  to  Receptor-Wells 
input  screen.  The  flow  rate  of  the  supply  well  was  also  entered  in  Section  1 . 

KEY  POINT: 


The  Mass  Flux  Toolkit  indicates  that  the  total  groundwater  mass  flux  across  transect  1  is  105 
g/day  which  results  in  a  concentration  of  1.94  mg/L  at  the  control  point.  This  value  assumes  the 
plume  is  at  a  steady-state  condition.  Some  corrective  action  measures  may  be  required  to 
reduce  the  control  point  concentration. 
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EXAMPLE  2:  MTBE  CONCENTRATION  AT  CONTROL  POINT 


Receptor  Impact  Worksheet  -  Wells  Version  1.0 Be to 


Site  Location  and  I.D.: 
Description: 


Texas 

MTBE 


Data  input  Instructions: 

1 10.80  | - ►  E/7ter  value  directly. 


10.80 


Value  ca/ctdafedb{/rj7ocfe/. 
(Don  'f  etifer  any  detaj 


1.  CONSTITUENT  CONCENTRATION 


9.  MTBE  C  Constituent  B 


Transect  of  Interest 
Time  Period  of  Interest 

Discharge  Rate  for  Each  Supply  Well 

Total  Mass  Flux  (Calculated  From  Mass  Flux  Transect  Method) 
Maximum  Constituent  Concentration  in  Water 
Extracted  from  Supply  Well 


2.  CAPTURE  ZONE  (OPTIONAL) 

Number  of  Supply  Wells 
Aquifer  Thickness 

Darcy  Velocity  for  Transect 


j  Transect  1 


h 


1.00E+Qlj|  gpm  ▼  | 


1.05E+02 


fg/diyl 


1.94E+00 


fmgLJ 


6.40E-05 


m 

fcmteec/ 


Capture  Zone  of  Suplly  Well(s) 


\ffl 


Return  To  Main  Screen 

Clear  Screen 

Paste  Example 

Save  File 

L  cad  File 

Print 

HELP 

Figure  2-1.  Control  Point  Concentration. 
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EXAMPLE  2:  MTBE  CONCENTRATION  AT  CONTROL  POINT 


CASE  STUDY 
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CASE  STUDY:  MTBE  AND  TBA 


Data  for  the  following  case  study  was  provided  by  LFR  Levine-Fricke  (personal  communication 
with  T.  L.  Roth  and  E.M.  Nichols  at  LFR  Levine-Fricke,  February  2006).  They  compared  methyl 
tert-butyl  ether  (MTBE)  and  tertiary  butyl  alcohol  (TBA)  mass  flux  calculated  using  the  Toolkit  to  a 
manual  calculation  based  on  a  site-specific  grid  using  an  interpolation  of  concentration  and  a 
detailed  sequential  indicator  kriging  interpolation  of  hydraulic  conductivity. 

This  case  study  documents  the  comparison  of  mass  flux  calculations.  Note  that  the  following  site 
figure  is  for  demonstration  only  and  does  not  represent  the  actual  site. 
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CASE  STUDY:  MTBE  AND  TBA 


Mass  Flux  Transect  Calculator  Input 


Data  Type 

Parameter 

Value 

Source  of  Data 

Transect  Data 

•  Transect  1  From  Source: 

1700  (ft) 

•  Plume  map 

Hydrogeology 

•  Hydraulic  Conductivity: 

See  table  below 

•  Slug-tests 

•  Hydraulic  Gradient: 

See  table  below 

results/estimations 

•  Static  water  level 

measurements 

End  of  Transect 

•  Transect  1  From  Source  (start  of  transect  assumed  to  be  at  the 

1700  (ft) 

•  Plume  map 

right  facing  downstream) 

Data 

Name 

Distance 

of  Point 

from  Start 

of  Transect 

Sampling 

Interval 

Plume 

Top 

Plume 

Bottom 

K 

i 

MTBE 

Cone 

TBA 

Cone 

•  Plume  top  based 
on  top  of  water 
table 

•  Plume  bottom 
based  on  top  of 
clay  layer  below 

(ft) 

(ft  bgs) 

(ft  bgs) 

(ft  bgs) 

(ft/d) 

(ft/ft) 

(mg/L) 

(mg/L) 

aquifer 

MW-1 

100 

5-20 

7.14 

81.30 

150 

0.00071 

0.5 

10 

MW-1 

100 

33-38 

7.14 

81.30 

22 

0.00071 

0.5 

10 

MW-1 

100 

51-56 

7.14 

81.30 

27 

0.00071 

0.5 

10 

MW-1 

100 

66-71 

7.14 

81.30 

1 

0.00071 

0.5 

10 

MW-2 

325 

7-22 

11.23 

92.00 

150 

0.0021 

2.8 

55 

MW-2 

325 

40-45 

11.23 

92.00 

26 

0.0021 

86 

310 

MW-2 

325 

62-67 

11.23 

92.00 

114 

0.0021 

7.2 

890 

MW-2 

325 

77-82 

11.23 

92.00 

1 

0.0021 

0.5 

10 

MW-3 

500 

10-25 

15.81 

92.00 

31 

0.0021 

21 

10 

MW-3 

500 

38-43 

15.81 

92.00 

66 

0.0021 

14 

10 

MW-3 

500 

58-63 

15.81 

92.00 

150 

0.0021 

17 

10 

MW-3 

500 

77-82 

15.81 

92.00 

1 

0.0021 

4.6 

30 

MW-4 

665 

15-30 

21.15 

95.00 

21 

0.0021 

11 

0.01 

MW-4 

665 

44-49 

21.15 

95.00 

50 

0.0021 

3.7 

0.01 

MW-4 

665 

63-68 

21.15 

95.00 

87 

0.0021 

3 

0.01 

MW-4 

665 

80-85 

21.15 

95.00 

1 

0.0021 

0.5 

0.01 
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CASE  STUDY:  MTBE  AND  TBA 


Distance 

of  Point 

from  Start 

Name  of  Transect 

Sampling 

Interval 

Plume 

Top 

Plume 

Bottom 

K 

i 

MTBE 

Cone 

TBA 

Cone 

(ft) 

(ft  bgs) 

(ft  bgs) 

(ft  bgs) 

(ft/d) 

(ft/ft) 

(mg/L) 

(mg/L) 

MW-5 

805 

20-35 

24.91 

96 

33 

0.0021 

12 

0.01 

MW-5 

805 

47-52 

24.91 

96 

57 

0.0021 

11 

0.01 

MW-5 

805 

65-70 

24.91 

96 

150 

0.0021 

4.2 

0.01 

MW-5 

805 

81-86 

24.91 

96 

1 

0.0021 

0.5 

0.01 

MW-6 

1005 

24-39 

29.55 

77 

150 

0.0019 

6.2 

0.01 

MW-6 

1005 

47-52 

29.55 

77 

58 

0.0019 

5.2 

0.01 

MW-6 

1005 

62-67 

29.55 

77 

1 

0.0019 

0.5 

0.01 

MW-7 

1185 

30-45 

35.63 

82 

44 

0.0019 

6.9 

0.01 

MW-7 

1185 

52-57 

35.63 

82 

51 

0.0019 

7.4 

0.01 

MW-7 

1185 

67-72 

35.63 

82 

1 

0.0019 

0.5 

0.01 

MW-8 

1430 

29-44 

32.54 

81 

44 

0.0019 

4.2 

0.01 

MW-8 

1430 

51-56 

32.54 

81 

51 

0.0019 

6.3 

0.01 

MW-8 

1430 

66-71 

32.54 

81 

1 

0.0019 

0.5 

0.01 

Mass  Flux  Transect  Calculator  Output 


Data  Type 

Parameter 

Value 

LFR  Levine-Fricke 
Independent 

Estimate 

Total  Mass  Flux  for  Transect  1 

MTBE  TBA 

1.7g/day  9.4  g/day 

Toolkit  Estimates 

Total  Mass  Flux  for  Transect  1 

MTBE:  Figures  3-8  and  3-9;  3-19  and  3-20;  and  3-30  and  3-31 

TBA:  Figures  3-37  and  3-38;  3-44  and  3-45;  and  3-51  and  3-52 

MTBE  TBA 

Natural  Neighbor  Interpolation  2.5  g/day  21.2  g/day 

Linear  Interpolation  2.9  g/day  20.4  g/day 

Log  Transformation  1.7  g/day  9.3  g/day 
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CASE  STUDY:  MTBE  AND  TBA 


Mass  Flux  Toolkit  Modeling  Summary: 

•  The  Mass  Flux  Toolkit  was  used  to  estimate  the  total  mass  flux  and  uncertainties 
associated  with  that  estimate  for  MTBE  and  TBA  at  a  petroleum  refinery. 

•  The  Mass  Flux  Transect  Calculator  was  used  as  the  primary  model  to  predict  the  total 
mass  flux  because  of  the  availability  of  site  hydrogeological  and  transect  concentration 
data. 

•  The  Distance  of  Transect  from  Source  was  entered  in  Section  6.1 ,  hydraulic  conductivity 
in  Section  6.4,  hydraulic  gradient  in  Section  6.5,  and  individual  transect  data  in  the  table. 
Because  the  data  was  collected  at  various  intervals,  the  Sampling  Interval  option  was 
selected  in  Section  6.6. 

•  A  grid  refinement  of  2  x  2  (to  match  the  LFR  Levine-Fricke  grid)  was  used  in  Section  7. 

•  Mass  flux  was  calculated  using  each  of  the  three  interpolation  techniques  available  in  the 
Toolkit  for  comparison  purposes. 

•  Figure  3-1  shows  the  input  data  for  the  Toolkit. 

•  Figures  3-2  and  3-3,  3-13  and  3-14,  and  3-24  and  3-25  show  the  final  MTBE 
concentration  girds  calculated  using  the  nearest  neighbor  interpolation,  linear 
interpolation,  and  log  transformation,  respectively. 

•  Figures  3-4  and  3-5,  3-15  and  3-16,  and  3-26  and  3-27  show  the  final  hydraulic 
conductivity  girds  calculated  using  the  nearest  neighbor  interpolation,  linear  interpolation, 
and  log  transformation,  respectively. 

•  Figures  3-6  and  3-7,  3-17  and  3-18,  and  3-28  and  3-29  show  the  final  hydraulic  gradient 
girds  calculated  using  the  nearest  neighbor  interpolation,  linear  interpolation,  and  log 
transformation,  respectively. 

•  Figures  3-8  and  3-9,  3-19  and  3-20,  and  3-30  and  3-31  show  the  total  MTBE  mass  flux 
calculated  using  the  nearest  neighbor  interpolation,  linear  interpolation,  and  log 
transformation,  respectively. 

•  A  simple  uncertainty  analysis  was  performed  on  the  data  to  determine  potential  areas  of 
high  uncertainty.  Figures  3-10  and  3-11,  3-21  and  3-22,  and  3-32  and  3-33  show  the 
results  of  the  uncertainty  analysis  for  MTBE  for  the  nearest  neighbor  interpolation,  linear 
interpolation,  and  log  transformation,  respectively. 

•  Figures  3-12,  3-23,  and  3-34  show  the  MTBE  Mass  Flux  Summary  results  for  the  nearest 
neighbor  interpolation,  linear  interpolation,  and  log  transformation,  respectively. 

•  Figures  3-35  and  3-36,  3-42  and  3-43,  and  3-49  and  3-50  show  the  final  TBA 
concentration  girds  calculated  using  the  nearest  neighbor  interpolation,  linear 
interpolation,  and  log  transformation,  respectively. 

•  Figures  3-37  and  3-38,  3-44  and  3-45,  and  3-51  and  3-52  show  the  total  TBA  mass  flux 
calculated  using  the  nearest  neighbor  interpolation,  linear  interpolation,  and  log 
transformation,  respectively. 
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CASE  STUDY:  MTBE  AND  TBA 


•  Figures  3-39  and  3-40,  3-46  and  3-47,  and  3-53  and  3-54  show  the  results  of  the 
uncertainty  analysis  for  TBA  for  the  nearest  neighbor  interpolation,  linear  interpolation, 
and  log  transformation,  respectively. 

•  Figures  3-41,  3-48,  and  3-55  show  the  TBA  Mass  Flux  Summary  results  for  the  nearest 
neighbor  interpolation,  linear  interpolation,  and  log  transformation,  respectively. 

KEY  POINTS: 


The  total  mass  groundwater  flux  across  transect  1  obtained  using  the  log  transformation  for  all 
concentration,  conductivity,  and  gradient  grids  of  1.7  g/day  for  MTBE  and  9.3  g/day  for  TBA  was 
in  agreement  with  the  LFR  Levine-Fricke  manual  calculation  of  1 .7  g/day  for  MTBE  and  9.4  g/day 
for  TBA. 

For  MTBE  log  transformation,  an  uncertainty  analysis  on  the  input  data  shows  that  of  all  the 
sampling  points  collected,  concentration  at  MW-2  sampling  interval  7-22  ft  has  the  largest 
influence  on  the  total  mass  flux  (87  %  contribution  to  total  mass  flux).  Hydraulic  conductivity  at 
MW-2  sampling  interval  40-45  ft  has  the  second  largest  influence  on  the  total  mass  flux  (68  % 
contribution  to  total  mass  flux). 

For  the  MTBE  nearest  neighbor  interpolation,  hydraulic  conductivity  at  MW-2  sampling  interval 
40-45  ft  has  the  largest  influence  on  the  total  mass  flux  (62  %  contribution  to  total  mass  flux) 
while  concentration  at  MW-2  sampling  interval  7-22  ft  has  the  second  largest  influence  on  the 
total  mass  flux  (53  %  contribution  to  total  mass  flux). 

For  MTBE  linear  interpolation,  concentration  at  MW-2  sampling  interval  7-22  ft  has  the  largest 
influence  on  the  total  mass  flux  (65  %  contribution  to  total  mass  flux)  while  hydraulic  conductivity 
at  MW-2  sampling  interval  40-45  ft  has  the  second  largest  influence  on  the  total  mass  flux  (55  % 
contribution  to  total  mass  flux). 

For  TBA  log  transformation,  an  uncertainty  analysis  on  the  input  data  (log  transformation)  shows 
that  of  all  the  sampling  points  collected,  hydraulic  conductivity  at  MW-2  sampling  interval  40-45  ft 
has  the  largest  influence  on  the  total  mass  flux  (68  %  contribution  to  total  mass  flux). 
Concentration  at  MW-2  sampling  interval  62-67  ft  also  has  a  large  influence  on  the  total  mass  flux 
(56  %  contribution  to  total  mass  flux). 

For  TBA  nearest  neighbor  interpolation,  concentration  at  MW-2  sampling  interval  62-67  ft  has  the 
largest  influence  on  the  total  mass  flux  (69  %  contribution  to  total  mass  flux)  while  hydraulic 
conductivity  at  the  same  location  has  the  second  largest  influence  on  the  total  mass  flux  (55  % 
contribution  to  total  mass  flux). 

For  TBA  linear  interpolation,  concentration  at  MW-2  sampling  interval  62-67  ft  has  the  largest 
influence  on  the  total  mass  flux  (55  %  contribution  to  total  mass  flux)  while  hydraulic  conductivity 
at  the  same  location  has  the  second  largest  influence  on  the  total  mass  flux  (49  %  contribution  to 
total  mass  flux). 

In  summary,  different  methods  of  interpolation  can  result  in  different  total  mass  flux  calculations. 
In  this  case,  a  difference  of  a  factor  of  2  was  observed  for  TBA  calculations. 
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CASE  STUDY:  MTBE  AND  TBA 


Input  Data  and  Grid 


Site  Location  and  I.D.: 
Description: 

4  CHOOSE  TRANSECT 


LFR  Site 
MTBE  and  TBA 


5.  CHOOSE  TIME  PERIOD  n~3 


G.  ENTER  TRANSECT  DATA 

G.1  Distance  of  Transect  1  from  Source 

G.2  C  Darcy  Velocity  Hydraulic  Conductivity 


- ►  ll  mol //t/ 

G.G  $  Sampling  Interval  C  Mid  Point  of  Sampling  Interval 


G.  3  Hydraul  i  c  Conducti  vi  ty  Uni  ts 

G.4  Uniform  Hydraulic  Conductivity? 
G.  5  Uni  form  Hydraul  i  c  G  radi  ent? 


Data  Input  Instructions 

I  1  I  E  nter  val  ue  di  recti  y. 

Value  calculated  by  model 
(Don't  enter  any  data) 


Monitoring  Point 

Distance  of 
Monitoring 
Point  from  Start 

of  Transect 

(fu 

Sampling  Interval 
(ft  bgs) 

Plume  Top 
(ft  bgs) 

Plume  Bottom 
(ft  bgs) 

Hydraulic 

Conductivity 

(ft/d) 

Hydraulic 

Gradient 

(ftfft) 

Concentration  (ugfL) 

Constituent  A 

Constituent  B 

1  I°P  1 

Bottom 

MTBE 

TBA 

1 

Start  of  Transect 

0 

0 

0 

2 

End  of  Transect 

1700 

0 

0 

3 

MW-1 

100 

5 

20 

7.14 

81.3 

150 

0.00071 

0.5 

10 

4 

MW-1 

100 

33 

38 

7.14 

81.3 

22 

0.00071 

0.5 

10 

5 

MW-1 

100 

51 

56 

7.14 

81.3 

27 

0.00071 

0.5 

10 

6 

MW-1 

100 

66 

71 

7.14 

81.3 

1 

0.00071 

0.5 

10 

7 

MW-2 

325 

7 

22 

11.23 

92 

150 

0.002133 

2.8 

55 

8 

MW-2 

325 

40 

45 

11.23 

92 

26 

0.002133 

86 

310 

9 

MW-2 

325 

62 

67 

11.23 

92 

114 

0.002133 

7.2 

890 

10 

MW-2 

325 

77 

82 

11.23 

92 

1 

0.002133 

0.5 

10 

11 

MW-3 

500 

10 

25 

15.81 

92 

31 

0.002133 

21 

10 

12 

MW-3 

500 

38 

43 

15.81 

92 

66 

0.002133 

14 

10 

13 

MW-3 

500 

58 

63 

15.81 

92 

150 

0.002133 

17 

10 

14 

MW-3 

500 

77 

82 

15.81 

92 

1 

0.002133 

4.6 

30 

15 

MW-4 

G65 

15 

30 

21.15 

95 

21 

0.002133 

11 

0.01 

7.  CHOOSE  GRID  (OPTIONAL) 


Current  Grid 


Number  of  rows 
Number  of  columns 


10 

10 


Ref  i  ne  G  ri  d  By  Ref  i  ned  G  ri  d 


20 

19 

8.  SELECT  CONSTITUENT  FOR  CALCULATIONS 

®  MTBE  c  TBA 


Next  Step: 

Back  to  Transect  Calculator  Screen  I 

import  MW'  Data 

Export  MW'  Data 

1  See  ConctF/ux  Grids  1 

Continue  Data  Input  1 

Paste  Exempt e  |  Be  store  Teb/e  Formatting  | 

Print 

HELP 

Figure  3-1.  Input  Data  and  Grid  Screen. 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Concentration 

Transect  1  Time  Period  1.  Nearest  Neighbor  Interpolation 


||  Step  10:  Interpolate  Cone  Horizontally 

Export  Cone  Grid  1 

Import  Cone  Grid 

_ 

Back  to  Data  Input  Jj 

| 

Back  to  K  Grid 

[  Print 

HELP 

Next  Step: 

Interpolate  Hydraulic 
Conductivity  Grid 


Data  Input  Instructions 

Enter  value  directly. 

110.80  1 

^  (Values  in 

italics  represent 

interpolated  values.) 

Top/Bottom  of  plume. 

r  (Don't  enter  any 

data). 

MTBE  Concentrations  (ug/L) 

Distance  from  Edge  of  Transect  (ft) 


Start  of  Transect 
0.0 

50.0 

MW  -1 

100.0 

212.5 

MW-2 

325.0 

412.5 

MW -3 

500.0 

582.5 

MW -4 

665.0 

735.0 

MW -5 

805.0 

7.1 

0 

2.50E-01 

5.00E-0  1  •-  -.  ;  - . 

11.6 

0 

2.50E-01 

5.00E-01 

1.65E+00 

2.80E+00  u  ; 

16.0 

0 

2.50E-01 

5.00E-01 

1.65E+00 

2.80E+00 

1. 19E+01 

2.10E+01 

20.5 

0 

0.00E+00 

5.00E-01 

1.65E+00 

2.80  E+00 

2.10E+01 

2.10E+01 

24.9 

0 

0.00E+00 

5.00E-01 

1.27E+00 

2.80E+00 

2.10E+01 

2.10E+01 

1.10E+01 

1.10  E+01 

1.15E+01 

1.20  E+01 

29.4 

0 

0.00E+00 

5.00E-01 

5.00E-01 

8.60E+01 

4.03E+01 

1.75E+01 

1.10E+01 

1.10E+01 

1.20E+01 

1.20  E+01 

33.8 

0 

2.50E-01 

5.00E-01 

5.00E-01 

8.60E+01 

5.00E+01 

1.40  E+01 

1.40E+01 

7.35E+00 

1.20E+01 

1.20E+01 

38.2 

0 

0.00E+00 

5.00E-01 

8.60  E+01 

8. 60 E+01 

5.00E+01 

1.40  E+01 

1.40E+01 

3.70E+00 

3.70E+00 

1.15E+01 

42.7 

0 

0.00E+00 

5.00E-01 

8.60  E+01 

8. 60 E+01 

8.60E+01 

1.40E+01 

3.70E+00 

3.70  E+00 

3.70E+00 

1.10E+01 

47.1 

0 

0.00E+00 

5.00E-01 

4.33  E+01 

8.60  E+01 

8.60E+01 

1.55E+01 

3.70E+00 

3.70E+00 

1.10E+01 

1.10  E+01 

51.6 

0 

2.50E-01 

5.00E-01 

5.00E-01 

4.66  E+01 

1.70E+01 

1.70E+01 

1.70E+01 

3.70E+00 

1.10E+01 

1.10E+01 

56.0 

0 

0.00E+00 

5.00E-01 

3.85E+00 

7. 20  E+00 

1.70E+01 

1.70E+01 

1.70E+01 

3.00E+00 

6. 07 E+00 

7.60E+00 

60.5 

0 

0.00E+00 

5.00E-01 

7.20E+00 

7.20E+00 

1.21  E+01 

1.70E+01 

1.70E+01 

3.00E+00 

3.60E+00 

4.20E+00 

64.9 

0 

2.50E-01 

5.00E-01 

5.00E-01 

7. 20  E+00 

1.21  E+01 

1.70E+01 

3.00E+00 

3. 00  E+00 

3. 60  E+00 

4. 20  E+00 

69.3 

0 

0.00E+00 

5.00E-01 

5.00E-01 

3.85E+00 

7.33E+00 

1.08E+01 

3.00E+00 

3.00E+00 

3.60E+00 

4.20E+00 

73.8 

0 

0.00E+00 

5.00E-01 

5.00E-01 

5.00E-01 

2.55E+00 

4.60  E+00 

2.55E+00 

5.00E-01 

5.00E-01 

2.35E+00 

78.2 

0 

0.00E+00 

5.00E-01 

5.00E-01 

5.00E-01 

2.55E+00 

4.60E+00 

2.55E+00 

5.00E-01 

5.00E-01 

5.00E-01 

82.7 

,  5.00E-Q1 

5.00E-01 

2.55E+00 

4.60  E+00 

5.00E-01 

5.00E-01 

5.00E-01 

5.00E-01 

87.1 

5  £>  oe -o ; 

5.00E-01 

2.55E+00 

4.60  E+00 

5.00E-01 

5.00E-01 

5.00E-01 

5.00E-01 

91.6 

5.00E-01 

5.00E-01 

5.00E-01 

5.00E-01 

Figure  3-2.  MTBE  Nearest  Neighbor  Interpolation:  Concentration  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


MTBE  Concentrations  (ug/L) 

Distance  from  Edge  of  Transect  (ft) 

MW -6  MW -7  MW -8  End  of  Transect 

905.0  1005.0  1095.0  1185.0  1307.5  1430.0  1565.0  1700.0 


7.1 

11.6 

16.0 

20.5 

24.9 

3.  294 

|  33.8 

38.2 

S'  42.7 

a 

47.1 

51.6 
56.0 
6  0.5 

6 4.9 

69.3 
73.8 

78.2 

82.7 
87.1 
91.6 
96.0 


6.20E+00 


6.20E+00 


8.10E+00 


8.10E+00 


8.10E+00 


5.00E-01 


5.00E-01 


4.20E+00 


4.20E+00 


2.35E+00 


6. 20  E+00 


4.20E+00 


6.20E+00 


5.20E+00 


5.20  E+00 


5.20E+00 


5.00E-01 


5.00E-01 


5  00E-01 


5.00E-01 


5.00E-01 


6.90E+00 


6.90E+00 


5. 20E+00 


7.40E+00 


3.95E+00 


5.00E-01 


5.00E-01 


5.00E-01 


5.00E-01 


6.90E+00 


6.90E+00 


7.1 5E+00 


7.40E+00 


7. 40  E+00 


5.00E-01 


5.00E-01 


5.00E-01 


5.00E-01 


5.55E+00 


6.90E+00 


6.87  E+00 


6.85E+00 


6. 85  E+00 


5.00E-01 


5  00E-01 


5.00E-01 


5.00E401 


4.20E+00 


4.20E+00 


4.20E+00 


6. 30  E+00 


6.30E+00 


6.30E+00 


5.00E-01 


5.00E-01 


5.00E-01 


5.00E-01 


5.00E-01 


2.10E+00 


2.10E+00 


0.00E+00 


0  00E+00 


3.15E+00 


0.00E+00 


0.00E+00 


2.50E401 


0.00E+00 


0.00E+00 


0.00E+00 


Figure  3-3.  MTBE  Nearest  Neighbor  Interpolation:  Concentration  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Hydraulic  Conductivity 

Transect  1  Time  Period  1.  Nearest  Neighbor  interpolation 


0  Step  10:  Interpolate  K  Horizontally 

Export  K  Grid  1  Import  K  Grid 

. 

Back  to  Data  Input 

^  Back  to  Cone  Grid  J 

1 1  Print  | 

[  HELP 

Next  Step: 

Interpolate  Hydraulic 
Gradient  Grid 


3 


Data  Input  Instructions 


Enter  value  directly. 
(Values  in 
italics  represent 
interpolated  values.) 

Top/Bottom  of  plume 
(Dont  enter  any 
data). 


Hydraulic  Conductivity  (ft/d) 

Distance  from  Edge  of  Transect  (ft) 


Start  of  Transect 
0.0 

50.0 

MW-1 

100.0 

212.5 

mt-2 

325.0 

412.5 

MW-3 

500.0 

582.5 

MW -4 

665.0 

735.0 

MW -5 

805.0 

7.1 

1.50  E+02 

1.50E+02 

1.50E+02  .._v  L 

11.6 

1 .50  E+02 

1.50E+02 

1.50E+02 

1.50E+02 

1.50  E+02 

16.0 

1 .50  E+02 

1.50E+02 

1.50  E+02 

1.50E+02 

1.50  E+02 

9.05E+01 

3.10E+01 

20.5 

1.50  E+02 

1.50E+02 

1.50  E+02 

1.50E+02 

1.50E+02 

3.10E+01 

3.10E+01 

24.9 

1 .50  E+02 

1.50E+02 

8.60E+01 

1.07E+02 

1.50  E+02 

3.10E+01 

3.10E+01 

2.10E+01 

2.10E+01 

2.70E+01 

3.30  E+01 

29.4 

1.50  E+02 

1.50E+02 

2.20E+01 

2.20E+01 

2.60  E+01 

4.10E+01 

4.85E+01 

2.10E+01 

2.10E+01 

3.30E+01 

3.30E+01 

33.8 

1 .50  E+02 

8.60E+01 

2.20E+01 

2.20E+01 

2.60E+01 

4.60E+01 

6.60E+01 

6.60E+01 

3.55E+01 

3.30E+01 

3.30E+01 

38.2 

1.50  E+02 

1.50E+02 

2.20E+01 

2.60E+01 

2.60E+01 

4.60E+01 

6.60E+01 

6.60E+01 

5.00E+01 

5.00E+01 

4.50E+01 

42.7 

1.50  E+02 

1.50E+02 

2.45E+01 

2.60E+01 

2.60E+01 

2.60E+01 

6.60E+01 

5.00E+01 

5.00E+01 

5.00E+01 

5. 70E+01 

47.1 

1.50  E+02 

1.50E+02 

2.70E+01 

2.65E+01 

2.60E+01 

2.60E+01 

1.08E+02 

5.00E+01 

5.00E+01 

5.70E+01 

5. 70  E+01 

51.6 

1.50  E+02 

8.85E+01 

2.70E+01 

2.70E+01 

7.00E+01 

1.50E+02 

1.50  E+02 

1.50E+02 

5.00E+01 

5.70E+01 

5.70E+01 

56.0 

1.50  E+02 

1.50E+02 

2.70E+01 

7.05E+01 

1.14E+02 

1.50E+02 

1.50  E+02 

1.50E+02 

8. 70E+01 

9.80E+01 

1.04E+02 

60.5 

1.50  E+02 

1.50E+02 

1.00E+00 

1.14E+02 

1.14E+02 

1.32E+02 

1.50  E+02 

1.50E+02 

8. 70E+01 

1.19E+02 

1.50E+02 

64.9 

1.50  E+02 

7.55E+01 

1.00E+00 

1.00E+00 

1.14E+02 

1.32E+02 

1.50  E+02 

8.70E+01 

8.70E+01 

1.1 9  E+02 

1.50  E+02 

69.3 

1.50  E+02 

1.50E+02 

1.00E+00 

1.00E+00 

5.75E+01 

6.65E+01 

7.55E+01 

8.70E+01 

8.70E+01 

1.19E+02 

150E+02 

73.8 

1.50  E+02 

1.50E+02 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

7.55E+01 

78.2 

1.50  E+02 

1.50E+02 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

82.7 

V  7  00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

100E+00 

1.00E+00 

1.00E+00 

1.00E+00 

87.1 

o  - + 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

91.6 

1. 00E+00 

1.00E+00 

1.00E+00 

1.00E+00 

Figure  3-4.  Nearest  Neighbor  Interpolation:  Hydraulic  Conductivity  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


w> 

U) 

-Q 

C 


Q 


Hydraulic  Conductivity  (ft/d) 

Distance  from  Edge  of  Transect  (ft) 

MW -6  MW-7  MW-8  End  of  Transect 

905.0  1005.0  1095.0  1185.0  1307.5  1430.0  1565.0  1700.0 


7.1 

11.6 

16.0 

20.5 
24.9 

29.4 

33.8 

38.2 

42.7 

47.1 

51.6 
56.0 

60.5 

64.9 

69.3 

73.8 

78.2 

82.7 
87.1 

91.6 
96.0 


1.50E+02 


1.50E+02 


1.50E+02 


4.40E+01 


1.50E+02 


4.40E+01 


4.40E+01 


4.40E+01 


4.40E+01 


4.40E+01 


2.25E+01 


2.25E+01 


1.00E+00 


1.00E+00 


5.75E+01 


5.80E+01 


4.40E+01 


4.40E+01 


4.40E+01 


4.40E+01 


1.00E+00 


1.00E+00 


5.75E+01 


5 .80  E+0 1 


5.80E+01 


4.75E+01 


4.97  E+01 


5.10E+01 


1.00E+00 


1.00E+00 


S.7SE+01 


5.80E+01 


5.10E+01 


5.10E+01 


5. 1 0E+01 


5.10E+01 


2. 60  E+01 


1.00E+00 


1.00E+00 


1.00E+00 


2.60E+01 


5.10E+01 


5.10E+01 


5.10E+01 


1.00E+00 


1.00E+00 


1.00E+00 


1 .00  E+00 


1.00E+00 


1.00  E+00 


1.00  E+00 


1.00  E+00 


1.00E+00 


1.00  E+00 


1.50E+02 


1.00E+00 


1.00E+00 


1.00E+00 


1.00E+00 


1.00E+00 


1.00  E+00 


1.00  E+00 


1.50E+02 


1.00E+00 


1.00E+00 


1.00  E+00 


1.00E+00 


1.00E+00 


1.00E+00 


1.00  E+00 


7.55E+01 


1.00E+00 


1.00E+00 


1.00E+00 


1.00  E+00 


1.00  E+00 


1.00E+00 


1.00  E+00 


1.00E+00 


1.00E+00 


1.00E+00 


Figure  3-5.  Nearest  Neighbor  Interpolation:  Hydraulic  Conductivity  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Hydraulic  Gradient 

Transect  1  Time  Period  1.  Nearest  Neighbor  Interpolation 


H  Step  10:  Interpolate  i  Horizontally 

Export  i  Grid 

Import  i  Grid 

Next  Step: 

Back  to  Data  Input  £  Back  to  Cone  Grid 

Back  to  K  Grid 

|(  Print  [j  HELP 

Calculate  Mass  Flux 

Data  Input  Instructions 


1 10.80  |  . 


Enter  value  directly. 
(Values  in 
italics  represent 
interpolated  values.) 

Top/Bottom  of  plume. 
(Dont  enter  any 
data). 


Hydraulic  Gradient  (ft/ft) 

Distance  from  Edge  of  Transect  (ft) 


Start  of  Transect 

MW  -1 

MW -2 

MW-3 

MW -4 

MW -5 

0.0 

50.0 

100.0 

212.5 

325.0 

412.5 

500.0 

582.5 

665.0 

735.0 

805.0 

7.1 

7.10E-04 

7.10E-04 

7.I0E-04  .  j 

11.6 

7.10E-04 

7.10E-04 

7.10E-04 

1.42E-03 

2.13E-03  ,  : 

16.0 

7.10E-04 

7.10E-04 

7.10E-04 

1.42E-03 

2.13E-03 

2. 13E-03 

2.13E-03 

20.5 

7.10E-04 

7.10E-04 

7.10E-04 

1.42E-03 

2.13E-03 

2.13E-03 

2.13E-03 

24.9 

7.10E-04 

7.10E-04 

7.10E-04 

1.18E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

29.4 

7.10E-04 

7.10E-04 

7.10E-04 

7.10E-04 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

33.8 

7.10E-04 

7.10E-04 

7.10E-04 

7. 1 0E-04 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

38.2 

7.10E-04 

7.10E-04 

7.10E-04 

2.13E-03 

2.13E-03 

2. 13E-03 

2.13E-03 

2. 13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

42.7 

7.10E-04 

7.10E-04 

7.10E-04 

2.13E-03 

2.13E-03 

2. 13E-03 

2.13E-03 

2. 1 3E-03 

2.13E-03 

2.13E-03 

2.13E-03 

47.1 

7.10E-04 

7.10E-04 

7.10E-04 

1.42E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

51.6 

7.10E-04 

7.10E-04 

7.10E-04 

7.10E-04 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 
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CASE  STUDY:  MTBE  AND  TBA 
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Data  Representation 

1.  Bold  values  represent  calculations  based  on  given  values. 

2.  Values  in  italics  represent  calculations  based  on  interpolation. 

3.  Black  shaded  cells  represent  the  top  and  bottom  of  the  plume. 
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CASE  STUDY:  MTBE  AND  TBA 


Interpolation  Error  Results 
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Next  Step: 
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Back  to  Mass  Flux  Result 
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Back  to  Data  Input 
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1.  Tte  table  below  shows  all  concEntration  and,  if  applicable, 
non-uniform  condjctivity/Darcy  velocity  and  gadient 
input  data. 

2.  During  the  uncertainty  analysis,  each  of  the  points  in  the 
table  are  removed  cne  at  a  time  and  the  mass  flux 
calculated  without  that  point  using  the  interpolation 
scheme  specified  in  the  Data  Input  section  (note  that 
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Figure  3-10.  MTBE  Nearest  Neighbor  Interpolation:  Simple  Uncertainty  Analysis  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


_ 

RESULTS 

End  of  Transect 

Top  of 
Sampling 
Interval 

Bottom  of 
Sampling 
Interval 

Param  eter 

Examined 

Parameter 

Value 

Removed  For 
Analysis 

Param  eter 

Units 

Total  Mass  Flux 
Excluding  Selected 
Point  and 
Interpolating 
(g/day) 

Contribution  of 

Selected  Point  to 

Total  Mass  Flux 

to 
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Figure  3-11.  MTBE  Nearest  Neighbor  Interpolation:  Simple  Uncertainty  Analysis  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Transect  Calculator:  Mass  Flux  Summary 


Site  Location  and  I.D.: 
Description: 


LFR  Site 
MTBE  and  TBA 


Mass  Flux  Summary  for  MTBE 


Distance  from 
Source  (ft) 

Mass  Flux  (g/day) 

Time  Period  1  Time  Period  2  Time  Period  3  Time  Period  4  Time  Period  5 

Transect  1 
Transect  2 

Transect  3 

Transect  4 
Transect  5 

1700 

2.53E+00 

85th  percentile  from  Advanced  Uncertainty  Analysis 
Maximum  of  the  three  Toolkit  interpolation  schemes 

Value  using  selected  interpolation  schemes 

Minimum  of  the  three  Toolkit  interpolation  schemes 
15th  percentile  from  Advanced  Uncertainty  Analysis 


Next  Step: 
Save  File 


Li 


View  Flux  Result  j  |  View  Cone  Grid  jj 

Print 

a  help 

^  Return  To  Toolkit  Main  Screen 

Return  To  Transect  Main  Screen 

Figure  3-12.  MTBE  Nearest  Neighbor  Interpolation:  Mass  Flux  Summary 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Concentration 

Transect  1  Time  Period  1.  Vertical  Interpolation:  Linear  Horizontal  Interpolation:  Linear 


1 

Step  10:  Interpolate  Cone  Horizontally 

J|  Export  Cone  Grid  | 

Import  Cone  Grid 

Back  to  Data  Input 

^  Back  to  Cone  Grid  j| 

Lzl 

helJ 

Next  Step: 

Interpolate  Hydraulic 
Conductivity  Grid 


Start  of  Transect 
0.0 


7.1 

11.6 

16.0 

20.5 
24.9 

29.4 

33.8 

38.2 

42.7 

47.1 

51.6 
56.0 

60.5 

64.9 

69.3 

73.8 

78.2 
82. 7 1 
87.1 

91.6 
96.0 


Data  Input  Instructions 


Enter  value  directly. 
(Values  in 
italics  represent 
interpolated  values.) 

Top/Bottom  of  plume. 
(Dont  enter  any 
data). 


MTBE  Concentrations  (ug/L) 

Distance  from  Edge  of  Transect  (ft) 


MW-1 

100.0 


MW -2 
325.0 


MW -3 
500.0 


MW -4 
665.0 


MW -5 
805.0 


0 

2.50E-01 

5.00E01  Z. 

0 

2.50E-01 

5.00E-01 

1.65E+00 

2.80E+00 

0 

2.50E-01 

5.00E-01 

1.65E+00 

2. 80  E+00 

1. 19E+01 

2.10E+01 

0 

2.50E-01 

5.00E-01 

9.97E+00 

1.94E+01 

2.02E+01 

2.10E+01 

0 

2.50E-01 

5.00E-01 

1.83E+01 

3.61E+01 

2.7  7  E+01 

1.93  E+01 

1.51E+01 

1.10  E+01 

1.15E+01 

1.20  E+01 

0 

2.50E-01 

5.00E-01 

2.66  E+01 

5.27E+01 

3.51  E+01 

1.75E+01 

133E+01 

9.18E+00 

1.06E+01 

1.20  E+01 

0 

2.50E-01 

5.00E-01 

3.49E+01 

6.94E+01 

4.26E+01 

1.58E+01 

1.16E+01 

7.35E+00 

9.55E+00 

1.18E+01 

0 

2  50E-01 

5.00E-01 

4.33  E+01 

8. 60 E+01 

5.00E+01 

1.40  E+01 

9. 76E+00 

5.52E+00 

8.51  E+00 

1.15E+01 

0 

2.50E-01 

5.00E-01 

4.33  E+01 

8. 60  E+01 

5.04E+01 

1.48E+01 

9.23E+00 

3.70  E+00 

7.48E+00 

1.13E+01 

0 

2.50E-01 

5.00E-01 

3.34E+01 

6.63  E+01 

4.09E+01 

1.55E+01 

9.53E+00 

3.56E+00 

7.28E+00 

1.10  E+01 

0 

2.50E-01 

5.00E-01 

2.35E+01 

4.66E+01 

3.14E+01 

1.63  E+01 

9.84E+00 

3.42E+00 

6.36E+00 

9.30E+00 

0 

2.50E-01 

5.00E-01 

1.37E+01 

2.69  E+01 

2.19E+01 

1.70E+01 

1.01  E+01 

3.28E+00 

5.44E+00 

7.60E+00 

0 

2.50E-01 

5.00E-01 

3.85E+00 

7. 20  E+00 

1.21  E+01 

1.70E+01 

1.01  E+01 

3.14E+00 

4. 52  E+00 

5.90E+00 

0 

2.50E-01 

5.00E-01 

3.01  E+00 

5.52E+00 

9.71  E+00 

1.39E+01 

8.45E+00 

3. 00  E+00 

3. 60  E+00 

4. 20  E+00 
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2  50E-01 

5.00E-01 

2.17  E+00 

3. 85  E+00 
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2.17E+00 

2.72E+00 

3.27  E+00 
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_  ■*  sooe-oi 

5.00E-01 

2.55E+00 

4.60  E+00 
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Figure  3-13.  MTBE  Linear  Interpolation:  Concentration  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


in 

O) 

-Q 

£ 

C 


a 


MTBE  Concentrations  (ug/L) 

Distance  from  Edge  of  Transect  (ft) 

MW -6  MW -7  MW -8  End  of  Transect 

905.0  1005.0  1095.0  1185.0  1307.5  1430.0  1565.0  1700.0 


7.1 

11.6 

16.0 

20.5 
24.9 

29.4 
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38.2 
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56.0 
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96.0 
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6.47E+00 


3.63E+00 


5.52E+00 


7.40E+00 


6.85E+00 


6.30E+00 


3.15E+00 


4. 8  3 E+00 


2.07E+00 


3.58E+00 


5.10E+00 


4.73  E+00 


4.37E+00 


2  18E+00 


3.20E+00 


5.00E-01 


1.65E+00 


2. 80  E+00 


2. 62  E+00 


2.43E+00 


1.22E+00 


2.35E+00 


5.00E-01 


5.00E-01 


5.00E-01 


5.00E-01 


5.00E-01 


2.50E-01 


1.89E+00 


5  OOE-OI 


5.00E-01 
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5.00E-01 
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Figure  3-14.  MTBE  Linear  Interpolation:  Concentration  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Hydraulic  Conductivity 

Transect  1  Time  Period  1.  Vertical  Interpolation:  Linear  Horizontal  Interpolation:  Linear 


Step  10:  Interpolate  K  Horizontally 

Export  K  Grid  j| 

Import  K  Grid 

Back  to  Data  Input 

f  Back  to  Cone  Grid 

[  Back  to  K  Grid 

[  Print 

|  HELP 

Next  Step: 

Interpolate  Hydraulic 
Gradient  Grid 


Data  Input  Instructions 


Enter  value  directly. 
(Values  in 
italics  represent 
interpolated  values.) 

Top/Bottom  of  plume . 
(Dont  enter  any 
data). 


Start  of  Transect 
0.0 


7.1 

11.6 

16.0 

20.5 
24.9 

29.4 

33.8 

38.2 

42.7 

47.1 

51.6 
56.0 

60.5 

64.9 

69.3 

73.8 

78.2 

82.7 
87.1 

91.6 
96.0 


Hydraulic  Conductivity  (ft/d) 

Distance  from  Edge  of  Transect  (ft) 


50.0 

MW-1 

100.0 

212.5 

MW -2 

325.0 

412.5 

MW -3 

500.0 

582.5 

MW -4 

665.0 

735.0 

MW -5 

805.0 

1.50E+02 

1.50E+02 

1.50E+02  ■  ■—  •  '  ;  -  — . .  "  —•  -  . 

1 .50  E+02 

1.50E+02 

1.50E+02 

1.50E+02 

1.50  E+02 

1.50  E+02 

1.50E+02 

1.50E+02 

1.50E+02 

1.50  E+02 

9.05E+01 

3.10E+01 

1 .50  E+02 

1.34E+02 

1.18E+02 

1.22E+02 

1.25E+02 

7. 81  E+01 

3.10E+01 

1.50  E+02 

1.18E+02 

8. 60  E+01 

9.32E+01 

1.00  E+02 

7.01E+01 

3.98E+01 

3.04E+01 

2.10E+01 

2. 70  E+01 

3. 30  E+01 

1.50  E+02 

1.02E+02 

5.40E+01 

6.48E+01 

7. 56  E+01 

6.21E+01 

4.85E+01 

3.84E+01 

2.83E+01 

3.06E+01 

3. 30  E+01 

1.50  E+02 

8.60E+01 

2.20E+01 

3.64E+01 

5.08E+01 

5.40E+01 

5.73E+01 

4  64E+01 

3.55E+01 

3.73E+01 

3.90E+01 

1.50  E+02 

8.66E+01 

2.33  E+01 

2.46E+01 

2. 60 E+01 

4.60E+01 

6.60E+01 

5.44E+01 

4.28E+01 

4.39E+01 

4.50E+01 

1.50  E+02 

8.73E+01 

2  45E+01 

2.53E+01 

2.60E+01 

5.65E+01 

8.70E+01 

6.85E+01 

5. 00  E+01 

5.05E+01 

5. 10E+01 

1.50  E+02 

8.79E+01 

2.58E+01 

3.69E+01 

4.80  E+01 

7.80E+01 

1.08E+02 

8.27  E+01 

5.74E+01 

5.72E+01 

5. 70 E+01 

1.50  E+02 

8.85E+01 

2.70E+-01 

4.85E+01 

7.00E+01 

9.95E+01 

1.29E+02 

9.69E+01 

6.48E+01 

7.25E+01 

8.03E+01 

1.50  E+02 

8.42E+01 

1.83  E+01 

5.52E+01 

9.20  E+01 

1.21  E+02 

1.50  E+02 

1.11  E+02 

7.22E+01 

8.78E+01 

1.04E+02 

1.50  E+02 

7.98E+01 

9.67  E+00 

6.18E+01 

1.14E+02 

1.32E+02 

1.50E+02 

1. 15E+02 

7.96E+01 

1.03  E+02 

1.27 E+02 

1.50  E+02 

7.55E+01 

1.00E+00 

4.34E+01 

8.58E+01 

9.93E+01 

1.13E+02 

9.99E+01 

8. 70  E+01 

1.19E+02 

1.50  E+02 

1.50  E+02 

7.55E+01 

1.00  E+00 

2.93E+01 

5.75E+01 

6.65E+01 

7.55E+01 

6.69E+01 

5.83E+01 

8.55E+01 

1.13E+02 

1.50  E+02 

7.55E+01 

1.00E+00 

1.51  E+01 

2.93  E+01 

3.38E+01 

3.83  E+01 

3.40E+01 

2.97E+01 

5.26E+01 

7.55E+01 

1.50  E+02 

7.55E+01 

1.00E+00 

1.00E+00 

1.00  E+00 

1.00E+00 

1.00  E+00 

1.00E+00 
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1.00  E+00 
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. ;  vt+iA' 
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1.00E+00 
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Figure  3-15.  Linear  Interpolation:  Hydraulic  Conductivity  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


in 

U) 

JQ 

& 

C 


O 


Hydraulic  Conductivity  (ft/d) 

Distance  from  Edge  of  Transect  (ft) 

MW -6  MW -7  MW -8  End  of  Transect 

905.0  1005.0  1095.0  1185.0  1307.5  1430.0  1565.0  1700.0 


7.1 

11.6 

16.0 

20.5 
24.9 

29.4 

33.8 

38.2 

42.7 

47.1 

51.6 
56.0 

60.5 

64.9 

69.3 

73.8 

78.2 

82.7 
87.1 

91.6 
96.0 


8.22E+01 


1.19E+02 


8.17E+01 


4.40E+01 


4. 40E+01 


4.40E+01 


2.25E+01 


1.00E+00 


6.98E+01 


8.87E+01 


6.63E+01 


4.40E+01 


4.52E+01 


4.63E+01 


2.37E+01 


1.00E+00 


5.75E+01 


5. 80  E+01 


5. 28E+01 


4.75E+01 


4.81  E+01 


4.87  E+01 


2.48E+01 


1.00E+00 


5.96E+01 


3.90E+01 


4.50E+01 


5.10E+01 


5.10E+01 


5.10E+01 


2. 60  E+01 


1.00E+00 


6.18E+01 


2.00E+01 


2.72E+01 


3. 43  E+01 


3. 43  E+01 


3.43E+01 


1.77  E+01 


1.00E+00 


6.39E+01 


I.OOE+OO 


9.33E+00 


1.77  E+01 


1.77  E+01 


1.77  E+01 


9.33E+00 


1.00E+00 


7.55E+01 


1.00E+00 


1.00E+00 


1.00E+00 


1.00E+00 


1.00E+00 


1.00E+00 


1.00E+00 


5 69E+01 


1.00E+00 


1.00E+00 


1.00E+00 


1.00E+00 


1.00E+00 


1.00E+00 


1.00E+00 


3.83E+01 


I.OOE+OO 


1.00E+00 


I.OOE+OO 


I.OOE+OO 


I.OOE+OO 


I.OOE+OO 


I.OOE+OO 


Figure  3-16.  Linear  Interpolation:  Hydraulic  Conductivity  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Hydraulic  Gradient 

Transect  1  Time  Period  1.  Vertical  Interpolation:  Linear  Horizontal  Interpolation:  Linear 


Step  10:  Interpolate 

Horizontally 

|  Export  i  Grid 

Import  i  Grid 

Next  Step: 

Back  to  Data  Input  J 

Back  to  Cone  Grid  j| 

Back  to  K  Grid 

j|  Print  |  HELP 

Calculate  Mass  Flux 

Data  Input  Instructions 

Enter  value  directly. 

110.80  1 

^  (Values  in 

italics  represent 

interpolated  values.) 

Top/Bottom  of  plume 

HI 

^  (Don't  enter  any 

data). 

Start  of  Transect 
0.0 


7.1 

11.6 

16.0 

20.5 
24.9 

29.4 

33.8 

38.2 

42.7 

47.1 

51.6 
56.0 

60.5 

64.9 

69.3 

73.8 

78.2 

82.7 
87.1 

91.6 
96.0 
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CASE  STUDY:  MTBE  AND  TBA 
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CASE  STUDY:  MTBE  AND  TBA 
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3.  Black  shaded  cells  represent  the  top  and  bottom  of  the  plume. 
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CASE  STUDY:  MTBE  AND  TBA 
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CASE  STUDY:  MTBE  AND  TBA 
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End  of  Transect 


End  of  Transect 


Top  of 
Sampling 
Interval 


40 


Bottom  of 
Sampling 
Interval 


45 


67 


Parameter 

Examined 


Parameter 

Value 

Removed  For 
Analysis 


2.13E-03 


2.13E-03 


Parameter 

Units 


RESULTS 


Total  Mass  Flux 
Excluding  Selected 
Point  and 
Interpolating 
(g/day) 

Contribution  of 

Selected  Point  to 

Total  Mass  Flux 

(%) 

2.96  E+00 

-3.8 

2.96  E+00 

•3.8 

2.96  E+00 

-3.8 

2.95  E+00 

•3.3 

2.96  E+00 

-3.8 

2.96  E+00 

-3.8 

2.96  E+00 

-3.8 

2.95  E+00 

-3.5 

2.96  E+00 

-3.8 

2.96  E+00 

-3.8 

2.96  E+00 

-3.8 

2.95  E+00 

-3.5 

2.96  E+00 

-3.8 

2.96  E+00 

-3.8 

2.96  E+00 

-3.8 

2.92  E+00 

-2.3 

2.96  E+00 

•3.8 

2.96  E+00 

-3.8 

2.95  E+00 

-3.2 

2.96  E+00 

-3.8 

2.96  E+00 

-3.8 

2.96  E+00 

-3.6 

2.96  E+00 

-3.8 

2.96  E+00 

•3.8 

2.96  E+00 

-3.8 

Figure  3-22.  MTBE  Linear  Interpolation:  Simple  Uncertainty  Analysis  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Transect  Calculator:  Mass  Flux  Summary 


Site  Location  and  I.D.: 
Description: 


LFR  Site 
MTBE  and  TBA 


Mass  Flux  Summary  for  MTBE 


Distance  from 

Mass  Flux  (g/day) 

Source  (ft) 

Time  Period  1 

Time  Period  2  Time  Period  3  Time  Period  4 

Time  Period  5 

Transect  1 
Transect  2 

Transect  3 

Transect  4 
Transect  5 

1700 

2.85E+00 

85th  percentile  from  Advanced  Uncertainty  Analysis 
Maximum  of  the  three  Toolkit  interpolation  schemes 

Value  using  selected  interpolation  schemes 
Minimum  of  the  three  Toolkit  Interpolation  schemes 
15th  percentile  from  Advanced  Uncertainty  Analysis 


Next  Step: 

View  Flux  Result  X  View  Cone  Grid  X  Print  If  HELP  J 

Save  File 

Return  To  Toolkit  Main  Screen  ‘if  Return  To  Transect  Main  Screen 

Figure  3-23.  MTBE  Linear  Interpolation:  Mass  Flux  Summary 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Concentration 

Transect  1  Time  Period  1.  Vertical  Interpolation:  Log  Horizontal  Interpolation:  Log 


Step  10:  Interpolate  Cone  Horizontally 

Export  Cone  Grid 

Import  Cone  Grid 

_ 

Back  to  Data  Input  J 

^  Back  to  Cone  Grid  | 

Back  to  K  Grid 

[  Print  | 

(  HELP 

Next  Step: 
Interpolate  Hydraulic 
Conductivity  Grid 


□ 


Data  Input  Instructions 


Enter  value  directly 
(Values  in 
italics  represent 
interpolated  values.) 

Top/Bottom  of  plume. 
(Dont  enter  any 
data). 


MTBE  Concentrations  (ug/L) 

Distance  from  Edge  of  Transect  (ft) 


Start  of  Transect 

MW  -1 

MW -2 

MW -3 

MW -4 

MW -6 

0.0 

50.0 

100.0 

212.5 

325.0 

412.5 

500.0 

582.5 

665.0 

735.0 

805.0 

7.1 

0 

7.07E-04 

5.00E-01  ^  .  ■;*  - ( 

11.6 

0 

7.07E-04 

5.00E-01 

1.18E+00 

2.80E+00  -  -  '  ■  -  ; 

16.0 

0 

7.07E-04 

5.00E-01 

1.18E+00 

2. 80  E+00 

7.67E+00 

2.10E+01 

20.5 

0 

7.07E-04 

5.00E-01 

1.67E+00 

5.55E+00 

1.08E+01 

2.10E+01 

24.9 

0 

7.07E-04 

5.00E-01 

2.35E+00 

1.10E+01 

1.45E+01 

1.90E+01 

1.44E+01 

1.10E+01 

1.15E+01 

1.20E+01 

29.4 

0 

7.07E-04 

5.00E-01 

3.31  E+00 

2.19E+01 

1.94E+01 

1.71E+01 

120E+01 

8.38E+00 

1.00E+01 

1.20E+O1 

33.8 

0 

7.07E-04 

5.00E-01 

4.66E+00 

4.34E+01 

2.59E+01 

1.55E+01 

9.94E+00 

6.38E+00 

8.65E+00 

1.17E+01 

38.2 

0 

7.07E-04 

5.00E-01 

6.56  E+00 

8.60E+01 

3.47E+01 

1.40E+01 

8.25E+00 

4.86E+00 

7. 47  E+00 

1.15E+01 

42.7 

0 

7.07E-04 

5.00E-01 

6.56  E+00 

8.60E+01 

3.56E+01 

1.47E+01 

7. 37  E+00 

3.70  E+00 

6.45E+00 

1.12E+01 

47.1 

0 

7.07E-04 

5.00E-01 

4.81  E+00 

4.63E+01 

2.67E+01 

1.54E+01 

7.40E+00 

3.55E+00 

6.25E+00 

1.10E+01 

51.6 

0 

7.07E-04 

5.00E-01 

3. 5  3  E+00 

2.49E+01 

2.01E+01 

1.62E+01 

7.42E+00 

3.40E+00 

5.42E+00 

8.65E+00 

56.0 

0 

7.07E-04 

5.00E-01 

2. 59  E+00 

1.34E+01 

1.51E+01 

1.70E+01 

7.45E+00 

3.26E+00 

4.71  E+00 

6.80E+00 

60.5 

0 

7.07E-04 

5.00E-01 

1.90E+00 

7. 20  E+00 

1.11E+01 

1.70E+01 

7.29E+00 

3.13E+00 

4.09E+00 

5.34E+00 

64.9 

0 

7.07E-04 

5.00E-01 

1.36E+00 

3.70E+00 

6.73E+00 

1.23E+01 

6.06E+00 

3. 00  E+00 

3.55E+00 

4.20  E+00 

69.3 

0 

7.07E-04 

5.00E-01 

9.74E-01 

1.90  E+00 

4.10E+00 

8. 84  E+00 

3.82E+00 

1.65E+00 

2.02E+00 

2.47  E+00 

73.8 

0 

7.07E-04 

5.00E-01 

6.98E-01 

9.74E-01 

2.49E+00 

6.38E+00 

2.41  E+00 

9.09E-01 

1.15E+00 

1.45E+00 

78.2 

0 

7.07E-04 

5.00E-01 

5.00E-01 

5.00E-01 

1.52E+00 

4. 60  E+00 

1.52E+00 

5.00E-01 

6.52E-01 

8.51  E-01 

82.7 

5  0QE-01 

5.00E-01 

1.52E+00 

4.60  E+00 

1.52E+00 

5.00E-01 

5.00E-01 

5.00E-01 

87.1 

: .  ....  5  00E-0 > 

5.00E-01 

1.52E+00 

4.60  E+00 

1.52E+00 

5.00E-01 

5.00E-01 

5.00E-01 

91.6 

5.00E-01 

5.00E-01 

5.00E-01 

5.00E-01 

Figure  3-24.  MTBE  Log  Transformation:  Concentration  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


MTBE  Concentrations  (ug/L) 

Di stan ce  from  Edge  of  Transect  (ft) 

MW -6  MW -7  MW -8  End  of  Transect 

905.0  1005.0  1095.0  1185.0  1307.5  1430.0  1565.0  1700.0 


7.1 

11.6 

16.0 

20.5 

24.9 

S,  29.4 

£  33.8 


S'  42.7 

a 

47.1 

51.6 
56.0 

60.5 
64.9 
69.3 
73.8 

78.2 

82.7 
87.1 

91.6 
96.0 


8.53E+00 

6.20  E+00 

4.20E+00 

4.20E+00 

2.05E-03 

0 

8.20E+00 

5.85E+00 

6.35E+00 

6.90E+00 

5.38E+00 

4. 20  E+00 

2. 05E-03 

0 

7.87E+00 

5.5 1E+00 

6. 17E+00 

6.90E+00 

5.76E+00 

4.81E+00 

2.19E-03 

0 

7  56E+00 

5.20E+C0 

6. 10E+00 

7.15E+00 

6.27E+00 

5. 50  E+00 

2.35E4)3 

0 

4.54E+00 

2.38E+00 

4.20E+00 

7.40E+00 

6.83E+00 

6. 30  E+00 

2.51  E-03 

0 

2.72E+00 

1.09E+00 

1.81E+00 

3.01  E+00 

2.86  E+00 

2.71  E+00 

1.65E-03 

0 

1.63E+00 

5.00E-01 

7.83E-01 

1.23E+00 

1.20E+00 

1 16  E+00 

108E-03 

0 

1.45E+00 

5.00E-01 

5.00E-01 

5.00E-01 

5.00E-01 

5.00E-01 

7.07E-04 

0 

1.11E+00 

5.00E-01 

5.00E-01 

5.00E-01 

5.00E-01 

5.00E-01 

7.07E-04 

0 

8.51E-01 

5.00E-01 

5.00E-01 

5 00E-01 

5.00E-01 

5.00E401 

7.07E434 

0 

5.00E-01 

5.00E-01 

5.00E-01 

5.00E-01 

7.07E-04 

0 

Figure  3-25.  MTBE  Log  Transformation:  Concentration  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Hydraulic  Conductivity 

Transect  1  Time  Period  1.  Vertical  Interpolation:  Log  Horizontal  Interpolation:  Log 


Step  10:  Interpolate  K  Horizontally 

Export  K  Grid 

Import  K  Grid 

1 1  Back  to  Data  Input 

|  Back  to  Cone  Grid  j 

| 

(  Print  1 

[  HELP 

Next  Step: 

Interpolate  Hydraulic 
Gradient  Grid 


Start  of  Transect 
0.0 


Hydraulic  Conductivity  (ft/d) 

Distance  from  Edge  of  Transect  (ft) 


MW  -1 
100.0 


MW -2 
325.0 


MW-3 

500.0 


Data  Input  Instructions 

Enter  value  directly. 

110.80  1  _ 

(Values  in 

1  .  1  V 

italics  represent 

interpolated  values.) 

Top/Bottom  of  plume. 

— ► 

(Don't  enter  any 

data). 

MW -4 

665.0 

735.0 

MW -5 
805.0 


7.1 

1.50  E+02 

1.50E+02 

1.50E  +  02  ..v:-..-  I,  :  . „ 

11.6 

1.50  E+02 

1.50E+02 

1.50E+02 

1.50E+02 

1.50E+02 

16.0 

1.50  E+02 

1.50E+02 

1.50E+02 

1.50E+02 

1.50  E+02 

6.82E+01 

3.10E+01 

20.5 

1.50  E+02 

1. 1 8E+02 

9.28E+01 

9. 90  E+01 

1.06  E+02 

5.72E+01 

3.10E+01 

24.9 

1.50  E+02 

9.28E+01 

5.74E+01 

6.54E+01 

7. 44  E+01 

5.28E+01 

3.74E+01 

2.80E+01 

2.10E+01 

2. 63  E+01 

3. 30  E+01 

29.4 

1.50  E+02 

7.30E+01 

3.56E+01 

4.32E+01 

5.24E+01 

4. 87  E+01 

4.52E+01 

3.44E+01 

2.61E+01 

2.93E+01 

3. 30  E+01 

33.8 

1.50  E+02 

5.74E+01 

2.20E+01 

2.85E+01 

3.69E+01 

4.49E+01 

5.46  E+01 

4.21  E+01 

3.24E+01 

3.50E+01 

3.78E+01 

38.2 

1.50  E+02 

5.89E+01 

2.32E+01 

2.45E+01 

2. 60  E+01 

4. 14E+01 

6. 60  E+01 

5. 15E+01 

4.03E+01 

4.18E+01 

4.34E+01 

42.7 

1.50  E+02 

6.05E+01 

2.44E+01 

2.52E+01 

2. 60  E+01 

4.59E+01 

8.10E+01 

6.37E+01 

5. 00  E+01 

4. 99  E+01 

4.97E+01 

47.1 

1.50  E+02 

6.20E+01 

2.57  E+01 

3.11  E+01 

3.76E+01 

6.12E+01 

9.95E+01 

7.46E+01 

5.59E+01 

5.64E+01 

5. 70  E+01 

51.6 

1.50  E+02 

6.36E+01 

2.70E+01 

3.83E+01 

5.44E+01 

8.16E+01 

1.22  E+02 

8.73E+01 

6.24E+01 

6.73E+01 

7.26E+01 

56.0 

1.50  E+02 

3.67E+01 

9.00E+00 

2.66E+01 

7.88E+01 

1.09E+02 

1.50  E+02 

1.02E+02 

6.97E+01 

8.03  E+01 

9.25E+01 

60.5 

1.50  E+02 

2.12E+01 

3.00E+00 

1.85E+01 

1.14E+02 

1.31  E+02 

1.50  E+02 

1  08E+02 

7.79E+01 

9.58E+01 

1.18E+02 

64.9 

1.50  E+02 

1.22E+01 

1.00E+00 

5.91  E+00 

3.49E+01 

3. 87  E+01 

4.29E+01 

6.11  E+01 

8. 70  E+01 

1.14E+02 

1.50  E+02 

69.3 

1.50  E+02 

1.22E+01 

1.00E+00 

3.27E+00 

1.07  E+01 

1.14E+01 

1.22  E+01 

1.55E+01 

1.96E+01 

2.90E+01 

4.29E+01 

73.8 

1.50  E+02 

1.22E+01 

1.00  E+00 

1.81  E+00 

3.27  E+00 

3.38E+00 

3. 50  E+00 

3.94E+00 

4. 4  3 E+00 

7.37E+00 

1.22E+01 

78.2 

1.50  E+02 

1.22E+01 

1.00E+00 

0.00E+00 

1.00  E+00 

0.00E+00 

1.00  E+00 

0. 00 E+00 

1.00  E+00 

1.87  E+00 

3. 50  E+00 

82.7 

'  j  00E+G0 

1.00E+00 

0.00E+00 

1.00  E+00 

0.00E+00 

1.00  E+00 

0.00  E+00 

1.00  E+00 

87.1 

.  1  0QE+00 

1.00E+00 

0.00E+00 

1.00  E+00 

0.00E+00 

1.00E+00 

0.00E+00 

1.00E+00 

91.6 

100E+00 

1.00E+00 

0.00  E+00 

1.00E+00 

Figure  3-26.  Log  Transformation:  Hydraulic  Conductivity  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


O) 

-C 


£ 


Q 


Hydraulic  Conductivity  (ft/d) 

Distance  from  Edge  of  Transect  (ft) 

MW -6  MW -7  MW-8  End  of  Transect 

905.0  1005.0  1095.0  1185.0  1307. 5  1430.0  1565.0  1700.0 


7.1 

11.6 

16.0 

20.5 
24.9 

29.4 

33.8 

38.2 

42.7 

47.1 

51.6 
56.0 

60.5 

64.9 

69.3 

73.8 

78.2 

82.7 
87.1 

91.6 
96.0 


7.5  3E+01 


6.88E+01 


1.50E+02 


4.40E+01 


1.09E+02 


6.93E+01 


4.40E+01 


4.40E+01 


4.40E+01 


4.40E+01 


6.63E+00 


6.63E+00 


1.00E+00 


1.00E+00 


6.29E+01 


7  96E+01 


5.92E+01 


4.40E+01 


4.51E+01 


4.62E+01 


6.80E+00 


1.00E+00 


5.75E+01 


5.80E+01 


5.24E+01 


4.74E+01 


4.80E+01 


4.86E+01 


6.97E+00 


1.00E+00 


3.30E+01 


1  50E+01 


2  76E+01 


5.10E+01 


5.10E+01 


5.10E+01 


7  14E+00 


1.00E+00 


1.89E+01 


387E+00 


7.30E+00 


1.38E+01 


1.38E+01 


1.38E+01 


3.71E+00 


1.00E+00 


1.09E+01 


I.OOE+OO 


1.93E+00 


3.71  E+00 


3.71  E+00 


3.71  E+00 


1.93E+00 


1.00E+00 


1.22E+01 


100E+00 


0.00E+00 


1.00E+00 


O.OOE+OO 


1.00E+00 


0.00E+00 


1.00E+00 


6.55E+00 


1.00E+00 


O.OOE+OO 


1.00E+00 


O.OOE+OO 


1.00E+00 


O.OOE+OO 


1.00E+00 


3.50E+00 


I.OOE+OO 


O.OOE+OO 


I.OOE+OO 


O.OOE+OO 


I.OOE+OO 


O.OOE+OO 


I.OOE+OO 


O.OOE+OO 


I.OOE+OO 


O.OOE+OO 


I.OOE+OO 


Figure  3-27.  Log  Transformation:  Hydraulic  Conductivity  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Hydraulic  Gradient 

Transect  1  Time  Period  1.  Vertical  Interpolation:  Log  Horizontal  Interpolation:  Log 

3 


Step  10:  Interpolate  i  Horizontally 


Export  i  Grid 


Import  i  Grid 


Back  to  Data  Input 


- 


Back  to  Cone  Grid 


Back  to  K  Grid 


I 


HELP 


Next  Step: 
Calculate  Mass  Flux 


□ 


Data  Input  Instructions 


10.80 


Enter  value  directly. 
(Values  in 
italics  represent 
interpolated  values.) 

Top/Bottom  of  plume 
(Dont  enter  any 
data). 


Hydraulic  Gradient  (ft/ft) 

Distance  from  Edge  of  Transect  (ft) 


Start  of  Transect 
0.0 

50.0 

MW-1 

100.0 

212.5 

MW -2 

325.0 

412.5 

MW-3 

500.0 

582.5 

MW -4 

665.0 

735.0 

MW -5 

805.0 

7.1 

7.10E-04 

7.10E-04 

7.10E-04  L*  .  . *  —  ,*  ..  .  :  •  ...l 

11.6 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

16.0 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

2. 13E-03 

2.13E-03 

20.5 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

2.13E-03 

2.13E-03 

24.9 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

29.4 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

33.8 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

38.2 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

2. 13E-03 

2.13E-03 

2. 13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

42.7 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

2 13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

47.1 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

51.6 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

56.0 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2. 13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

60.5 

7.10E-04 

7.10E-04 
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2.13E-03 
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2.13E-03 

2.13E-03 

2.13E-03 

69.3 

7.10E-04 

7.10E-04 

7.10E-04 

1.23E-03 

2.13E-03 
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2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 

2.13E-03 
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2.13E-03 
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2. 13E-03 
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2.13E-03 
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Figure  3-28.  Log  Transformation:  Hydraulic  Gradient  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


u> 

U3 

SI 


Hydraulic  Gradient  (ft/ft) 

Distance  from  Edge  of  Transect  (ft) 

MW -6  MW-7  MW -8  End  of  Transect 

905.0  1005.0  1095.0  1185.0  1307.5  1430.0  1565.0  1700.0 


7.1 

11.6 

16.0 

20.5 
24.9 

29.4 

33.8 

38.2 

42.7 

47.1 

51.6 
56.0 

60.5 

64.9 

69.3 

73.8 

78.2 

82.7 
87.1 

91.6 
96.0 


2.02E-03 


2.02E-03 


1.90E-03 


1.90E-03 


1.91E-03 


1.90E-03 


1.9  IE-03 


1.90E-03 


1.90E-03 


1.91E-03 


1.91E-03 


2.02E-03 


1.90E-03 


1.90E-03 


1.91E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.91E-03 


2.02E-03 


1.91E-03 


1. 90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.91E-03 


2.02E-03 


1.90E-03 


1.90E-03 


1.91E-03 


1.90E-03 


1.91E-03 


1. 90E-03 


1.91E-03 


2.02E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.91E-03 


2.02E-03 


1.91E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.91E-03 


2.02E-03 


1.90E-03 


1.90E-03 


1.91E-03 


1  90E-03 


1.9  IE-03 


1.90E-03 


1.91E-03 


2.02E-03 


1.90E-03 


1.90E-03 


1.91E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.91E-03 


2.02E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.91E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.90E-03 


1.91E-03 


Figure  3-29.  Log  Transformation:  Hydraulic  Gradient  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Mass  Flux  Result 


TOTAL  MASS  FLUX 


|(g/day) 


|  (kg/yr) 


Data  Representation 

1.  Bold  values  represent  calculations  based  on  given  values. 

2.  Values  in  italics  represent  calculations  based  on  interpolation. 

3.  Black  shaded  cells  represent  the  top  and  bottom  of  the  plume. 


|Next  Step:  Mass  Flux 

Run/View  Uncertainty  1 

Via*/  Final  Concentration  Grid 

Summary 

Analysis  (Optional) 

Back  to  Data  Grid  j|  Print  HELP  ]  1 

Start  of  Transect 
0.0 


7.1 

11.6 

16.0 

20.5 
24.9 
294 

33.8 

38.2 

42.7 

47.1 

51.6 
56.0 

60.5 

64.9 

69.3 

73.8 

78.2 

82.7 
87.1 

91.6 
96.0 


MTBE  Mass  Flux  (g/day) 

Distance  from  Edge  of  Transect  (ft) 
MW-1 

100.0  212.5 


rvrw-2 

325.0 


SELECT  TRANSECT  TO  VIEW 
SELECT  TIME  PERIOD  TO  VIEW 


MW -3 
500.0 


MAM 

665.0 


5.45E-04 
5.45E-04  3.09E-03 


MW -6 
805.0 


Figure  3-30.  MTBE  Log  Transformation:  Mass  Flux  Results  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


■9 

e 


3. 

a> 

o 


MTBE  Mass  Flux  (g/day) 

Distance  from  Edge  of  Transect  (ft) 


42.7 


47.1 


51.6 

56.0 


60.5 


64.9 


905.0 


1.43E-02 


1.26E-02 


1.10E-02 


3.80E-03 


1.31  E-03 


4.50E-04 


4.50E-04 


1.84E-04 


7.56E-05 


rvwv-6 

1005.0 


1.46E-02 


1.00E-02 


6.87E-03 


8. 13E-04 


9.62E-05 


1. 14E-05 


1. 14E-05 


1. 14E-05 


9.50E-03 


7.88E-03 


6.90E-03 


2.50E-03 


2.86E-04 


3.26E-05 


0.00E+00 


0. 00E+00 


O.OOE+OO 


7.73E-03 


7.73E-03 


8.62E-03 


1. 06E-03 


1. 16E-04 


1.27E-05 


1.27E-05 


6.96E-03 


7.63E-03 


8.83E-03 


1.02E-02 


1. 15E-03 


1.30E-04 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


5.71  E-03 


5.71  E-03 


6.86E-03 


8.25E-03 


1. 15E-03 


1.33E-04 


1.54E-05 


1.54E-05 


1.54E-05 


1565.0 


4.40E-07 


4.40E-07 


4.83E-07 


5. 29E-07 


5.80E-07 


1.98E-07 


6.72E-08 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


1700.0 


End  of  Transect 


O.OOE+OO 


0. 00E +00 


O.OOE+OO 


O.OOE+OO 


O.OOE+OO 


O.OOE+OO 


O.OOE+OO 


Figure  3-31.  MTBE  Log  Transformation:  Mass  Flux  Results  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Interpolation  Error  Results 


Site  Location  and  I  D.: 
Description: 


LFR  Site 
MTBE  and  TBA 


unsmutntM 

1.  the  table  below  shorn  all  oorcene  anon  and,  if  apciicabie, 
ncrHjnfarm  ccnductnnty/Darcy  veloatv  and  gadent 
incut  data. 

2.  During  the  unoertaintr  analysis.  ea«Ti  of  tie  cents  in  the 
table  sre  remewed  ere  at  a  time  and  the  mass  flux 
catenated  Mthout  that  pont  using  the  inter  potation 


Next  Step:  1 

Back  to  Mass  Flux  Result  1 

Print 

SELECT  TRANSECT  TO  VIEW 

transect  I  w 

Mass  Flux  Summary  [j 

Back  to  Data  Input 

HELP 

SELECT  TIME  PERIOD  TO  VIEW 

>  * 

MTBE  Interpolation  Methods 

Hydraulic  ConiAjciivty  1)  Vertical  Log  Transformation  2)  Horizontal,  Lug  Transformation 
Concentration.  1)  Vertical  Log  Transformation  2)  Horizontal  Log  Transformation 
I  lydrmihr  Gradient  1)  Vertical  Log  Transformation  2)  Hortrontal  l  og  Transformation 


Total  Mass  Flux  Including  All  Points 


I  (9<toy) 


End  of  Transect 

Top  of 
Sampling 
Interval 

Bottom  of 
Sampling 
Interval 

Parameter 

Examined 

Parameter 

Value 

Removed  For 
Analysis 

Parameter 

Units 

Total  Ma«»  Flux 
Excluding  Selected 
Point  and 
Interpolating 

Contribution  of 
Selected  Point  to 

Total  Mass  Flux 

IS 

Start  of  Transect 

Conductivity 

1  50E-02 

ft/d 

1.81  E+00 

<4 

MW- 1 

50 

20  0 

Conductivity 

1  50E  *02 

ft/d 

1.80000 

-5.9 

MW  1 

330 

380 

Conductivity 

2.20E+01 

ft/d 

182000 

-6.8 

MW-1 

51  0 

56  0 

Conductivity 

2  70E  *01 

ft/d 

1.80000 

■6.0 

MW-1 

66  0 

710 

Conductivity 

1  OOC'OO 

ft/d 

1.81  E+00 

■6.6 

MW  2 

7.6 

22  0 

Conductivity 

1  COE *02 

ft/d 

1.68E+00 

1.0 

MW-2 

40  0 

45  0 

Conductivity 

2  60001 

ft/d 

2.85  000 

-67.6 

MW-2 

62  0 

67  0 

Conductivity 

1  14C  *02 

ft/d 

1.61000 

5.3 

MW  2 

77.0 

82  0 

Conductivity 

1  OOE+OO 

ft/d 

1.84000 

•8.4 

MW-3 

10.0 

25  0 

Conductivity 

3 1 0C  *01 

ft/d 

1.88000 

-10.7 

MW-3 

30  0 

430 

Conductivity 

6  60001 

ft/d 

1.81000 

-6.7 

MW  3 

G8  0 

030 

Conductivity 

1  50E+02 

ft/d 

144E+00 

16.2 

MW-3 

77.0 

02  0 

Conductivity 

1 00C  *00 

ft/d 

1.96E+00 

•15.1 

MW-4 

15.0 

300 

Conductivity 

2.10001 

ft/d 

1.83000 

-7.7 

MW  4 

44  0 

400 

Conductivity 

5.00E+01 

ft/d 

179000 

-6  6 

MW-4 

63.0 

68  0 

Conductivity 

8  706 >01 

ft/d 

1.73000 

•1.7 

MW-4 

80  0 

85  0 

Conductivity 

>00E-00 

ll/d 

1.82E+00 

•T2 

Mwr, 

200 

35  0 

Conductivity 

3 306*01 

11/d 

1.84E+00 

-8.0 

MW-5 

47  0 

52  0 

Conductivity 

5  706>01 

ft/d 

1.83  E+00 

-7.5 

MW-5 

65  0 

TOO 

Conductivity 

1 10E-02 

ft/d 

1.74000 

•2.4 

MW  f> 

010 

06  0 

Conductivity 

1006*00 

ft/d 

183000 

•7  6 

MW-6 

24  0 

39  0 

Conductivity 

1  50002 

ft/d 

1.73000 

-1.6 

MW -6 

47.0 

5 2  0 

Conductivity 

5  S0E-0I 

ft/d 

1.76000 

•3.7 

Mwr, 

620 

67  0 

Conductivity 

1  006*00 

ft/d 

183  E+00 

-7.8 

MW-7 

30.0 

45  0 

Conductivity 

4  406*01 

ft/d 

1.80000 

-6.6 

MW -7 

5  2  0 

5  TO 

1 Conductivity 

5 10E*01 

ft/d 

1.78000 

-1.6 

MW  7 

67  0 

72  0 

Conductivity 

1  006*00 

ft i  d 

1.82000 

-6.9 

MW-3 

29 

44 

Conductivity 

4  406*01 

ft/d 

1.81000 

-6.3 

MW-6 

51 

56 

Conductivity 

5  WE-01 

ft/d 

1.78000 

■4.4 

MWfi 

66 

71 

Conductivity 

1  006*00 

ff/d 

1.82  E+00 

43 

6nd  of  Transect 

Conductivity 

1006*00 

ft/d 

1.81000 

-6.4 

Start  of  Transect 

iConcvntration 

0 00E*0Q 

uoi 

1.81000 

■6.7 

MW  1 

5 

20 

Concentration 

5  0064)1 

ug/L 

1,81000 

«4 

MW-1 

33 

38 

Concentration 

5  0000/ 

ug/L 

1.81000 

4.4 

MW-1 

51 

56 

Concentration 

5  00E4I 

ug/L 

1.81000 

■6.4 

MW  r 

66 

71 

Concentration 

5  00601 

1.81E+00 

-6  4 

MVS-2 

7 

22 

Concentration 

2  80000 

ug/L 

3.19E+00 

■67.4 

MVS-2 

40 

45 

Concentration 

8  60007 

ug/L 

1.36000 

20.2 

MVS  2 

6? 

67 

Concentration 

--00 

tig/) 

1.80000 

■6.7 

MW-2 

17 

82 

Concentration 

5  006-01 

ug/L 

1.82000 

-7.0 

MW-3 

id 

25 

Concentration 

2  10E*01 

ug/L 

1.76E+0Q 

•3.6 

MW .? 

3a 

43 

Cnncentratim 

1  406*01 

ug/t 

1.87  E+00 

-9.8 

MVS-3 

58 

63 

Concentration 

1  7000/ 

ug/L 

1.65000 

3.2 

MVS-3 

7? 

82 

Concentration 

4 60E*00 

ug/L 

1.82000 

-7.1 

MW  4 

15 

30 

Cnncentraton 

1  106*01 

ua/L 

1.77000 

-4.3 

MVS-4 

44 

49 

Concentration 

3  706*00 

ug/L 

1.84000 

•6.4 

MW-4 

63 

68 

Concentration 

3  006*00 

uy.'L 

1.77000 

-4.1 

MW  4 

a: 

or, 

Concentration 

5  006  01 

1.81E+00 

-6  5 

MVS-5 

20 

35 

Concentration 

7  2000/ 

ug/L 

1.79000 

■6.6 

MW-5 

47 

52 

^Concentration 

1  106*01 

ug/L 

1.78000 

-4.4 

MW  5 

66 

70 

Concentration 

4  206*00 

ug/) 

1.78E+00 

-6.0 

MW-5 

3: 

86 

Concentration 

5  006-01 

ug-L 

1.81000 

■6.6 

MW-6 

24 

39 

Concentration 

6  206*00 

ug-l 

1.76E+00 

•3.3 

Mwr, 

47 

52 

Concentration 

5  206*00 

ug/) 

1.77E+00 

-4.3 

MVS-6 

62 

67 

Concentration 

5.0000/ 

ugft- 

1.82E+00 

-6.9 

MW-7 

30 

45 

Concentration 

0 906*00 

ug/L 

1.80  E+00 

-5.6 

MW7 

52 

57 

Concentration 

7 406*00 

tig/) 

1,79000 

MVS-7 

67 

72 

Concentration 

5  006-01 

ug/L 

1.81000 

-6.6 

MW-6 

53 

44 

Concentroton 

4  206*00 

ug/L 

1.82  E+00 

43 

MW-3 

51 

56 

Concentration 

6  306*00 

ug* 

1,78000 

-4  8 

MVS-6 

66 

71 

Concentration 

5  006-01 

ug/L 

1.81E+00 

-6.6 

End  of  Transect 

Concentration 

0  006*00 

ug/L 

1.61  E+00 

■6.7 

Start  of  Transect 

Gradient 

7  106  04 

ftm 

161E+0C 

4.4 

MW-1 

5 

20 

Gradient 

7  106-04 

ft/ft 

1.81000 

4.4 

MW-1 

33 

38 

Gradient 

7  106-04 

ft/n 

1.81  E+00 

4.4 

MW  1 

51 

56 

Gradient 

7  106  04 

m 

1.81E+00 

44 

MVS-1 

66 

7 1 

Gradient 

7  106-04 

m 

1.81000 

4.4 

MW-2 

1 

22 

Gradient 

2  t3E43 

m 

1.80  E+00 

_ -91 _ 

Figure  3-32.  MTBE  Log  Transformation:  Simple  Uncertainty  Analysis  -  Page  1 
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69  ' 

70  ' 

71  ' 

72  ' 

73  ' 

74  ' 

75 

75  ' 

77  " 

78 

79  ' 

80  ' 

87 

82  ' 

83  ' 

84 

85  ' 

88  ' 

87 

88  ' 

89  ' 

90 

97  ' 

92  ' 
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CASE  STUDY:  MTBE  AND  TBA 


RESULTS 


End  of  Transect 


End  of  Transect 


Top  of 
Sampling 
Interval 


40 


Bottom  of 
Sampling 
Interval 


45 


Parameter 

Examined 


Parameter 

Value 

Removed  For 
Analysis 


2.13E-03 


2.13E-03 


Parameter 

Units 


Total  Mass  Flux 
Excluding  Selected 
Point  and 
Interpolating 
(g/day) 

Contribution  of 

Selected  Point  to 

Total  Mass  Flux 
(%! 

1.81E+00 

■6.4 

1.81  E+00 

-6.4 

1.81  E+00 

-6.4 

1.79  E+00 

■6.6 

1.81  E+00 

-6.4 

1.81  E+00 

-6.4 

1.81  E+00 

■6.4 

1.80  E+00 

-6.8 

1.81  E+00 

-6.4 

1.81  E+00 

-6.4 

1.81  E+00 

-6.4 

1.80  E+00 

-6.8 

1.81  E+00 

-6.4 

1.81  E+00 

-6.4 

1.81  E+00 

•6.4 

1.77  E+00 

-3.9 

1.81  E+00 

-6.4 

1.81  E+00 

•6.4 

1.79  E+00 

-6.5 

1.81  E+00 

-6.4 

1.81  E+00 

-6.4 

1.80  E+00 

-6.0 

1.81  E+00 

■6.4 

1.81  E+00 

■6.4 

1.81  E+00 

-6.4 

Figure  3-33.  MTBE  Log  Transformation:  Simple  Uncertainty  Analysis  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Transect  Calculator:  Mass  Flux  Summary 


Site  Location  and  I.D.: 
Description: 


LFR  Site 
MTBE  and  TBA 


Mass  Flux  Summary  for  MTBE 


Distance  from 
Source  (ft) 

Mass  Flux  (g/day) 

Time  Period  1  Time  Period  2  Time  Period  3  Time  Period  4  Time  Period  5 

Transect  1 

Transect  2 

Transect  3 

Transect  4 
Transect  5 

1700 

"  1.70E+00 

85th  percentile  from  Advanced  Uncertainty  Analysis 
Maximum  of  the  three  Toolkit  interpolation  schemes 

Value  using  selected  interpolation  schemes 

Minimum  of  the  three  Toolkit  interpolation  schemes 
15th  percentile  from  Advanced  Uncertainty  Analysis 


I 


Next  Step: 
Save  File 


View 


- - -  ^ 

Flux  Result  I  View  Cone  Grid  |  Print  I  HELP 


Return  To  Toolkit  Main  Screen 


I 


Return  To  Transect  Main  Screen 


Figure  3-34.  MTBE  Log  Transformation:  Mass  Flux  Summary 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Concentration 

Transect  1  Time  Period  1.  Nearest  Neighbor  Interpolation 


Step  10:  Interpolate  Cone  Horizontally 

Export  Cone  Grid  1 

Import  Cone  Grid 

. 

Back  to  Data  Input 

j 

| 

|  Print  J 

[  HELP 

Next  Step: 

Interpolate  Hydraulic 
Conductivity  Grid 


1 


Start  of  Transect 
0.0 


7.1 

11.6 

16.0 

20.5 
24.9 

29.4 

33.8 

38.2 

42.7 

47.1 

51.6 
56.0 

60.5 

64.9 

69.3 

73.8 

78.2 
82.7 1 
87.1 

91.6 
96.0 


Data  Input  Instructions 


Enter  value  directly. 
(Values  in 
italics  represent 
interpolated  values.) 

Top/Bottom  of  plume. 
(Dont  enter  any 
data). 


TBA  Concentrations  (ug/L) 

Distance  from  Edge  of  Transect  (ft) 


MW  -1  MW-2  MW -3  MW -4  MW -5 

50.0  100.0  212.5  325.0  412.5  500.0  582.5  665.0  735.0  805.0 

0 

5.00E+00 

i.ooe+01  ^  -  - ... , .  ’  - 1 .  . . .. \  . . . . . 

0 

5.00E+00 

1.00E+01 

3.25E+01 

5.50E+01 

0 

5.00E+00 

1.00E+01 

3.25E+01 

5.50E+01 

3.25E+01 

1.00E+01 
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Figure  3-35.  TBA  Nearest  Neighbor  Interpolation:  Concentration  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 
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Figure  3-36.  TBA  Nearest  Neighbor  Interpolation:  Concentration  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Mass  Flux  Result 
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Data  Representation 

1 .  Bold  values  represent  calculations  based  on  given  values. 

2.  Values  in  italics  represent  calculations  based  on  interpolation. 

3.  Black  shaded  cells  represent  the  top  and  bottom  of  the  plume. 
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Figure  3-37.  TBA  Nearest  Neighbor  Interpolation:  Mass  Flux  Results  Grid  -  Page  1 
(For  conductivity  and  gradient  grids  see  Figures  3-6  through  3-9) 
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CASE  STUDY:  MTBE  AND  TBA 


o 


TBA  Mass  Flux  (g/day) 
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Figure  3-38.  TBA  Nearest  Neighbor  Interpolation:  Mass  Flux  Results  Grid  -  Page  2 
(For  conductivity  and  gradient  grids  see  Figures  3-6  through  3-9) 
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CASE  STUDY:  MTBE  AND  TBA 


Interpolation  Error  Results 


Site  Location  and  I.D.: 
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Figure  3-39.  TBA  Nearest  Neighbor  Interpolation:  Simple  Uncertainty  Analysis  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


_ 
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Figure  3-40.  TBA  Nearest  Neighbor  Interpolation:  Simple  Uncertainty  Analysis  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Transect  Calculator:  Mass  Flux  Summary 


Site  Location  and  I.D.: 
Description: 


LFR  Site 
MTBE  and  TBA 


Mass  Flux  Summary  for  TBA 

85th  percentile  from  Advanced  Uncertainty  Analysis 

Distance  from 
Source  (ft) 

Mass  Flux  (g/day) 

Time  Period  1  Time  Period  2  Time  Period  3  Time  Period  4 

Time  Period  5 

Maximum  of  the  three  Toolkit  interpolation  schemes 

Transect  1 

Transect  2 

Transect  3 

Transect  4 
Transect  5 

1700 

2.12E+01 

Value  using  selected  interpolation  schemes 

Minimum  of  the  three  Toolkit  interpolation  schemes 

15th  percentile  from  Advanced  Uncertainty  Analysis 

Distance  of  Transect  From  Source  (ft) 


Next  Step: 
Save  File 


I 


View  Flux  Result 


View  Cone  Grid 


).[  Print 


Return  To  Toolkit  Main  Screen  Return  To  Transect  Main  Screen 


Figure  3-41.  TBA  Nearest  Neighbor  Interpolation:  Mass  Flux  Summary 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Concentration 

Transect  1  Time  Period  1.  Vertical  Interpolation:  Linear  Horizontal  Interpolation:  Linear 


Step  10:  Interpolate  Cone  Horizontally 

Export  Cone  Grid  1 
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Figure  3-42.  TBA  Linear  Interpolation:  Concentration  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 
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Figure  3-43.  TBA  Linear  Interpolation:  Concentration  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Mass  Flux  Result 


TOTAL  MASS  FLUX 


(g/day) 


I  (kg/yr) 


Data  Representation 

1 .  Bold  values  represent  calculations  based  on  given  values. 

2.  Values  in  italics  represent  calculations  based  on  interpolation. 

3.  Black  shaded  cells  represent  the  top  and  bottom  of  the  plume. 


jNext  Step:  Mass  Flux 

Run/View  Uncertainty 
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Transect  1  ▼ 

Summary 

Analysis  (Optional) 

Back  to  Data  Grid  ||  Print  J  ^  HELP  j 

SELECT  TIME  PERIOD  TO  VIEW 
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Figure  3-44.  TBA  Linear  Interpolation:  Mass  Flux  Results  Grid  -  Page  1 
(For  conductivity  and  gradient  grids  see  Figures  3-15  through  3-18) 
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CASE  STUDY:  MTBE  AND  TBA 
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Figure  3-45.  TBA  Linear  Interpolation:  Mass  Flux  Results  Grid  -  Page  2 
(For  conductivity  and  gradient  grids  see  Figures  3-15  through  3-18) 
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CASE  STUDY:  MTBE  AND  TBA 


Interpolation  Error  Results 


Sits  Location  and  I.D.: 


LFR  Site 
MTBE  and  TBA 


Next  Step: 
Mass  Flux  Summary 


Back  to  Mass  Flux  Result 

Print 

Back  to  Data  Input 

HELP 
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nmuibm  mrd ir B 'll ty/Uvry  vrtonty  .vrl  ij.vlmt 
input  ititi. 

?.  Dunjtli-imil.Mrty  jnJy,!',,  initial  UirfinN  in  III- 

table  are  rancvea  ore  at  a  we  ana  re  mass  flur 
(.alLiialrO  wiffrjj*.  Vet  puwrt  using  theinte(d«tKjn 
scheme  specified  in  the  Data  Input  sector  (nose  that 


SELECT  TRANSECT  TO  VIEW  Transect  1  w 

SELECT  TIME  PERIOD  TO  VIEW  l  » 


TBA  Interpolation  Methods 

Hydraulic  Conductivity  f)  Vertical  t.inear  2)  1  tomontal  l  inear 
Concentration  1)  Vertical  linear  2)  Horrronlat  linear 
Hydraulic  Gradient  1)  Vertical  Linear  2)  Momenta!  Linear 

Total  Mass  Flux  Including  All  Points 


Total  Matt  Flux 

Paramatar 

Excluding  Salactad 

Contribution  of 

Top  of 

Bottom  of 

Value 

Point  and 

Salactad  Point  to 

Sampling 

Sampling 

Parameter 

Removed  For 

Paramatar 

Intarpolatlng 

Total  Matt  Flux 

End  of  Transact 

Interval 

Interval 

Examined 

Analytic 

Units 

(jldiy) 

W 

Sturt  u!  transect 

Conductivity 

1  60b+02 

ll/d 

2.08E+01 

-1.0 

MW-I 

60 

mg 

Conductivity 

1  60b  *02 

ft/d 

2.03E+01 

02 

MW  1 

230 

300 

Conductivity 

?  ?6F*0l 

ft/d 

2.10E+01 

-3.0 

MW-I 

51  0 

56  0 

Conductivity 

270001 

ft/d 

2.04E+01 

■0.3 

MW-I 

660 

710 

Conductivity 

i  ooc*oo 

ft/d 

2.09E+01 

-2.6 

MW-2 

7.0 

22  0 

Conductivity 

1  50b  *02 

U/d 

1.846+01 

9.7 

MW? 

40  n 

450 

Conductivity 

?6flF*0t 

ft/d 

280E+01 

-37.3 

MW? 

G?n 

(17  0 

Conductivity 

1  14F+0? 

fr/d 

1.05E+01 

48  6 

MW-2 

77  0 

02  0 

Conductivity 

1  OOE-OO 

ft/d 

2.33E+01 

-14.3 

MW-3 

10  u 

25  U 

Conductivity 

U  1  Ufc-01 

ft/d 

2.07E+01 

•1.7 

MW-3 

38.0 

43  U 

Conductivity 

8.8U6+Q1 

tt/d 

2.08E+01 

■22 

MW  3 

58  0 

63  0 

Conductivity 

i snF+n? 

ft/d 

1.82E+01 

10.8 

MW-3 

77.0 

U2  U 

Conducovity 

I.UUL  *00 

ft/d 

2.17E+01 

■6.7 

MW-4 

15  U 

300 

Conduct)  vi  tv 

2  1 0L  *01 

m 

2.06E+01 

•1.1 

MW-4 

44.0 

49  0 

Conductivity 

5  OOE+Ol 

ft/d 

2.08E+01 

-1.1 

MM  4 

(12  0 

68  0 

Conductivity 

8 70F+0I 

ft/d 

2.06E+01 

•1.0 

MW-4 

80.0 

050 

Conductivity 

1. OOL  *00 

It/d 

206E+O1 

•1.4 

MW-5 

20  0 

35  0 

Conductivity 

3  30L*01 

It/d 

2.06E+01 

•1.1 

MW-5 

47  0 

520 

•Conductivity 

5.70E-0, 

ft/d 

2.08E+01 

-1.1 

IWV.1 

650 

70  0 

Conductivity 

1 506+02 

ft/d 

2.08E+01 

-1.1 

MM', ! 

010 

860 

Conductivity 

t  00F*00 

ft/d 

2. 06 E+01 

-1.1 

MM-6 

24.0 

39.0 

Conductivity 

1  50L*02 

ll/d 

2.06E+01 

•1.1 

MM-C 

470 

520 

Conductivity 

5.5oe*01 

ft/d 

2.08  E+01 

-1.1 

MM-0 

62.0 

67  0 

Conductivity 

1  006+00 

ft/d 

2.06E+01 

-1.1 

MM  7 

300 

45  0 

Conductivity 

4  40001 

ft/d 

2 06 E+01 

-1.1 

MM -f 

520 

57.0 

Conductivity 

b  106+01 

ll/d 

2.06E+01 

-1.1 

MM-/ 

6/0 

720 

Conductivity 

1.006+00 

ll/d 

2.08  E+01 

-1.1 

MMQ 

29 

44 

Conductivity 

4.406-01 

ft/d 

2.06E+01 

-1.1 

MM-8 

51 

56 

Conductivity 

5106*01 

ft/d 

2.06E+01 

■1.1 

MM -fl 

66 

71 

Conductivity 

i  non  oo 

fVd 

2.06E+01 

•1.1 

End  otf  J  ransed 

Conductivity 

1  006+00 

It/d 

2.08E+01 

-1.1 

Start  of  Transect 

Concentration 

000  6-00 

,.y, 

2.07E+01 

-1.8 

MM  1 

5 

20 

Concentration 

1 OOF+OI 

ug/L 

2.06E+01 

•1.1 

MM-I 

33 

38 

Concentration 

1.006*0/ 

ug/L 

2.06E+01 

•1.1 

MM-1 

51 

56 

Concentration 

1  OOb+OI 

UgA. 

2.06E+01 

•1.1 

MM-1 

66 

71 

Concentration 

1  OOb+OI 

ug/L 

2.08E+01 

-1.1 

MM2 

7 

Concentration 

5  506+01 

ug/L 

2.57E+01 

-264 

MM-2 

40 

45 

Concentration 

3  I0C>02 

UQ/L 

2.35E+01 

-155 

MM-2 

62 

67 

Concentrator? 

0.906*02 

uq/L 

9.24  E+00 

64.6 

MM-2 

77 

82 

Concentration 

1.006+01 

ug/L 

2.38  E+01 

-18.1 

MM2 

10 

25 
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1  00F+01 

•&> 

2.08E+01 

•1.1 

MM2 

3fl 

43 

Conrentralnin 

1 OOF+Of 

ug/t 

2.08E+01 

•1.1 

MM-3 

58 

63 

Concentrator? 

1  00C>01 

ug/L 

2.06E+01 

-2.3 

MM-3 

77 

82 

Concentrator? 

3  006+01 

ug/L 

2.05E+01 

■0.7 

MM-4 

1b 

30 

Concentration 

1  out- 02 

ug/L 

2.08E+01 

-1.1 

MM  4 

44 

49 

Concentration 

1  OOF  02 

2  06E+01 

•1.1 

MM-4 

63 

68 

Concentration 

1  00C-O2 

ugA. 

2.06E+01 

•1.1 

MM-4 

80 

85 

Concentration 

1  OOC-02 

ug/L 

2.06E+01 

•1.1 

MM-b 

20 

3b 

Concentration 

1  006-02 

ug/L 

2.08  E+01 

-1.1 

MM  5 

47 

52 

Concentration 

1  OOF  02 

14,0 

2.08E+01 

-1.1 

MM  5 

65 

70 

Concentration 

1  OOF  02 

Uf}/I 

2.06E+01 

-1.1 

MM-b 

01 

86 

Concentration 

1. 006-02 

ug/L 

2.06E+01 

•1.1 

MM-6 

24 

39 

Concentration 

1  00E-O2 

ug/L 

2.06E+01 

•1.1 

MM -6 

47 

52 

Concentration 

1  OOE-02 

ug/L 

2.08E+01 

-1.1 

MM  f, 

62 

67 

Concentration 

t  OOF  02 

"90 

2.06E+01 

•1.1 

MM-/ 

30 

45 

Concentration 

1  006-02 

USM 

2.06E+01 

-1.1 

MM-/ 

52 

57 

Concentration 

1  006-02 

ug/L 

2.06E+01 

•1.1 

MM-7 

67 

72 

Concentration 

1  0OE-O2 

ug/L 

2.08E+O1 

-1.1 

MM  l\ 

29 

44 

Concentration 

1  OOF  02 

iig/L 

2.08E+01 

•1.1 

MM  8 

51 

56 

Concentration 

1  OOF  02 

>&} 

2.06  E+01 

-1.1 

MM-8 

66 

71 

Concentration 

1.006-02 

ug/L 

2.08  E+01 

-1.1 

End  of  Tranced 

Conccntraton 

0  !WE*0G 

ug/L 

2.06  E+01 

-1.1 
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Gradient 

7  I OE-04 

m 

2.06  E+01 
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MM  1 

5 

20 

Gradient 
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2.06E+01 

-1.1 
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33 

38 
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56 

Gradient 
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71 
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Figure  3-46.  TBA  Linear  Interpolation:  Simple  Uncertainty  Analysis  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


Parameter 


Top  of 
Sampling 


Bottom  of 


Value 


Sampling 


Parameter 


Removed  For 


Param  eter 


End  of  Transect 


Interval 


Interval 


Examined  A 


nalysis 


Units 


MW-2 
MW-2 
MW-2 
MW-3 
MW-3 
MW-3 
MW-3 
MW-4 
MW-4 
MW-4 
MW-4 
MW-5 
MW-5 
MW-5 
MW-5 
MW-6 
MW-6 
MW-6 
MW -7 
MW -7 
MW -7 
MW-8 
MW-8 
MW-8 

End  of  Transect 


40 

62 

77 

10 

38_ 

58 

77 

15 

44 

63 

80 

20 

47 

65 
81 
24 
47 
62 
30 
52 
67 
29 
51 

66 


45 

67 
82 
25 

43 
63 
82 
30 
49 

68 

85 
35 
52 

70 

86 
39 
52 
67 
45 
57 
72 

44 
56 

71 


Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 

Gradient 


2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
2.13E-03 
1.9  IE-03 
1.9 1E-0  3 
1.91E-03 
1.91E-03 
1.91E-03 
1.91E-03 
1.91E-03 
1.91E-03 
1.91E-03 
1.91E-03 


ft/ft 

m 

ft/ft 

m 

ft/ft 

fVft 

m 

m 

m 

m 

m 

ft/ft 

m 

m 

ft/n 

m 

ft/n 

m 

m 

m 

m 

m 

m 

m 

m 


RESULTS 


Total  Mass  Flux 
Excluding  Selected 
Point  and 
Interpolating 
(g/day) 

Contribution  of 

Selected  Point  to 

Total  Mass  Flux 

(%) 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

i  2.06E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06E+01 

-1.1 

2.06  E+01 

-1.1 

2.06E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

2.06  E+01 

-1.1 

Figure  3-47.  TBA  Linear  Interpolation:  Simple  Uncertainty  Analysis  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Transect  Calculator:  Mass  Flux  Summary  vmo 


Site  Location  and  I.D.: 
Description: 


LFR  Site 
MTBE  and  TBA 


Mass  Flux  Summary  for  TBA 

85th  percentile  from  Advanced  Uncertainty  Analysis 

Distance  from 
Source  (ft) 

Mass  Flux  (g/day) 

Time  Period  1  Time  Period  2  Time  Period  3  Time  Period  4 

Time  Period  5 

Maximum  of  the  three  Toolkit  interpolation  schemes 

Transect  1 
Transect  2 

Transect  3 

Transect  4 
Transect  5 

1700 

2.04E+01 

Value  using  selected  interpolation  schemes 

Minimum  of  the  three  Toolkit  interpolation  schemes 

15th  percentile  from  Advanced  Uncertainty  Analysis 

Figure  3-48.  TBA  Linear  Interpolation:  Mass  Flux  Summary 
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CASE  STUDY:  MTBE  AND  TBA 


Grid  Completion:  Concentration 

Transect  1  Time  Period  1.  Vertical  Interpolation:  Log  Horizontal  Interpolation:  Log 


I 


Step  10:  Interpolate  Cone  Horizontally 


Export  Cone  Grid 


Back  to  Data  Input 


Back  to  Cone  Grid 


Back  to  K  Grid 


I 


Import  Cone  Grid 


HELP 


Next  Step: 

Interpolate  Hydraulic 
Conductivity  Grid 


] 


Data  Input  Instructions 


Enter  value  directly. 
(Values  in 
italics  represent 
interpolated  values.) 

Top/Bottom  of  plume 
(Don't  enter  any 
data). 


TBA  Concentrations  (ug/L) 

Distance  from  Edge  of  Transect  (ft) 


Start  of  Transect 
0.0 

50.0 

MW  -1 

100.0 

212.5 

MW-2 

325.0 

412.5 

MW-3 

500.0 

582.5 

MW -4 

665.0 

735.0 

MW -5 

805.0 

7.1 

0 

3.16E-03 

1 .00E+01  *•  Tr-,. 

11.6 

0 

3.16E-03 

1.00E+01 

2.35E+01 

5.50E+01 

16.0 

0 

3.16E-03 

1.00E+01 

2.35E+01 

5.50E+01 

2.35E+01 

1.00  E+01 

20.5 

0 

3.16E-03 

1. 00E+01 

2.79E+01 

7.77E+01 

2.79E+01 

1. 00E+01 

24.9 

0 

3.16E-03 

1.00E+01 

3.31E+01 

1.10E+02 

3.31E+01 

1.00  E+01 

3.16E-01 

1.00E-02 

1.00E-02 

1.00E-02 

29.4 

0 

3.16E-03 

1.00E+01 

3.94E+01 

1.55E+02 

3.94E+01 

1.00  E+01 

3.16E-01 

1.00E-02 

1.00E-02 

1.00E-02 

33.8 

0 

3.16E-03 

1.00E+01 

4.68E+01 

2.19E+02 

4.68E+01 

1. 00E+01 

3.16E-01 

1.00E-02 

1.00E-02 

1.00E-02 

38.2 

0 

3.16E-03 

1.00E+01 

5.57E+01 

3.10E+02 

5.57E+01 

1. 00E+01 

3. 16E-01 

1.00E-02 

1.00E-02 

1.00E-02 

42.7 

0 

3.16E-03 

1.00E+01 

5.57E+01 

3.10E+02 

5.57E+01 

1. 00E+01 

3.16E-01 

1.00E-02 

1.00E-02 

1.00E-02 

47.1 

0 

3.16E-03 

1. 00E+01 

6.35E+01 

4.04E+02 

6.35E+01 

1.00  E+01 

3.16E-01 

1.00E-02 

1.00E-02 

1.00E-02 

51.6 

0 

3.16E-03 

1.00E+01 

7.25E+01 

5.25E+02 

7.25E+01 

1.00  E+01 

3.16E-01 

1.00E-02 

1.00E-02 

1.00E-02 

56.0 

0 

3.16E-03 

1.00E+01 

8.27E+01 

6.84E+02 

8.27E+01 

1.00  E+01 

3.16E-01 

1.00E-02 

1.00E-02 

1.00E-02 

60.5 

0 

3.16E-03 

1.00E+01 

9.43E+01 

8.90E+02 

9.43E+01 

1. 00E+01 

3. 16E-01 

1.00E-02 

1.00E-02 

1.00E-02 

64.9 

0 

3.16E-03 

1.00E+01 

5.38E+01 

2.90E+02 

6.18E+01 

1.32E+01 

3.63E-01 

1.00E-02 

1.00E-02 

1.00E-02 

69.3 

0 

3.16E-03 

1.00E+01 

3.07E+01 

9.43  E+01 

4.04E+01 

1.73E+01 

4.16E-01 

1.00E-02 

1.00E-02 

1.00E-02 

73.8 

0 

3.16E-03 

1.00  E+01 

1.75E+01 

3.07  E+01 

2.65E+01 

2.28E+01 

4.77E-01 

1.00E-02 

1.00E-02 

1.00E-02 

78.2 

0 

3.16E-03 

1.00E+01 

1. 00E+01 

1. 00E+01 

1. 73E+01 

3. 00  E+01 

5.48E-01 

1.00E-02 

1.00E-02 

1.00E-02 

82.7 

j  ooe+oi 

1. 00E+01 

1.73E+01 

3.00  E+01 

5.48E-01 

1.00E-02 

1.00E-02 

1.00E-02 

87.1 

'  . :  i  ooe+oi 

1. 00E+01 

1.73E+01 

3. 00  E+01 

5.48E-01 

1.00E-02 

1.00E-02 

1.00E-02 

91.6 

1.00E-02 

1.00E-02 

1.00E-02 

1.00E-02 

Figure  3-49.  TBA  Log  Transformation:  Concentration  Grid  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 


TBA  Concentrations  (ug/L) 

Distance  from  Edge  of  Transect  (ft) 

MW -6  MW -7  MW -8  End  of  Transect 

905.0  1005.0  1095.0  1185.0  1307.5  1430.0  1565.0  1700.0 


7.1 

11.6 

16.0 

20.5 
24.9 

29.4 

33.8 

38.2 

42.7 

47.1 

51.6 
56.0 

60.5 

64.9 

69.3 

73.8 

78.2 

82.7 
87.1 

91.6 
96.0 
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1.00E-02 


1.  OOE-02 
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1. 00E-02 


1. OOE-02 


1. 00E-02 
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1.00E-02 


1. OOE-02 


1.00E-02 


1. 00E-02 


1.00E-02 


1. 00E-02 
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1.00E-02 


1.00E-02 


1. 00E-02 


1.00E-02 


1. 00E-02 


1. OOE-02 


1.00E-04 


1.00E-02 


1.00E-02 


1.00E-02 


1.00E-02 


1. 00E-02 


1. 00E-02 
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1.00E-02 


1. 00E-02 


1. 00E-02 
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Figure  3-50.  TBA  Log  Transformation:  Concentration  Grid  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Mass  Flux  Result 


TOTAL  MASS  FLUX 


I  (g/day) 


I  (kg/yr) 


Data  Representation 

1.  Bold  values  represent  calculations  based  on  given  values. 

2.  Values  in  italics  represent  calculations  based  on  interpolation. 
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Figure  3-51.  TBA  Log  Transformation:  Mass  Flux  Results  Grid  -  Page  1 
(For  conductivity  and  gradient  grids  see  Figures  3-26  through  3-29) 
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CASE  STUDY:  MTBE  AND  TBA 
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Figure  3-52.  TBA  Log  Transformation:  Mass  Flux  Results  Grid  -  Page  2 
(For  conductivity  and  gradient  grids  see  Figures  3-26  through  3-29) 
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CASE  STUDY:  MTBE  AND  TBA 
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Figure  3-53.  TBA  Log  Transformation:  Simple  Uncertainty  Analysis  -  Page  1 
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CASE  STUDY:  MTBE  AND  TBA 
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Figure  3-54.  TBA  Log  Transformation:  Simple  Uncertainty  Analysis  -  Page  2 
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CASE  STUDY:  MTBE  AND  TBA 


Transect  Calculator:  Mass  Flux  Summary 
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Figure  3-55.  TBA  Log  Transformation:  Mass  Flux  Summary 
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